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Abstract 


Quantum hadrodynamics (QHD) is a framework for describing the nuclear 
many-body problem as a relativistic system of baryons and mesons. Motiva¬ 
tion is given for the utility of such an approach and for the importance of bas¬ 
ing it on a local, Lorentz-invariant lagrangian density. Calculations of nuclear 
matter and finite nuclei in both renormalizable and nonrenormalizable, effec¬ 
tive QHD models are discussed. Connections are made between the effective 
and renormalizable models, as well as between relativistic mean-field theory 
and more sophisticated treatments. Recent work in QHD involving nuclear 
structure, electroweak interactions in nuclei, relativistic transport theory, nu¬ 
clear matter under extreme conditions, and the evaluation of loop diagrams 
is reviewed. 
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1. INTRODUCTION 


The study of atomic nuclei plays an important role in the development of many-body the¬ 
ories. Early experimental probes of the nucleus were limited to energy scales considerably less 
than the nucleon mass M ^ 939MeV/c^, and the nucleus has traditionally been described 
as a collection of nonrelativistic nucleons interacting through an instantaneous two-body 
potential, with the dynamics given by the Schrodinger equation. The two-body potential is 
fitted to the empirical properties of the deuteron and to low-energy nucleon-nucleon (NN) 
scattering data, and one then attempts to predict the properties of many-nucleon systems. 
This is a difficult problem, because the NN potential is strong, short ranged {R ~ Ifm), 
spin dependent, and has a very repulsive central core. Nevertheless, over a period of many 
years and with the advent of more and more powerful computers, reliable methods have 
been developed for solving the nonrelativistic nuclear many-body problem. 

A new generation of accelerators will allow us to study nuclei at higher energies, at shorter 
distances, and with greater precision than ever before. For example, electron-nucleus scat¬ 
tering at CEBAF (now known as the Thomas Jefferson National Accelerator Facility) will 
sample distance scales down to tenths of a Fermi, and ultra-relativistic heavy-ion collisions 
at RHIC may produce nuclear densities of 10 times equilibrium density and temperatures 
of several hundred MeV. These experiments will clearly involve physics that goes beyond 
the Schrodinger equation, such as relativistic motion of the nucleons, dynamical meson 
exchanges, baryon resonances, modifications of hadron structure in the nucleus, and the dy¬ 
namics of the quantum vacuum, which will include the production of a quark-gluon plasma. 
The challenge to theorists is to develop techniques that describe this new physics, while 
maintaining the important general properties of quantum mechanics, Lorentz covariance, 
electromagnetic gauge invariance, and microscopic causality. 

Since quantum chromodynamics (QCD) of quarks and gluons is the fundamental theory 
of the strong interaction, it is natural to look to QCD as the means to describe this new 
physics. While this may be desirable in principle, there are many difficulties in practice, 
primarily because the QCD coupling is strong at distance scales relevant for the vast ma¬ 
jority of nuclear phenomena. Although significant progress has been made in performing 
strong-coupling lattice calculations, actual QCD predictions at nuclear length scales with 
the precision of existing (and anticipated) data are not presently available. This situation 
will probably persist for some time, particularly with regard to many-nucleon systems. 

In contrast, a description based on hadronic degrees of freedom is attractive for several 
reasons. First, these variables are the most efficient at normal densities and low tempera¬ 
tures, and for describing particle absorption and emission, because these are the degrees of 
freedom actually observed in experiments. Second, hadronic calculations can be calibrated 
by requiring that they reproduce empirical nuclear properties and scattering observables; we 
can then extrapolate to the extreme situations mentioned earlier. As an example, accurate 
microscopic meson-exchange models have been constructed to describe the NN interaction 
|[M a73| , pF7S| , pr’a79| , |Zu8 l| , |lVla^jTe87| , [lVla89| , piU^ , [St94|| . These contain several mesons, the most 
important of which are the 7r(0“,l), ct( 0’'',0), a;(l“,0), and p(l“,l), where the indicated 
quantum numbers denote spin, parity, and isospin, together with both N and A(1232) de¬ 
grees of freedom. Moreover, relativistic held theories of hadrons have been successful in 
describing the bulk and single-particle properties of nuclei and nuclear matter in the mean- 
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field and Dirac-Brueckner-Hartree-Fock approximations |^e86|JTb87|,[Ma89| , ^e92 


Our basic goal is therefore to formulate a consistent microscopic treatment of nuclear 
systems using hadronic (baryon and meson) degrees of freedom. In principle, one could 
derive the form of the hadronic theory from the underlying QCD lagrangian, using the ideas 
of modern Effective Field Theory ||We90a| , [We91| , [We92a| , |Ge93|| . For example, one could dehne 
the low-energy effective theory by requiring scattering amplitudes computed using hadrons 
to “match” the corresponding amplitudes derived from the underlying theory. Unfortu¬ 
nately, unlike certain cases where this matching can be done (the Standard Model of weak 
interactions, QED descriptions of atomic physics, and interactions in heavy-quark systems) 
[Ma96|| , the derivation of a low-energy hadronic theory directly from QCD is intractable at 
present. Thus we must rely on other properties of QCD to constrain hadronic lagrangians. 

One constraint is that the effective hadronic theory should embody the symmetries of 
QCD: Lorentz invariance, parity conservation, isospin symmetry, and spontaneously broken 
chiral symmetry.]^ These symmetry constraints severely restrict the forms of the allowed 
interactions, but they are insufficient to completely specify the low-energy dynamics. Thus 
we will use existing phenomenology to guide us in the construction of the effective lagrangian, 
in order to hnd the relevant degrees of freedom and the most efficient ways to structure the 
interactions. 

For example, meson-exchange models of the NN interaction tell us which mesons and 
baryons are the most important; these models have been useful both for scattering and for the 
nuclear matter problem. Moreover, experiments on the electrodisintegration of the deuteron 
show unambiguously the presence of pion-exchange currents, which arise when the incoming 
(virtual) photon couples to a pion being exchanged between two nucleons. Medium-energy 
pion-nucleus and proton-nucleus scattering indicate the importance of baryon resonances, 
such as the A, in nuclear reactions. In addition, a decomposition of the low-energy NN 
scattering amplitude using Lorentz invariants reveals that the empirical Lorentz scalar, vec¬ 
tor, and pseudoscalar amplitudes are much larger than the amplitudes deduced from a 
nonrelativistic decomposition containing Galilean invariants. These large amplitudes have 
important consequences for the spin, velocity, and density dependence of the NN interaction 
and are at the heart of relativistic descriptions of nucleon-nucleus scattering that reproduce 
spin observables in a very economical fashion ||Mc83|,ph83|,|C183|JWa87|| . Finally, relativis¬ 


tic mean-held models using classical Lorentz scalar and vector helds show that these are 
useful degrees of freedom for describing bulk and single-particle nuclear properties and for 
elucidating the important density dependence in the NN interaction. 

A Lorentz-covariant description is important for extrapolation to astrophysical objects 
and for describing processes at large energy-momentum transfer, as will be observed with the 
new accelerators. On the other hand, recasting the nuclear many-body problem in nonrela¬ 
tivistic form (he., with two-component rather than four-component nucleon spinors) leads to 
interactions (both two-body and few-body) that are similar to those used in successful cal¬ 
culations of few-nucleon systems. Many years of study within the nonrelativistic framework 


^The final two symmetries are only approximate. We also want to maintain the usual gauge 
invariance in electromagnetic interactions. 
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have produced a quantitative description of the structure of light nuclei (here both NN and 
NNN potentials are important) ||Wi92| , Ff% , Pu95 |. Moreover, we have a qualitative under¬ 
standing of some general features of nuclear structure, such as the importance of the Pauli 
principle in reducing NN correlations |[Wa95|| , which justihes both the shell model and the 
single-nucleon optical potential, and the interplay of single-particle and collective degrees 
of freedom, which determines the shape of the nucleus |[Ei87|| . Although the equilibrium 


properties of nuclear matter cannot be precisely reproduced with modern NN potentials 
Pa83| , pa85| , |Ja92^ , the addition of a phenomenological, density-dependent interaction to 


the free NN potential leads to excellent results for the single-particle structure and charge 
and mass densities of a large number of nuclei po79|| . One of the goals of the effective 
hadronic theory is to provide a deeper understanding of these successes and a more concrete 
link between these results and the underlying QCD. 

We desire a microscopic treatment of the nuclear many-body problem that is consistent 
with quantum mechanics, special relativity, unitarity, causality, cluster decomposition, and 
the intrinsic symmetries mentioned earlier. The modern viewpoint is that relativistic quan¬ 
tum held theory based on a local, Lorentz-invariant lagrangian density is simply the most 
general way to parametrize an S matrix (or other observables) consistent with these desired 
properties ||We95|| . Thus there is no reason that relativistic quantum held theory should be 
reserved for “elementary” particles only. We will refer to relativistic held theories based on 
hadrons as quantum hadrodynamics or QHD |^e86|JSe92|| . In principle, the held-theoretic 
formulation allows for the construction of “conserving approximations” ||Ba61|,pa62|| that 


maintain the general properties mentioned above. Analyses based on QHD, as dehned here, 
can hopefully provide a correct description of many-baryon systems at large distances and at 
energy scales that are not too high. One must examine, however, the limitations of hadronic 
held theory when one attempts to extrapolate away from this regime. 

Although spontaneously broken chiral symmetry implies that tttt and vrN interactions are 
weak for small momenta, the NN interaction is too strong to be treated perturbatively in the 
couplings. (This is also true in nonrelativistic formulations of the nuclear many-body prob¬ 
lem.) It is therefore necessary to develop reliable nonperturbative approximations, so that 
unambiguous comparisons between theory and experiment can be made. The formulation 
of practical, reliable techniques for hnite-density calculations in strong-coupling relativistic 
quantum held theories is still basically an unsolved problem |Pa86| , |Ka87| , |hu89| , [Ko95d| , |Ka96a| . 
The development of such tools in a hadronic held theory is not only useful in its own right, 
but it may also provide insight into similar approaches for QCD. Nuclear many-body theory 
has had such inhuence on other areas of physics in the past. 

Historically, QHD models ||Wa74| , pe79|| were conhned to the class of renormalizable held 
theories, which can be characterized by a hnite number of coupling constants and masses 
Ca82| , po84| . The motivation was that these parameters could be calibrated to observed 


nuclear properties, and one could then extrapolate into regions of extreme density, temper¬ 
ature, or energy-momentum transfer without the appearance of new, unknown parameters 
Wa74||. The hope was that nonrenormalizable and vacuum ehects would be small enough 


that they could be adequately described by the long-distance structure of a renormalizable 
hadronic theory, through the systematic evaluation of quantum loops |[Se92||. These models 


had several successes, some of which will be described below, but also difficulties in carrying 
out the renormalization program, as well as concrete indications that renormalizability is too 


4 






































restrictive. Thus it is important to generalize these models to the more modern viewpoint 
of nonrenormalizable, effective held theories, which can still provide a consistent treatment 
of the nuclear many-body problem. Reviewing the progress made toward this generalization 
in the past few years is a central theme of this article. 

Another goal of this article is to review the developments in QHD that have transpired 
since 1992.| We begin by discussing simple renormalizable models: QHD-I |[Wa74|| , which 
contains neutrons, protons, and the isoscalar, Lorentz scalar and vector mesons a and a;; 
the linear a model ||Sc57| , pe^JLe72|| , which contains neutrons, protons, pions, and neutral 
scalar mesons interacting in a chirally invariant fashion; and the extension of these models to 
include the isovector p meson, which we generically call QHD-II |[Se79|,[CTa82|| . These mod¬ 


els serve as pedagogical tools for introducing the relativistic mean-held and Dirac-Hartree 
approximations, and their application to both inhnite nuclear matter and the ground states 
of atomic nuclei. At the relativistic mean-held theory (RMFT) level with classical, isoscalar, 
Lorentz scalar and vector helds, QHD-I provides a simple picture of the equilibrium prop¬ 
erties of nuclear matter. We focus on the new features that arise in a relativistic framework 
and emphasize the important concepts of Lorentz covariance and self-consistency. When ex¬ 
tended to hnite nuclei through Dirac-Hartree calculations with a classical, isovector, Lorentz 
vector held and a few parameters htted to the properties of nuclear matter, one derives the 
nuclear shell model and can obtain, through the relativistic impulse approximation (RIA), 
an efficient description of the scattering of medium-energy nucleons from nuclei. We also 
show that the bulk and single-particle properties of nuclei provide stringent enough con¬ 
straints to distinguish between various RMFT models. These successful QHD results are 
discussed in Section 

In Section we examine the linear sigma model. We show how the constraints of linear 
chiral symmetry and spontaneous symmetry breaking, when realized in the usual fashion, 
preclude a successful description of nuclei at the mean-held level. We then discuss how the 
exchange of two correlated pions between nucleons produces a strong mid-range attraction in 
the scalar-isoscalar channel; this allows for an alternative (nonlinear) implementation of the 
chiral symmetry, with an additional scalar held to simulate correlated two-pion exchange, 
that produces successful mean-held results. We conclude that the scalar held in QHD- 
I should be regarded as an effective held that includes the pion dynamics that is most 
important for describing bulk nuclear properties. 

In Section we consider the application of ehective hadronic theories to the two- and 
few-nucleon problems. Several meson-exchange models of the NN interaction have been 
developed and applied to nuclear systems, both relativistically and nonrelativistically. These 
models focus primarily on the precise reproduction of the two-nucleon data, and they are 
typically not renormalizable. We discuss the relativistic quasipotential formalism and its 
successful application to two-nucleon systems, as well as its extension to the study of two- 
nucleon correlations in nuclear matter. We then examine how the RMFT results are modihed 
when short-range correlations are included and also gain some insight into the relationship 
between relativistic and nonrelativistic approaches to the nuclear matter problem. The 


^We try to be as complete as possible with references from early 1992 through 1995. 
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successes of these nonrenormalizable models, as well as the introduction of an effective 
scalar held in Section give strong empirical hints that one must really view QHD as arising 
from an effective, nonrenormalizable lagrangian that yields an economical description of the 
underlying theory of QCD in the nuclear domain. 

Therefore, in Section we present a formulation of QHD in the modern context of effec¬ 
tive held theory. Since a nonrenormalizable QHD lagrangian generally contains an inhnite 
number of unknown constants, one must identify suitable expansion parameters and dehne 
a meaningful truncation scheme. Based on the results in Section we argue that although 
the mean meson helds are large on the scale of nuclear energies (several hundred MeV versus 
tens of MeV), and one has all the complexities of a strong-coupling theory, the ratios of the 
mean helds to the nucleon mass are small at moderate densities. These ratios thus provide 
suitable expansion parameters. This observation leads naturally into Section ^ where a wide 
range of applications of RMFT to nuclear structure throughout the periodic table and to a 
variety of nuclear phenomena is discussed. The numerous successful calculations since early 
1992 are considered within the general context of ehective held theory, and the importance 
of both renormalizable and nonrenormalizable couplings is considered. We also introduce 
some basic ideas of density functional theory ||Dr90|| to illustrate how RMFT models can in¬ 
clude many-body ehects that go beyond the simple Hartree approximation. When combined 
with the systematics of two-nucleon correlations described in Section ^ these ideas help to 
explain why RMFT models are so successful. 

The RMFT has also been widely applied to relativistic transport theory and heavy-ion 
reactions through the Vlasov-Uehling-Uhlenbeck equations for the phase-space distribution 
function. We briehy review the basic theoretical background and discuss developments 
since 1992 in Section |^. We also consider the extrapolation of RMFT results for the nuclear 
equation of state to high baryon density and high temperature. 

To have a consistent effective held theory, one must also consider loop corrections to 
the RMFT. These loops contain both familiar many-body contributions |[Fe71|| and quan¬ 
tum vacuum effects. The underlying vacuum dynamics of QCD is implicitly contained in 
the parameters of the effective hadronic lagrangian; nevertheless, these vacuum effects are 
modihed in the presence of valence nucleons at hnite density. Although chiral perturbation 
theory ||Ga84|,^e93a|| provides a systematic framework for including pion loops in meson- 
meson and meson-nucleon scattering amplitudes, the inclusion of vacuum loops that arise 
in the many-body theory is still an unsolved problem. We review the present status of loop 
calculations in QHD in Section Finally, Section contains a summary. 

The reader is assumed to have a working knowledge of relativistic quantum held the¬ 
ory, canonical quantization, and the use of path integrals at zero temperature. For gen¬ 
eral background on relativistic held theory, the reader can turn to a number of texts 
Bj64IB,i65| , |Na78| , r^JRa81| , pVla84|JRy85|JKa89|JSt93| , |Um93| , |Gr93| , |We95| , |We9^ ; for background 


on quantum hadrodynamics, two recent reviews exist that cover the basic formulation and 


developments up to 1992 ||Se86|,|Se92 
in more detail ||Wa95 


A recent text develops many of the theoretical tools 


6 

































2. QUANTUM HADRODYNAMICS (QHD) 


The purpose of these next three sections is to develop some insight into the ingredients 
necessary for a realistic relativistic description of few- and many-body nuclear systems. In 
this section, we review the basic ideas of QHD and some of its applications. Most of the 
concepts are developed in |^e86|JSe9^ , which discuss the literature up to 1992, and a recent 
text provides a more thorough treatment of the theoretical tools and background ||Wa95 


We will not repeat all of that material here. Moreover, since we are concentrating on basic 
results, we postpone until later the discussion of recent progress, where we will try to be as 
complete as possible with developments and references from early 1992 through 1995. While 
the literature in this area prior to 1992 is detailed in |^e86| , |Se9'^ , we would like to single out 
for special mention the important background papers ||Sc51| , |Jo53| , pu56|jWe67| , |Mi72|| . 

Quantum hadrodynamics as dehned above is a general framework for the relativistic nu¬ 
clear many-body problem. The detailed dynamics must be specified by choosing a particular 
lagrangian density. As we will discuss later, there is increasing evidence that QHD must be 
regarded as an effective held theory, where all types of hadronic couplings satisfying general 
symmetry requirements are to be included. Nonetheless, the original motivation for using 
a renormalizable lagrangian to formulate a consistent relativistic nuclear many-body theory 
based on hadronic degrees of freedom remains valid. Furthermore, as we shall see, there is 
now convincing empirical evidence that whatever the form of the effective theory, it must 
be dominated by strong isoscalar, Lorentz scalar and vector interactions. 

To introduce the relativistic formalism, we consider a simple model called QHD-I |[Wa-74|| , 

'V'p' 

The lagrangian density for this model (with h = c = 1) is given by Ubet 

- g^V^) - (M - - ml(j)b 

- + 6C , (2.1) 

where and 5C contains counterterms. The parameters M, g^, 

rus, rriv, k, and A are phenomenological constants that may be determined (in principle) 
from experimental observables. This lagrangian resembles massive QED with an additional 
scalar interaction, so the resulting relativistic quantum field theory is renormalizable |[Bo70 . 
The inclusion of the scalar self-interactions proportional to 0^ and 0^ make this the most 


which contains helds for baryons 0 = 


and neutral scalar (0) and vector {V^) 


mesons. 


general lagrangian consistent with renormalizability (for these degrees of freedom) |[Bo77 . 
The counterterms in SC are used for renormalization. 

The motivation for this model has evolved considerably since it was introduced. As 
discussed in the Introduction, when the empirical NN scattering amplitude is described 
in a Lorentz-covariant fashion, it contains large isoscalar, scalar and four-vector pieces 
|[Mc83|,|STi83|,|C183|,|Wa87|| , and the simplest way to reproduce these is through the exchange 


of neutral scalar and vector mesons. The neutral scalar and vector components are the most 
important for describing bulk nuclear properties, which is our main concern here. Other 
Lorentz components of the NN interaction average essentially to zero in spin-saturated nu¬ 
clear matter and may be incorporated as refinements to the present model; we will discuss 
pion dynamics in the following section. The important point is that even in more com¬ 
plete models, the dynamics generated by scalar and vector mesons will remain; thus, it is 


7 


































important to first understand the consequences of these degrees of freedom for relativistic 
descriptions of nuclear systems. 

The held equations for this model follow from the Euler-Lagrange equations and can be 
written as 


( 9 „ 9 '‘ + = g.irip , 

+ mlv*‘ — 5vV'7'‘^ ) 
l7'‘(®Sp - 9.U) - (V - = 0 . 


( 2 , 2 ) 

(2.3) 

(2.4) 

(The counterterms have been suppressed.) Equation (|2.2|) is a Klein-Gordon equation with 
a scalar source term and nonlinear scalar self-interactions. Equation (|2.3D looks like massive 
QED with the conserved baryon current 




= 0 


(2.5) 


rather than the (conserved) electromagnetic current as the source. Finally, Eq. ( |2.4|) is 
the Dirac equation with scalar and vector helds entering in a minimal fashion. These held 
equations imply that the canonical energy-momentum tensor is conserved = 

= 0 ). 

When quantized, Eqs. (^(^) - (^(^) become nonlinear quantum field equations, whose exact 
solutions (if they exist) are very complicated. In particular, they describe mesons and 
baryons that are not point particles, but rather objects with structure due to the implied 
(virtual) meson and baryon-antibaryon loops. Here the dynamical input of renormalizability 
is apparent, since we are assuming that the intrinsic structure (or at least the long-range 
part of it) can be described using hadronic degrees of freedom. Strictly speaking, the validity 
of this input and its limitations have yet to be tested conclusively within the framework of 
QHD; nevertheless, as we shall see, there are now strong indications that this assumption is 
too optimistic. 

We also expect the coupling constants in Eqs. 


- (|2.4|) to be large, so perturbative 
solutions are not useful. Fortunately, there is an approximate nonperturbative solution that 
can serve as a starting point for studying the implications of the lagrangian in Eq. ( 0 ). 
Consider a system of B baryons in a large box of volume V at zero temperature. Assume 
that we are in the rest frame of the matter, so that the baryon flux B = 0. As the baryon 
density BjV increases, so do the source terms on the right-hand sides of Eqs. (E3) and 
( p.3|) . If the sources are large enough, the meson held operators can be approximated by 
their expectation values, which are classical helds: 


(.#.) = ^ 0 , = 


( 2 , 6 ) 


For our stationary, uniform system, 0o and Vq are constants that are independent of space 
and time, and since the matter is at rest, the classical three-vector held V = 0. 

It is important to emphasize that the preceding “mean-held” theory (MET) serves only 
as a starting point for calculating corrections within the framework of QHD, using Feynman 


diagrams and path-integral methods, as discussed in |pe86| , |Se9^[Wa95|] . We will return later 


to decide at which densities this starting point is actually useful. 






A. The Nuclear Matter Equation of State 


When the meson helds in Eq. (p^.l|) are approximated by the constant classical helds of 
Eq. (|2.6|) , we arrive at the mean-held lagrangian density 

- g^Vo'^o - (M - 5's0o)]^/’ - . (2.7) 


(The counterterms have been suppressed.) The conserved baryon four-current remains as in 
Eq. (|2.5|) , and the canonical energy-momentum tensor becomes 

^MFT = ■ (2-8) 


As discussed by Freedman |[Fr78||, there is no need to symmetrize if we consider only 


uniform nuclear matter. This follows because the additional terms in the symmetrized tensor 
enter as a total four-divergence, whose diagonal matrix elements vanish in a uniform system. 

Since the meson helds are classical, only the fermion held must be quantized. The Dirac 
held equation follows from £mft: 


- SvToVi - (M - x) = 0 , 


(2.9) 


and since this equation is linear^ it can be solved exactly. The scalar held 0o shifts the baryon 
mass from M to M* = M — gs4>o, while the vector held Vq shifts the energy spectrum. We 
look for normal-mode solutions with both positive and negative energies, as is natural for 
the Dirac equation. These solutions can be used to dehne quantum held operators -0 and 
-01 in the usual fashion, and by imposing the familiar equal-time anticommutation relations, 
we can construct the baryon number operator B = J d^x (: :) and the four-momentum 

operators = {H, P) = /d^xT°^, with the results 


H - (0|.ff|0) = .^mft + 5H , 

Apt = E(k" + + At Aa) + JvVoB 

kA 

+ , 

6H = -J2 [(k^ + M*2) 1/2 _ (1^2 + m2)1/2] , 

kA 

B = EC^kA-'lkA - AaAa) . 

kA 

P = i: kC^kA^tA + Aa Aa) ■ 

kA 


( 2 . 10 ) 

( 2 . 11 ) 

( 2 . 12 ) 

(2.13) 

(2.14) 


Here Aa. and are creation and destruction operators for (quasi)baryons and 
(quasi)antibaryons with shifted mass and energy, and B is the “normal-ordered” baryon 
number operator, which clearly counts the number of baryons minus the number of an¬ 
tibaryons. (The index A denotes both spin and isospin projections.) The correction term 
5H arises from placing the operators in Hmft in normal order and includes the contribution 
to the energy from the hlled Dirac sea, where the baryon mass has been shifted by the 
uniform scalar field 0o |[8e8(j||. Since all energies are measured relative to the vacuum, we 
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must subtract the total energy of the Dirac sea in the vacuum state |0), where the baryons 
have their free mass M. We will return later to discuss this “zero-point energy” correction; 
for now, let us concentrate on the MFT hamiltonian dehned by Eq. ( p.ll|) . 

Since Hmft is diagonal, this model mean-held problem has been solved exactly once the 
meson helds are specihed; their determination is discussed below. The solution retains the 
essential features of QHD; explicit mesonic degrees of freedom, consistency with relativistic 
covariance ||Eu9(J|| , and the incorporation of antiparticles. Since B and P are also diag¬ 
onal, the baryon number and total momentum are constants of the motion, as are their 
corresponding densities and V, since the volume is hxed. 

For uniform nuclear matter, the ground state is obtained by Filing energy levels with 
spin-isospin degeneracy 7 up to the Fermi momentum k^. (The generalization to hnite 
temperature will be discussed at the end of this section.) The Fermi momentum is related 
to the baryon density pe = B/V by 


Pb — 


7 


(2ir) 


p/cp 


|3l _ 7 73 


= 


k-^ 


(2.15) 


where the degeneracy factor is 4 for symmetric {N = Z) matter and 2 for pure neutron 
matter {Z = 0). The constant vector held Vq can be expressed in terms of conserved 


quantities from the expectation value of the vector meson held equation (|2.3|): 


En = 




mt 


tPb ■ 


(2.16) 


The expressions for the energy density and pressure now take the simple forms |^e 86 

,2 .. ^ 




*44 


(M - M*) 


+ 


7 


f'k-p 


(27r)3 Jo 
2 


d^k E*{k) , 


A 


p = -^pl - ^(M - M*f - Z—iM - M*f - 
^ 2mf^ 2gf > dgf ^ 24p, 


{M -M 


*4 4 
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3 (27r)3 Jo 


d\ 


E*{k) ’ 


(2.17) 


(2.18) 


where E*{k) = (k^ -I- The hrst four terms in Eqs. (p.l7| ) and ( p.l8| ) arise from the 

classical meson helds. The hnal terms in these equations are those of a relativistic gas of 
baryons of mass M*. These expressions give the nuclear matter equation of state at zero 
temperature in parametric form: T(pb) and p(Pb)- 

The constant scalar held (pQ, or equivalently, the ehective mass M*, can be determined 
thermodynamically at the end of the calculation by minimizing £{M*) with respect to M*. 
This produces the self-consistency condition 
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M* = M- ^p, + 


K 


mt 


2g^m. 


{M -M*Y + 


A 


{M -M*f , 


dg^ml 


(2.19) 
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FIG. 1. Saturation curves for nuclear matter. These results are calculated in the relativistic 
mean-field theory with baryons and neutral scalar and vector mesons (QHD-I). The coupling 
constants are chosen to fit the value and position of the minimum. The solid curve uses the 
couplings in Eq. ( 2.21 ), while the dashed curve uses the parameter set NLC, as described later in 
the text. The predictions for neutron matter (7 = 2) are also shown. 


where the scalar density ps is dehned by 


Ps = (: ijip :) = 


7 


C^F 


(2vr) 


(Fk 


M* 


E*{k) ■ 


( 2 . 20 ) 


Equation (|2.19|) is equivalent to the MET scalar held equation for 0o- Note that the scalar 
density is smaller than the baryon density [Eq. ( p.lSj) ] due to the factor M*/E*{k), which 
is an effect of Lorentz contraction. Thus the contribution of rapidly moving baryons to 
the scalar source is signihcantly reduced. Most importantly, Eq. (|2.19| ) is a transcendental 
self-consistency equation for M* that must be solved at each value of kp. This illustrates 
the nonperturbative nature of the mean-held solution. 

To analyze these results, we initially set k = A = 0, as in the original version of the model 


Wa74|| . An examination of the analytic expression (p.l7|) for the energy density shows that 


the system is unbound {£ f Pb > M) at either very low or very high densities. At intermediate 
densities, the attractive scalar interaction will dominate if the coupling constants are chosen 
properly. The system then saturates. Nuclear matter with an equilibrium Fermi wavenumber 
and an energy/nucleon cq = {S/p^ — M) = —15.75 MeV is obtained if the 


fcp = 1.30 fm 


-1 
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FIG. 2. Effective mass as a function of density for nuclear (7 = 4) and 
based on Fig. |^. 


neutron (7 = 2 ) matter 


couplings are chosen asQ 


\ rn.^ / 


m. 


= 357.4, 




mt, 


= 273.8 


( 2 . 21 ) 


In this approximation, the nuclear compression modulus K is 545 MeV. Note that the meson 
masses enter only through the ratios gf/mj in Eqs. (p.l7|) , (|2.18|) , and (p.l9|) . (For a more 
complete discussion of the relevant dimensional coupling parameters and their specification 
from nuclear matter properties, see ||Fu96|| .) The resulting saturation curve is shown in 


Fig.i In this approximation, the relativistic properties of the scalar and vector fields are 
responsible for saturation; a Hartree-Fock variational estimate built on the nonrelativistic 
(Yukawa) potential limit of the interaction shows that such a system is unstable against 
collapse |[Fe71|,[Wa95|] . 

The solution of the self-consistency condition (|2.19|) for M* yields an effective mass 
that is a decreasing function of the density, as illustrated in Fig. |^. Evidently, M*/M 
becomes small at high density and is significantly less than unity at ordinary nuclear densities 
{M*/M = 0.541 at = kp). This is a consequence of the large scalar field which is 
approximately 400 MeV (at kp = kp) and which produces a large attractive contribution 


^The values = 267.1 and 


195.9 used in 


| Se 86 | yield A:)) = 1.42fm 
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FIG. 3. Predicted equation of state for neutron matter at all densities. The solid and dashed 
curves show the result for QHD-I based on Fig. ||. The dotted line represents the causal limit 
p = £. The density regime relevant for neutron stars is also shown. 

to the energy/baryon. There is also a large repulsive energy/baryon from the vector field 
gyVo ~ 350 MeV. Thus the Lorentz structure of the interaction leads to a new energy scale in 
the problem, and the small nuclear binding energy 16 MeV) arises from the cancellation 
between the large scalar attraction and vector repulsion. As the nuclear density increases, 
M* decreases, the scalar source ps becomes smaller than the vector source pe, and the 
attractive forces saturate, producing the minimum in the binding curve.^ Clearly, because 
of the sensitive cancellation involved near the equilibrium density, corrections to the MET 
must be calculated before the importance of this saturation mechanism can be assessed. 
Nevertheless, the Lorentz structure of the interaction provides a new saturation mechanism 
that is not present in the nonrelativistic potential limit, as this limit ignores the distinction 
between Lorentz scalar and vector fields. We will see later that part of this new saturation 
mechanism can be expressed in terms of repulsive many-body forces in a nonrelativistic 
formulation. 

The corresponding curves for neutron matter obtained by setting 7 = 2 are also shown 
in Figs. Ill and||, and the equation of state (pressure vs. energy density) for neutron matter 
at all densities is given in Fig. At high densities, the system approaches the “causal 
limit” p = £ representing the stiffest possible equation of state. Thus we have a simple, 
two-parameter model that is consistent with the equilibrium point of normal nuclear matter 


^Figure 15 in Se8ti | shows the decrease in relative to for 0 < fep < /Cp. 
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and that allows for a covariant, cansal extrapolation to any density. 

Nevertheless, the two-parameter (linear) model is £t only to the eqnilibrinm point of 
nnclear matter; the valnes of M* and K at eqnilibrinm are predictions. For example, for a 
given Fermi wavennmber and energy/nncleon cq at eqnilibrinm, M* at eqnilibrinm must 
satisfy ||Fu96|| 


(e„ + M+ - (M - M-)P° - 7^ E-(k) = 


n 


( 2 . 22 ) 


As will be shown in the next subsection, the properties of hnite nuclei place signihcant 
constraints on both M* and K, and the linear model predicts too small a value for the 
former and too large a value for the latter. Moreover, the absence of isovector mesons in 
QHD-I leads to a bulk symmetry energy that is too small, and consequently, the repulsive 
forces in neutron matter are underestimated. This can be corrected by introducing a mean 
held for the isovector p meson, as we describe shortly. Thus, for comparison, we also show 
in Figs. |l] through the more realistic MFT results obtained in a nonlinear model (NLC), 
with parameters fit to the properties of nuclei and given in Table below. The important 
point is that the additional parameters allow for small adjustments of the nuclear properties 
near equilibrium, but the basic features implied by the large scalar and vector fields remain 
intact ||Wa74|J5e^ , Po89| , |Bo91|JFu96|| . 


B. Finite Nuclei 


We now generalize the results of the preceding subsection to study atomic nuclei. We 
continue to work in the mean-held approximation to QHD-I, but since the system now has 
hnite spatial extent, these helds are spatially dependent .Q If we initially restrict consideration 
to spherically symmetric nuclei, the meson helds depend only on the radius, and since the 
baryon current is conserved, the spatial part of the vector held V again vanishes |Se86 
Thus the mean-held QHD-I lagrangian of Eq. ( p.7|) becomes 


■CSU = - SvToVo - (M - - IKVi^io)^ + mJiAi 

- + mlvS] . 

and the Dirac equation for the baryon held is 


(2.23) 


- g^joVoir) - [M - gMr)]}i^ix) = 0 . 


(2.24) 


Appropriate values for the scalar and vector couplings {ps and g^), masses (mg and m^), and 
nonlinear parameters {k and A) will be given below. 

Although the baryon held is still an operator, the meson helds are classical; hence 
Eq. (|2.24|) is linear, and we may again seek normal-mode solutions of the form 'if{x) = 
'(/'(x) e“*®*. This leads to the eigenvalue equation 


®As usual in discussions of nuclear structure, calculations will be carried out in a frame where the 
nucleus is at rest. 
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/iV’(x) = {-ioL • V + 5 'vK)(r) + (3[M - gs(j)o{r)]}^p{x) = E'ijj{x) , 


(2.25) 


which dehnes the single-particle Dirac hamiltonian h, with a and (3 the usual Dirac matrices. 
Equation (|2.25| ) has both positive- and negative-energy solutions W(x) and V(x), which allow 
the held operators to be constructed in the Schrodinger picture. The positive-energy spinors 
can be written as 


Tjl \ — ij I \ _ ( km \ 

(x) = ^nkmt (x) — I r f \ / ICt) 

\-[Fnkt{r)/r\^-kmJ 


(2.26) 


where n is the principal quantum number, ^km is a spin-1/2 spherical harmonic |Pjld57|| , and 
Q is a two-component isospinor labeled by the isospin projection t. (We take t = | for 
protons and t = — | for neutrons.) The phase choice in Eq. (|2.26|) produces real bound-state 
wave functions F and G for real potentials 0o and Vq, and the normalization is given by 


dr(|G„(r)p + 


Fa{r)\^) = l. 


(2.27) 


which ensures unit probability to hnd each nucleon somewhere in space. 

The classical meson held equations follow from Eq. (p.23|) and resemble Eqs. (|2.2| ) and 
( P^.3| ) restricted to static, spherically symmetric helds. With the general form for the spinors 
in Eq. (p.26|) , we can evaluate the nuclear densities, which serve as source terms in the meson 


held equations. Assume that the nuclear ground state consists of hlled shells up to some value 
of n and fc, which may be diherent for protons and neutrons; this is appropriate for doubly 
magic nuclei. In addition, assume that all bilinear products of baryon operators are normal 
ordered, which removes contributions from the negative-energy spinors Vq,(x). This amounts 
to neglecting the hlled Dirac sea of baryons and dehnes the mean-held approximation. The 
contributions from the Dirac sea to nuclear matter will be considered in Section ^ 

With these assumptions, the local baryon (pe) and scalar (ps) densities become 


Pb(x) 

Ps(x) 




7 




X = 


occ 

E 


'2 ja + 1 


G„(r)p±|F„(r)p) 


(2.28) 


which holds for hlled shells, as appropriate for spherically symmetric nuclei. The remaining 
quantum numbers are denoted by {a} = {a;m} = {n,and the nonzero integer k 
determines j and ^ through k = (2j -|- l)(f' — j)- Notice that since the shells are hlled, the 
sources are spherically symmetric. 

The sources produce the meson helds, which satisfy static Klein-Gordon equations: 


dr^ 


0 o(?^) + 


2 d 
r dr 








dr^ 


Vo{r) 



mlVo{r) 


-9sPs(,r) , 
-9vpB{r) . 


(2.29) 

(2.30) 


The equations for the baryon wave functions follow upon substituting Eq. ( p.26|) into 
Eq. (p^.25[), which produces 
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(2.31) 


^Ga{r) + -Ga{r) - [E^ - (?vVo(r) + M - gs<l>o{'^)]Ea{r) = 0 , 
dr r 

- ^Fa{r) + [Ea - gvVo{r) - M + gs(j)o{r)]Ga{r) = 0 . (2.32) 


Thus the spherical nuclear ground state is described by coupled, ordinary differential equa- 

|M| , |Fu8l] 


tions that may be solved by an iterative procedure, as discussed 


m 


They 


contain all information about the static ground-state nucleus in this approximation. 

The mean-field hamiltonian can be computed just as for infinite matter, and after normal 
ordering, the ground-state energy is given by 



+ ^ Wl(x)[ - icx-V + (3{M- + g.Vo] W„(x)} . (2.33) 


Here the meson fields are functions of the radial coordinate. Notice that if we interpret this 
expression as an energy functional for the Dirac-Hartree ground state, extremization with 
respect to the meson fields reproduces the field equations (|2.29|) and ( |2.3CI| ), with the densities 
from Eq. (|2.28|) . Moreover, extremization with respect to the baryon wave functions Wj,(x), 
subject to the constraint 


d^x Wj((x)W„(x) = 1 


(2.34) 


for all occupied states (which is enforced by Lagrange multipliers Ea), leads to the Dirac 
equation (|2.25|) . This alternative derivation of the Dirac-Hartree equations from an energy 
functional is useful for extensions of the simple model discussed above ||Eu96 . 


Once the solutions to the Dirac-Hartree equations have been found, the ground-state 
energy can be computed by using the Dirac equation ( p.25|) to introduce the eigenvalues Ea 
and by partially integrating the meson terms to introduce the densities. In the end, the 
total energy of the system is given by 

E = ^Ea{2ja + 1) - 2 y d^a;[-^s0o(r)ps(r) + gyVo{r)p^{r)] - d^x[K(j)l{r) + |A0^(r)] . 

(2.35) 


Before discussing the Dirac-Hartree solutions, let us generalize the equations to include 
some additional degrees of freedom and couplings. Although the isoscalar meson fields 
are the most important for describing general properties of nuclear matter, a quantitative 
comparison with actual nuclei requires the introduction of some additional dynamics. 

For example, it is necessary to include the electromagnetic interaction to account for the 
Coulomb repulsion between protons. Moreover, since hadronic interactions exhibit an almost 
exact SU(2) isospin symmetry, the nucleons can couple to isovector mesons in addition to 
the isoscalar (neutral) mesons of QHD-I. These isovector mesons, for example the p and 
TT, come in three charge states (-(-, 0, —) and couple differently to the proton and neutron. 
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TABLE I. Dirac-Hartree Parameter Sets. Note that mg and k are in MeV. 


Set 

al 

al 

9p 

rUg 

K 

A 

L2 

109.63 

190.43 

65.23 

520. 

0 

0 

NLB 

94.01 

158.48 

73.00 

510.0 

800 

10 

NLC 

95.11 

148.93 

74.99 

500.8 

5000 

-200 


Thus they affect the nuclear symmetry energy, which arises when there are unequal numbers 
of neutrons and protons. 

The construction of a renormalizable lagrangian containing charged, massive vector helds 
is somewhat complicated and is discussed at length in Abers and Lee [|Ab73|| ; applications 
to the present model can be found in ||Se 86 || . For our purposes, we require only the classical 
contributions from these helds, and in this case, the lagrangian simplihes considerably. In 
particular, since the nuclear ground state has well-dehned charge, only the neutral rho meson 
held (denoted by ho) enters, and since the ground state is assumed to have well-dehned parity 
and spherical symmetry, there is no classical pion held. Thus the mean-held lagrangian for 
this extended model, which we call QHD-II, is given by 


-^MPT = - QvloVo - 5 'p|r 37 o&o - e|(l + r 3 ) 7 oAo - (M - 

- ^[(Wo)^ + + mlV^] 

+ • 


(2.36) 


Here Aq is the Coulomb potential, e is the proton electromagnetic charge, Qp is the rho- 
nucleon coupling constant, and r* are the usual isospin Pauli matrices. For now, all of the 
boson helds are assumed to be functions of the radial coordinate only. 

The Dirac-Hartree equations for this extended model can be derived just as before. The 
Dirac equations for the baryon wave functions now contain and Aq, and because of the 
structure of the matrix, Bq couples with opposite sign to protons and neutrons, and Aq 
couples only to the protons. In addition to the source terms in Eq. ( p.28|) , which sum 
over both proton and neutron occupied states, the source term for the p meson involves 
the difference between proton and neutron densities, while the Coulomb source involves 
only protons. These diherent types of couplings allow for a more accurate reproduction of 
real nuclei, where the proton and neutron wave functions are not identical. The full set of 
equations are presented in ||l 8 e 86 || and are used to compute the results discussed below. 


1. Spherical Nuclei 

The solutions of the preceding equations depend on the parameters Ps, mg, and 
Pp (when the p meson is included); n and A will be set to zero in this subsection. We 
take the experimental values M = 939 MeV, = 783 MeV, rup = 770 MeV, and 

e^ = a = 1/137.036 (which determines the Coulomb potential) as hxed. The free 
parameters can be chosen by requiring that when the Dirac-Hartree equations are solved in 
the limit of inhnite nuclear matter, the empirical equilibrium density (p^ = 0.1484 fm“^). 
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FIG. 4. Charge density distribution for The experimental curve is from [ De87 ]. The 
Dirac-Hartree calculations for parameter set L2 yield the long-dashed curve, while those from set 
NLC yield the dot-dashed curve. 


energy/nucleon (cq = —15.75 MeV), and bulk symmetry energy (35 MeV) are reproduced.|^ 
The empirical equilibrium density is determined here from the density in the interior of 
^°®Pb and corresponds to kp = 1.30 fm“^. We also fit the empirical rms charge radius of 
^°Ca (rrms = 3.482 fm), which is determined primarily by rris- This procedure produces the 
parameters in the row labeled L2 in Table |, which are taken from ||Ho81|| . This parameter 
set yields the same values for G/ and as in Eq. (p.21|) , so that M*/M = 0.541 and 
K 545 MeV at equilibrium. 

Once the parameters have been specified, the properties of all closed-shell nuclei are 
determined in this approximation. For example. Figs. ^ through ^ show the Dirac-Hartree 
charge densities of “^^Ca, and ^osp]^ compared with the empirical distributions deter¬ 
mined from electron scattering |Pe87|| . The empirical proton charge form factor has been 
folded with the calculated “point proton” density to determine the charge density. 

In Fig. ^ the predicted energy levels in ^ospp compared with experimental values 
derived from neighboring nuclei |Bb69| , |nu79| . The relativistic calculations clearly reveal 
a shell structure; the level orderings and major shell closures of the nuclear shell model 
are correctly reproduced. This successful result arises from the spin-orbit interaction that 
occurs naturally when a Dirac particle moves in large, spatially varying classical scalar and 


®The number of significant digits in the empirical input values is not intended to indicate how 
accurately these quantities are known. We are merely reporting the precise values used in [Ho81| 
to determine the model parameters. 


18 













r (fm) 


FIG. 5. Charge density distribution for ^°Ca. The experimental curve is from |De87]. The 
Dirac-Hartree calculations for parameter set L2 yield the long-dashed curve, while those from set 
NLC yield the dot-dashed curve. 



FIG. 6. Charge density distribution for ^^®Pb. The solid curve is from | De87| . Dirac-Hartree 
results are indicated by the long-dashed curve (set L2) and the dot-dashed curve (set NLC). 
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FIG. 7. Predicted spectrum for occupied single-particle levels in ^''^Pb. Experimental values 
are taken from neighboring nuclei. 


vector fields ||Fu36|,^<7^ , |Kr74| , |Mi74|JMi75|| . Note that whereas (7s0o and g^Vo tend to cancel 


in the central potential that determines satnration, they add constructively in the spin- 
orbit potential. We emphasize that no parameters are adjusted specihcally to produce the 
spin-orbit interaction, as is usually the case in nonrelativistic calculations. Thus, with a 
minimal number of phenomenological parameters determined from bulk nuclear properties, 
one derives the level structure of the nuclear shell model. 

We emphasize that in this relativistic model of nuclear structure, the calculation of 
the ground state is self-consistent. The scalar and vector fields follow directly from the 
scalar and baryon densities, which are in turn determined by the solutions to the Dirac 
equation (|2.25 ) in the classical fields. Moreover, this relativistic shell model arises from a 
simple approximation to the underlying QHD lagrangian. Thus one has a consistent many- 
body framework to systematically investigate corrections. 


2. Deformed Nuclei 

To study the systematics of this relativistic model of nuclear structure, we extend the 
preceding equations to deal with deformed, axially symmetric nuclei. This allows us to 
calculate not only the ground states of nuclei with fully closed shells, but also those for even- 
even nuclei between closed shells. We will concentrate here on nuclei with 12 < B < 40, 
which includes the Ip and 2s-ld shells The restriction to azimuthal and reflection- 


symmetric deformations is reasonable for light, even-even nuclei. 
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These assumed symmetries of the ground state, together with the assumption of well- 
dehned parity, imply that the nonvanishing meson helds are the same ones that appear 
in spherical nuclei |[Pr87|| . Thus the Dirac-Hartree equations are essentially the same as 
those written earlier, except that all helds now depend on both a radial and an angular 
coordinate [for example, 6^)], and all the differential equations become partial differential 
equations. The source densities are still computed as in Eq. (|2.28|) , but they now depend 
on r and 6. There are several methods for solving the resulting set of coupled partial 
differential equations, and the interested reader is directed to the literature for a discussion 
|[Le8tj , PaS7| , [f 'u87|J'/h91a|| . The equilibrium deformation is obtained by choosing the occupied 
single-particle states to minimize the energy. 

The results of Furnstahl, Price, and Walker [[Fu87|| for quadrupole moments in the 2s-ld 
shell show how this observable can constrain the properties of nuclear matter near equilib¬ 
rium density. These authors use the full complement of parameters in the present model 
and constrain them to produce nuclear saturation at = 1.30 fm“^ with a binding energy 
of 15.75 MeV and a bulk symmetry energy of 35 MeV, as well as the correct rms radius 
for ^°Ca, just as for the calculations of spherical nuclei in the preceding subsection. An 
inhnite number of parameter sets will satisfy these constraints, and two examples (NLB 
and NEC) are shown in Table |. Both of these sets produce roughly equal values for M* 
at nuclear matter equilibrium {M*/M = 0.61 for NLB and M*/M = 0.63 for NLC), so 
they generate similar spin-orbit splittings and deformations. As seen in Figs. §-|^, NLC 
accurately reproduces the nuclear charge densities and splittings for spherical nuclei (results 
for NLB are similar). Moreover, the agreement with experimental quadrupole moments is 
excellent, particularly the systematic trends and the oscillation between oblate and prolate 
shapes around B = 32 (see Fig. 7 in |[Fu87|| ). Thus the successful description of spherical 
nuclei in this relativistic model can be extended to reproduce the observed systematics of 
light deformed nuclei with the same parameters. 

These authors also observe, however, that if one uses the set L2 to compute the 
quadrupole moments, the smaller value of M*/M = 0.54 at equilibrium leads to a prediction 
of spherical shapes for ^®Si, and which are inconsistent with B(E2; O’*" —>■ 2+) values 
derived from experiment ||Le75|| . Thus, the sensitive dependence of the deformation on the 
level density near the Fermi surface (which is determined essentially by the inverse of M*) 
allows one to conclude that the small adjustments in equilibrium properties afforded by the 
nonlinear scalar couplings k and A are significant. A similar analysis of the 2s-ld shell 
binding energies (see Fig. 10 in ||Fu87|]) shows that compression moduli greater than roughly 


350 MeV are too large, which is true of sets L2 {K ^ 545 MeV) and NLB {K ^ 420 MeV). 
Note that the favored parameter set NLC (with K k. 225 MeV) has a negative value for the 
quartic coupling A, which may lead to problems for large values of the scalar held, since the 
energy is no longer bounded from below. As we will see later, however, the mean-held energy 
is actually valid only as an expansion in powers of the helds, and thus global questions about 
stability are not particularly relevant. 
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C. Nucleon—Nucleus Scattering 


The scattering of medium-energy nucleons from nuclei can provide information about 
both nuclear structure and the NN interaction. Since the NN interaction has complex spin, 
isospin, momentum, and density dependence, nucleon-nucleus scattering exhibits a wide 
variety of phenomena. As a starting point for describing these phenomena, we use the Dirac- 
Hartree description of the nucleus, together with the relativistic impulse approximation 
(RIA), which assumes that the interaction between the projectile and target nucleons has 
the same form as the interaction between two nucleons in free space. This interaction is 
used to produce a nucleon-nucleus optical potential that incorporates the leading term in a 
multiple-scattering series. 

Although the simple QHD models discussed above are useful for studying the average 
properties of the nuclear interaction, they are less useful for describing the detailed quan¬ 
titative features (such as spin dependence) of the full NN scattering amplitude. These 
quantitative features are important for any reasonable description of the nucleon-nucleus 
scattering observables. The RIA allows us to combine the empirical free-space scattering 
amplitude with a relativistic calculation of the nuclear ground state. 

The RIA as originally formulated ||Mc83|JSh83|JCI83|| involves two basic procedures. First, 
the experimental NN scattering amplitude is represented by a particular set of five Lorentz- 
covariant functions |[Mc83a|] that multiply the so-called “Fermi invariant” Dirac matrices. 
The Lorentz covariant functions are then folded with the Dirac-Hartree target densities to 
produce a hrst-order optical potential for use in the Dirac equation for the projectile ||Ho91 


Here we briefly summarize the formalism and the results; a more complete discussion is 
given in [[Se86 . 

The constraints of Lorentz covariance, parity conservation, isospin invariance, and that 
free nucleons are on their mass shell imply that the invariant NN scattering operator JF can 
be written in terms of five complex functions for pp scattering and five for pn scattering. In 
the original RIA, JF was taken as 






A, 5 


(2.37) 


where the subscripts (0) and (1) refer to the incident and struck nucleons, respectively. Each 
amplitude JF^ is a complex function of the Lorentz invariants t (four-momentum transfer 
squared) and s (total four-momentum squared), or equivalently, of the momentum transfer 
q and incident energy E. It is found empirically that the amplitudes , JF^, and are 
much larger than any amplitudes obtained in a nonrelativistic decomposition, which uses 
Galilean-invariant operators. 

The RIA optical potential Uopt{(l,E) is dehned as 


U^,(q,E) = E , (2.38) 

n=l 

where JF is the scattering operator of Eq. ( |2.37|) , p is the magnitude of the projectile three- 
momentum in the imcleon-nudeus c.m. frame (where the scattering observables are calcu¬ 
lated), |T) is the A-particle nuclear ground state, and the sum runs over all nucleons in the 
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target. JF is a function of the momentum transfer q and collision energy E, which we take 
to be the proton-nucleus c.m. energy; this amounts to neglecting nuclear recoil. 

With these simplihcations, the Dirac optical potential is local, and only diagonal nuclear 
densities are needed. For a spin-zero nucleus, the only nonzero densities are the baryon and 
scalar densities of Eq. (|2.28|) , plus a tensor term computed by inserting cr°* between the 
spinors in Eq. (|2.28|) . Thus the optical potential takes the form 

Uopt = U^ + - 2ia • fU^ , (2.39) 

where = U^{r]E) for each component. The tensor contribution U'^ is small and is 
neglected in what follows. The RIA optical potential then has only scalar and vector contri¬ 
butions, and the Dirac equation for the projectile has precisely the same form as in Eq. ( p.25D , 
with replacing g^Vo and replacing {—gs(l>o)- 

hWo(x) = {-ia-V+ U^{r;E) +f3[M+ U^{r;E)]]Uo{^) = EUo{^) . (2.40) 


In practice, one includes in the Coulomb potential computed from the empirical nu¬ 
clear charge density; other electromagnetic contributions arising from the proton anomalous 
magnetic moment are of similar size to the tensor term U'^. 

Since representative RIA results have been shown many times in the literature (see, for 
example. Figs. 25-28 in ||Se86||), we will not reproduce the hgures here.[] The target densities 


in these calculations are taken from the results for spherical nuclei discussed above, with no 
further adjustment of parameters, and the RIA calculations agree remarkably well with the 
data. Moreover, when compared with nonrelativistic impulse-approximation calculations, 
the relativistic results are superior, particularly for the spin observables. Although the 
nonrelativistic results improve when higher-order corrections are included ||Co90| , pu90|] , one 
concludes that the important spin- and density-dependent effects are already contained in 
the relativistic impulse-approximation framework ||Pi84| , [EI}^ , |Lu87|| . 

The spin dynamics is inherent in the relativistic formalism and arises naturally from the 
large Lorentz scalar and vector potentials in the Dirac equation ( |2.40|) . This is precisely the 
same spin dynamics that produces the observed spin-orbit splittings in the bound single¬ 
particle levels. Thus the relativistic Hartree calculations provide a minimal unifying theoreti¬ 
cal basis for both the nuclear shell model and medium-energy proton-nucleus scattering—two 
essential aspects of nuclear physics. 


D. Nuclear Excited States 

We now turn to the calculation of nuclear response functions and the properties of nuclear 
excited states, as described in the random-phase approximation (RPA) built on the MFT 
ground state. These excitations arise from the consistent linear response of the ground state, 
in which the nucleons move coherently in varying classical meson helds that are in turn 


^Computer codes for performing RIA calculations are also readily available [Po91|. 
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determined by oscillatory nuclear sources. Many studies of the relativistic nuclear response 
have been carried out for both inhnite matter and hnite nuclei, beginning with the pioneering 
work of Chin ||Ch77|| , and here we will focus on two basic issues. The hrst is the importance 


of consistency, which implies that the particle-hole interaction in the excited states must 
be the same (and use the same parameter values) as the interaction in the ground state. 
This ensures that the occupied and unoccupied single-particle states are indeed orthogonal. 
Second, the relativistic response involves not only the familiar positive-energy particle-hole 
conhgurations, but also conhgurations that mix positive- and negative-energy states. These 
new conhgurations are crucial for the conservation of the electromagnetic current and the 
separation of the “spurious” = 1“ state. This emphasizes that the Dirac single-particle 
basis is complete only when both positive- and negative-energy states are included. 

The calculation of the linear response is basically the same as in nonrelativistic many- 
body theory |[Fe71j| . The principal idea is to compute the particle-hole (polarization) prop¬ 
agator and to extract the collective excitation energies and transition amplitudes from the 
poles and residues of this propagator. Several methods have been developed and applied to 
hnite nuclei |[hu85|JFu8^JJNi86|JWe87| , [We88| , |hu88a| , |B188|JSh89|JDa90|JHo90| , [Pi90| , [Pr92a|| . Initial 
calculations were restricted primarily to isoscalar excitations, because htting bulk nuclear 
properties constrains only the isoscalar particle-hole interaction signihcantly, and because 
the isovector response depends critically on pion dynamics, with the associated complexity 
that we discuss later. 

An excellent discussion of the role of consistency is contained in the work of Dawson and 
Furnstahl |Pa90|| , where results for the low-lying, negative-parity, isoscalar states in 

and ^‘^Ca are compared to several empirical levels that might be reasonably described 
as particle-hole excitations. The full RPA eigenvalues agree favorably with the empirical 
values, which is a nontrivial result because of the large cancellations between scalar and 
vector contributions. Moreover, it is found that the negative-energy states play an important 
role in determining the RPA spectrum, particularly for the spurious 1“ state. Lorentz 
covariance implies that in a consistent RPA calculation, this state should appear at zero 
excitation energy, which occurs only when the full Dirac basis is maintained. 

Moreover, RPA calculations by Furnstahl ||Fu89b|| of the ratio of the transition charge 
density to the longitudinal current for a particular (3“, 0) excitation in show that if only 
the positive-energy (particle-hole) conhgurations are retained, the electromagnetic current 
is not conserved. In contrast, with the full RPA calculation, including the contributions 
from negative-energy states, current conservation is restored. The conclusion is that it is 
essential to include all states in the Dirac basis to maintain the conservation of the current. 


E. Nuclear Matter at Finite Temperature 

The preceding discussion of the nuclear matter equation of state was restricted to zero 
temperature. The extension to hnite temperature is straightforward in the MFT, since the 
hamiltonian is diagonal and the mean-held thermodynamic potential D can be calculated 
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exactly.0 The results for the scalar density, baryon density, energy density, and pressure are 


given by (here k = A = (?„ = 0) |Se 86 | 


'J f Af* 

(2.41) 

PB = Uk) , 

(2.42) 


(2.43) 


(2.44) 

where the baryon and antibaryon distribution functions are 


\nk[T,v)J 1 J 

(2.46) 

and the reduced chemical potential is z/ = /r — (We set Boltzmann’s constant /cb = !•) 

The appropriate value of M* is determined by minimizing the thermodynamic potential 
with respect to that parameter: 


' dn 
dM* 


= 0 . 


(2.46) 




The nuclear matter equation of state at all densities and temperatures for this hadronic 
MFT model (QHD-I) is shown in ||!Se92|,|WaU5|] . The model has also been combined with a 


simple description of quark-gluon matter to describe the hadron/quark phase transition as a 


function of temperature and density |^e 86 || . A more extensive examination of the equation of 


state, as well as a discussion of the covariance of the hnite-temperature results, is contained 


in Fu90 . 


3. PION DYNAMICS AND CHIRAL SYMMETRY 

The relativistic neutral scalar and vector helds are the most important for determining 
the bulk properties of nuclear systems. Nevertheless, the lightest and most accessible meson 
is the pion, whose interactions with nucleons and nuclei have been extensively studied at the 
meson factories. It is therefore impossible to formulate a complete and quantitative hadronic 
theory without including pionic degrees of freedom. Here we briefly review some important 
aspects of pion dynamics in QHD. More complete discussions, together with references to 
the original literature, can be found in ||be86|,PeU^jWa95| . 


®We neglect the zero-point corrections from the Dirac sea in this section (see 
as well as thermal contributions from the massive isoscalar mesons. 


|Fr7S| and [Fu91|), 
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In the limit of massless u and d quarks, QCD possesses a global, chiral SU(2)i x SU(2)i^ 
symmetry. Chiral transformations may be written in terms of a set of vector {V) and axial- 
vector (A) generators, which produce corresponding isospin rotations. The symmetry is 
spontaneously broken, leading to the existence of pseudoscalar Goldstone bosons (pions). 
The vector (isospin) symmetry, which forms an SU(2)y subgroup of the original chiral group, 
remains unbroken. Our effective hadronic theories should respect these underlying symme¬ 
tries, which have important consequences for the way mesons interact with themselves and 
with each other; a thorough discussion is contained in |[Wa95|| . In nature, electromagnetic 


interactions and hnite quark masses imply that these symmetries are only approximate, but 
the symmetry-violating terms can be added as small perturbations. 


A. The Linear Sigma Model 


The simplest model illustrating these ideas is the linear sigma model |^c57| , |Ge6^ , [Le72|| , 
which contains a pseudoscalar (75) coupling between pions and nucleons, plus an auxiliary 
scalar field (denoted here by s) to implement the symmetry. A small symmetry-violating 
(SV) term is included to generate a hnite pion mass. We will also include a massive, neutral, 
isoscalar vector held to supply a repulsive interaction, as in QHD-I. 

By demanding that the theory be Lorentz covariant, parity invariant, isospin and chiral 
invariant, and renormalizable, one is led to the form 

= -^chiral + 'CsV 5 (3.1) 

>^^chirai = - 9vV^) - gn{s + • 7r)]i/> + ^{d^sd^s + • d^ir) 

-iA(s2 + 77^ - vY - + 5C , (3.2) 

= es • (3.3) 

Here ip, tt, and s are the nucleon, isovector pion, and neutral scalar meson helds, respectively, 
and is the pion-nucleon coupling constant. The parameters A and v describe the strength 
of the meson self-interactions, and e is a small chiral-symmetry-violating parameter related 
to the pion mass; the exact chiral limit is obtained by setting e = 0. The form of the meson 
self-interactions allow for spontaneous symmetry breaking, which is used to give the nucleon 
a hnite mass. 

The lagrangian (p.2|) is invariant under global vector and axial-vector isospin transfor¬ 
mations, which imply (by Noether’s theorem) the conserved isovector currents^ 


H 


77 X (9^77 , 


= |'0757^t-V^ -I- — sd^Tz 

in the chiral limit e = 0. When e 7^ 0, we obtain instead the PCAC relation 

df,A^ = 677 , 


(3.4) 

(3.5) 


(3.6) 


®The (isoscalar) baryon current is also conserved. 
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which follows from the held equations. Note that the chiral symmetry is realized linearly, 
which means that under a general chiral transformation, neutrons mix with protons, and 
the scalar mixes with the pions. (The vector held is a chiral scalar.) Moreover, the linear 
symmetry requires that the scalar and pion couple to the nucleon with equal strength ^ 
13.4. In applications, one often relaxes this condition (with no justihcation) and sets = 
g-n/gAi where the axial coupling gA ~ 1.26; this ad hoc procedure allows the Goldberger- 
Treiman relation to be satished at the tree level ||We67| , [Ce72|| . 

The baryon mass is generated by spontaneous symmetry breaking |[Le72| , [5e86| , which 
implies that the scalar held s has a nonzero vacuum expectation value Sq, and that the pion 
is a massless Goldstone boson (in the limit e = 0). In terms of the shifted scalar held 


a = So — s 


and the physical masses dehned by 


M = g^so , 

the lagrangian reads 


M 


e = 




A = 


m: 


mz 


2 M 2 


(3.7) 


(3,8) 


Cau = - 9.7(.W - (M - t/,tr) - it/,75T • Tr]i> + ^(di^ad^a - 

- + 7(9,.TT • d^TT - mln^) 


+ 971 


mz 


mz 


2M 


a(a^ + TT^) 


9 


2 ml-ml 


8 M 2 


'2 I ,^272 

a +7r 


(3.9) 


where the counterterms 6C will here and henceforth be suppressed. Note that the explicit 
chiral symmetry violation is now contained entirely in the parameter ml, the square of the 
pion mass. Moreover, if we consider the (renormalized) parameters M, m.^, m^, and g.,^ as 
experimentally known, then apart from the vector meson coupling g„, there is only one free 
parameter in this model, namely, the chiral scalar mass m^j. Both the signs and magnitudes 
of the nonlinear meson interactions proportional to and are determined by the chiral 
symmetry and its spontaneous breaking in the vacuum. 

A vector-isovector p meson can be added in a renormalizable manner to produce QHD- 
II, as discussed in [ Se86 ].H At the mean-held level of interest here, the rho meson enters 
as in Eq. ( [^.36|) . Given the similarity between the lagrangians in Eqs. ( |3.9|) and ( [^.ID , it is 
only natural (but incorrect!) to identify the a held with the scalar held 0 studied earlier. 
Since the pion held vanishes at the mean-held level (we assume our systems possess the same 
symmetries as those in Section |^), the dynamical equations are exactly the same as before, 
except that the nonlinear scalar interactions are now determined by the symmetry: 


Kn 


M 


^977 

'JP 


{ml - ml) , 


\ ‘^977 ( 2 2 ' 

A... = 


(3.10) 


^^Including the isovector p meson destroys the chiral symmetry. A fully chiral, renormaliz¬ 
able model (QHDMII) containing both the p and ai mesons is derived in [pe92b |. See also 
| Ko93a , Fo95a|| . 
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Unfortunately, the identification of a with 0 leads to serious difficulties. First, even 
though there are still two free parameters in symmetric nuclear matter, it is impossible to 
reproduce the empirical equilibrium point, as documented long ago by Kerman and Miller 


(Their result is reproduced in Fig. 5 of ||Fu96|| . 


The basic problem is that the 
nonlinear scalar interactions are large and have the wrong signs. [Compare naw < 0 and 
> 0 in Eq. ( p.lO|) with the favored NEC values in Table |.] This problem was solved 
in ||Bo83|J5aS5| , where interactions between the scalar and vector helds were introduced by 


making the replacement 






(3.11) 


in Eq. ( p.2|) . Now the vector held also acquires its mass through spontaneous symmetry 
breaking, when the scalar held is shifted as in Eq. (|3.7|) . This results in the replacement 








M 


_ glrnl 

2 M2 




(3.12) 


in Eq. ( p.9|) . The signs of the scalar-vector cubic and quartic interactions are opposite to 
those in Eq. (|3.9|) , which now allows the nuclear matter equilibrium point to be reproduced.]^ 

Nevertheless, an extensive mean-held analysis shows that it is impossible to generate 
realistic results for nuclear densities, single-particle spectra, and total binding energies within 
this framework |[Fu93|JFu93a|jlu9^ . The primary problem is still the existence of scalar 
nonlinearities with incorrect systematics; in addition, the large value of the (known) scalar 
coupling leads to a scalar mass that is much larger than that required by the phenomenology 
{nifj ^ 800 MeV compared to a desired value of rria ~ 500 MeV). 

We conclude that the standard form of spontaneous chiral symmetry breaking, imple¬ 
mented in a model with a linear realization of the symmetry, cannot produce successful 
nuclear phenomenology at the mean-held level when we identify the chiral scalar a with 
the scalar in QHD- IQ Note that the interaction potential (the term proportional to A) 
in Eq. (|3.2|) is the only form consistent with both chiral symmetry and renormalizability. 
This failure of the sigma model is therefore evidence, even at the mean-held level, that the 
simultaneous constraints of linear chiral symmetry and renormalizability are too restrictive. 

A diherent approach to the chiral symmetry is more compatible with observed nuclear 
properties and with successful relativistic mean-held models. Here one takes the chiral scalar 
mass rOfy to be large to eliminate the unphysical scalar nonlinearities and then generates 
the mid-range attractive force between nucleons dynamically through correlated two-pion 


^^The new scalar-vector interactions destroy the renormalizability of the model, but we will 
overlook that for the moment. 


^2lt has recently been shown |He94a JCa96[ that realistic nuclear systematics can be achieved 
in a linear realization by using a logarithmic, nonrenormalizable scalar potential to generate the 
spontaneous symmetry breaking. 
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exchange.0 Numerous calculations [|Ja75|jDu77|Pu8d|JLi89|,pi^ , [Ki94|| have shown that the 
exchange of two interacting pions in the scalar-isoscalar channel produces a strong attractive 
force with a range comparable to that of a scalar meson with a mass of 500 to 600 MeV. In 
principle, this approach can be realized within the linear models discussed above, by reject¬ 
ing the mean-held approximation (which is inadequate) and by retaining contributions from 
correlated two-pion exchange between nucleons. This clearly introduces much greater com¬ 
plexity, since one must hrst construct a boson-exchange kernel containing correlated two-pion 
exchange and then allow this kernel to act to all orders (for example, in a ladder approx¬ 
imation) to determine the NN interaction and the resulting nuclear matter energy density 
|[Ma89|| ■ This calculation is further complicated by the need to maintain chiral symmetry 
at hnite density, which is difficult to do when one uses the non-derivative (“pseudoscalar” 
or PS) ttN coupling implied by the linear realization of the symmetry 


These complexities motivate us to look for a more efficient way to implement the chiral 
symmetry. 


B. A Nonlinear Realization of Chiral Symmetry 


The chiral dynamics of pions and the successful mean-held picture described in Section 
can be combined if the chiral symmetry is realized in a nonlinear fashion. In a nonlinear 
realization, chiral transformations mix the nucleon and pion helds, and a scalar-isoscalar 
held is unchanged. There are many ways to implement a nonlinear realization; here we 
follow the approach of Weinberg |[We67|| and discuss a diherent method in Section 
First notice that a chiral rotation matrix can be written as 


exp(|a; • T 75 ) = 1 cos(a;/2) -f ifi • rys sin(ci;/2) , (3.13) 

where 1 denotes the unit matrix in spin-isospin space and uj = nuj is real. The key obser¬ 
vation is that the linear terms in a and tt in Eq. ( B.9D can be written as a chiral rotation 
matrix: 


M - g^a + ig^T • 7775 = [(M - g^af exp(i0 • t'^^) 

= [{M — g-ircry + gln^]^^^[l cos 6 + in • T75sin6'] , 

where tt = nn and 0 = nO. It now follows immediately that 

M - gnO- . „ g^Ti 


cos 9 = 


HM - ’ 


sin 9 = 




(3.14) 


(3.15) 


One can then eliminate the linear pion-nucleon interaction in the lagrangian C^uj by dehning 
a new baryon held iV as a unitary transformation of the old baryon held 


^^One can see that the scalar self-interactions become small by writing the lagrangian ( |3.9| ) in 
terms of the re-scaled field y = mo-cr. Note, however, that the cttt^ coupling becomes large when 
nia becomes large. 
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N = exp(|0-r75)V’ ^ ^J^{1 + cosO) (1 + , 


(3.16) 


where 


$, = n-. 


sin 6 


' 1 + cos 9 

and by defining new meson fields tt' and a' as 

^ _ QtT / _ / / ^ 

= M - g^a = -— 

2M cos 9 


(3.17) 


(3.18) 


Although the algebra is tedious, it is straightforward ||Wa95|| to rewrite C^ui in terms of 
the new fields as 


C.. = -{M- g.<y') 

+ 1 I /// \2 12 [(/M7r)7^75T • d^ir' - • tt' X a^7r']}iV 

+ ^(d^a'd^a' - 

+ {ml - ml)[{f /m^)a'^ - |(//m^)V'^] - + \mlV^V^ , (3.19) 


where we have defined the ratio 


R = R{a\ Tz') 


1 - 2{f /m^)a' 

1 + (//m„)27r'2 ■ 


(3.20) 


This somewhat imposing lagrangian has several important advantages. First, the linear PS 
coupling between the nucleon and pion has been replaced by a derivative or pseudovector 
(PV) coupling and a so-called “sea gull” term, where the nucleon couples to two pions 
simultaneously. In fact, there are an infinite number of such couplings, when one expands 
the prefactor [1 -|- (//m7r)^7r'^]“^ as a power series in tt'^. Nevertheless, all of these couplings 
involve at least one derivative acting on the pion field; thus, in the limit of vanishing pion 
momenta, the pions and nucleons decouple, as required by the chiral soft-pion theorems that 
follow from the SU(2)i x SU(2)/^ current algebra and PCAC |[Sa69|,[no^ . 

Second, the new pseudovector coupling constant is 




(3.21) 


which is much smaller than the pseudoscalar coupling constant gl/Fi^ 14.4. Moreover, the 
explicit derivative couplings eliminate the sensitive cancellations between Feynman diagrams 
that are necessary in the linear realization. The large coupling and sensitive cancellations 
show that the linear realization is an inefficient way to implement the symmetry. 

As noted earlier, m„ should be taken to be large, to avoid the unwanted nonlinear 
interactions. If desired, can be kept finite, so that it plays the role of a regulator that 

[|Ma82|| 


maintains the renormalizability of the model 


or m^j can be taken to infinity, so 


that the chiral scalar field decouples, resulting in the effective nonlinear model of Weinberg 
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||We67|JWa95|| .n If the process tt + tt—>7 r + 7 ris investigated in this model with a heavy 
a, and the chiral-invariant Born amplitude is unitarized, one observes a broad, low-mass, 
near-resonant amplitude in the (O’*", 0) channel, even though the chiral a has a large mass. 
When this model tttt scattering amplitude is included in the two-pion-exchange part of the 
NN interaction, the result is a dynamically generated, broad, low-mass (~ 600 MeV) peak 
that resembles the exchange of a light scalar meson ||Li89|,^n^ . 


This strong scalar-isoscalar two-pion exchange can be simulated by introducing a low- 
mass, effective scalar held (j) coupled directly to the nucleon, and scalar self-interactions can 
be added to include a density dependence in the mid-range NN attractive force. We empha¬ 
size that the purpose of the effective held is to parametrize the NN attraction, so that we 
can avoid the complicated calculation of scalar-isoscalar pion loops. All scalar propagation 
will be restricted to spacelike momenta, and thus scalar particles are always virtual.|^ At 
the mean-held level, the classical scalar held is an efficient way to incorporate the ehects of 
pion exchange that are the most important for describing bulk nuclear properties. 

If we now take rriu —>■ cx), the resulting nonlinear, chiral, ehective lagrangian is (primes 
on the helds are omitted) 






- (M - g^cf) 

1 


)27r2 


[{f • 9^77 - • TT X 9/,7r]}A^ 


1 -h {f/m^ 

+ \{d,cfd^ct> - ml4>^) - + \mlV,V^ 

+ 2(1 + [ l + (/M^ • a"- - 


‘77^ 


(3.22) 


The explicit symmetry-violating term proportional to has been included, and since 
f /m.,^ = gT^/2M from eq. (^.211) , the exact chiral limit (m^r —> 0) is sensible. 

We retain the notation N for the baryon held to remind us that it transforms nonlinearly, 
in contrast to the held used earlier. (Linear transformations are independent of the pion 
held, while nonlinear transformations depend on the pion held.) Although the nonlinear 
transformation law implied by Eq. (|3.16 ) is complicated, it has been shown [ We67 , Co69 
that under a general chiral transformation. 


N —y N' = h{x)N , 


(3.23) 


contrast to the lagrangian (^), if one rewrites the nonlinear lagrangian ( 3.19| ) in terms of 
y; = macr', it is easy to see that all interactions involving the scalar field vanish in the nia oo 
limit. 


^^The effective scalar-isoscalar field is introduced to simplify the description of the attractive NN 
interaction. Thus it plays a role in problems with B > 1 and cannot be directly interpreted as an on- 
shell particle. The idea of introducing virtual degrees of freedom to describe fermionic interactions 
is not revolutionary; phonon fields in metals and the Ginzburg-Landau field in quantum liquids 
have played such a role for many years with considerable success [Fe71]. 
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where h{x) G SU(2)y is an isospin rotation matrix that depends locally on the pion held. 
Thns the inclnsion of a Ynkawa conpling 

Cy = g^NNcj) (3.24) 


(as well as scalar self-conplings) in £eff leaves the chiral symmetry intact, becanse the light 
scalar is an isoscalar, and the nncleon held transforms as in Eq. (|3.23|) | Wa95|| .[^ Moreover, 


since the pion mean held vanishes, the chiral mean-held theory obtained from C^s prodnces 
precisely the same held eqnations and energy density as in Section ^ Finally, althongh the 
conpling strength Qs of the light scalar is comparable to (as verihed by explicit calcnlation 
of the correlated two-pion exchange [|Li89| , tLi90|] ), we no longer reqnire gs = gTT and are free to 
adjnst gs within a reasonable range. The nonlinear scalar interactions, which parametrize the 
density dependence of the correlated two-pion exchange, also contain adjnstable parameters. 

Thns, throngh these somewhat lengthy argnments, we conclnde that the mean-field 
QHD-I model studied earlier is consistent with chiral symmetry, provided we think in terms 
of a nonlinear realization of the symmetry. The importance of the resulting scalar-isoscalar 
mean held and optical potential in producing a successful nuclear phenomenology was illus¬ 
trated in Section ^ Although the one-pion-exchange interaction (and its iteration) produces 
a relatively small contribution in nuclear matter due to its spin dependence, the correlated 
two-pion-exchange contribution has a spin-independent, isoscalar part that is large and that 
must be included from the outset.]^ The discussion in Section ^ shows that the scalar- 
isoscalar held in QHD-I is an efficient and successful way to incorporate these important 
pionic ehects. A profound change has occurred, however, because in contrast to the original 
proposal of QHD-I as a renormalizable field theory, we are now forced to consider the scalar 
field as an effective degree of freedom and the lagrangian in Eg. (1^) as a nonrenormalizable 
effective lagrangian. We will consider this more general strategy for QHD in Section |^. 


C. The A Isobar 

The dominant phenomenological features of the low-energy ttN interaction are that low- 
momentum pious interact weakly with nucleons (they decouple as 0), and that the 

interaction is dominated by the hrst pion-nucleon resonance, the A(1232). This resonance 
represents the hrst excited state of the baryon, with (J’^,T) = (| , |). It is essential to 


^®Note that one cannot add a Yukawa interaction fnfcl) with an effective scalar field to the linear 


sigma model [Eq. (3.£)], as this destroys the chiral symmetry. 


^^In the few-nucleon problem, the situation is reversed. Because of the more complicated spin 
dependence of the wave function, matrix elements for one-pion exchange no longer average to zero. 
Moreover, the existence of an almost-bound state at essentially zero energy in the NN channel 
[|Fe71 1 implies that there is a nearly exact cancellation between the attractive (scalar) and repulsive 
(vector) parts of the NN force in few-body systems. Thus, the pion-exchange “tensor force” plays 
a more prominent role. 
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include this degree of freedom in the theory to generate realistic results for pion-nucleus 
interactions and for few-nucleon systems. It is impossible to put a held with these quantum 
numbers into a renormalizable lagrangian, but it has been shown that the A degree of free¬ 
dom can be produced dynamically within the model. The most efficient way to incorporate 
the A is as an effective degree of freedom. 

Here we simply summarize the results, as the arguments and references to the original 
literature are detailed in ||Se86|,|Se9^,|Wa95[1 . The sum of ttN ladder diagrams with nucleon 


exchange can be investigated within the framework of the chiral vrN theory discussed above. 
Partial-wave dispersion relations can be used, with the one-baryon-exchange mechanism 
considered as the driving term, and the resulting integral equations solved with the N/D 
method [[Li91|| . This is a relativistic extension of nonrelativistic Chew-Low theory. As in 
the Chew-Low theory, a resonance is found in the (| , |) channel. Evidently, the box 
diagram in the ladder sum involves a loop integral, which obtains signihcant contributions 
from large loop momenta or short distances; thus, the predicted position of the resonance is 
sensitive to the approximations made. In contrast, the predicted resonance width is much 
less sensitive. As a result of the work summarized in | Se92a , Wa95| , it is clear that the hrst 
excited state of the nucleon, the A(1232) with (I"*",!), which is the dominant feature of 
low-energy pion-nucleus interactions, is generated dynamically in QHD. Thus it can also 
be included as an effective degree of freedom in a nonrenormalizable hadronic lagrangian 
Pl8^Pe92a| , [We93i[Hi96| , nM| . 


4. FEW-NUCLEON SYSTEMS 


The purpose of this section is to present some additional evidence of the successes of 
hadronic held theory by considering the few-nucleon problem. As before, we want to focus 
on basic ideas and phenomenology that will guide us in the construction of an effective held 
theory for the many-body problem. We begin by summarizing modern models of the NN in¬ 
teraction and show that when applied to two-nucleon systems, accurate results are obtained 
both for electromagnetic observables and for threshold pion production. We also consider 
the extension of the two-nucleon problem into nuclear matter, by discussing the relativis¬ 
tic generalization of the independent-pair approximation ||Fe71|JWa95|| , which is often called 


Dirac-Brueckner-Hartree-Fock (DBHF) theory. This will allow us to compare nuclear mat¬ 
ter calculations involving two-nucleon correlations to the mean-held results discussed earlier 
and also to understand the relationship between the relativistic two-nucleon interaction and 
nonrelativistic many-body forces. 


A. The Nucleon—Nucleon Interaction 


Models of the NN interaction have been studied seriously for nearly 50 years. An excellent 
historical review is contained in ||Ma89|] . Modern meson-exchange models provide excellent 


hts to experimental observables up to and somewhat beyond pion-production threshold, 
which is at a laboratory kinetic energy of roughly 300 MeV. Once the models are calibrated 
to the NN data, they can be used to study both the interaction of few-nucleon systems 
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with experimental probes and the nuclear many-body problem. As mentioned in the In¬ 
troduction, to achieve accurate reproduction of the data, several mesons are needed [the 
most important are the 7r(0“, 1), (t( 0+,0), a;(l“,0), and p(l“, 1), although the 77 ( 0 “, 0) and 
5(0^, 1) are often included], and many models contain both N and A degrees of freedom. 
There are models based entirely on boson exchange ||Zu81| , |CTa87| , [fe87| , [IVla89|JGr90| , |P194|| , as 


well as models that are supplemented by dispersion relations ||Br76|,|Pa79|] or by Regge theory 
Ja73|,I8t94||. (A comparison of the different models is contained in |[lVfa94e|| .) These models 


are not renormalizable and contain form factors at the meson-nucleon vertices that are not 
usually expanded in powers of momentum, as one would do in a strict implementation of 
effective held theory ||We91| , |We92a| , |Qr9^JVa93| , |Or96|| . (In other words, the interactions are 
nonlocal.) The range parameters of the vertex functions are treated as adjustable, and they 
are important for producing realistic results for observables. 

The NN system is usually studied with a quasipotential approach, in which one starts 
with the two-nucleon Bethe-Salpeter equation and then reduces it to a three-dimensional 
integral equation by dehning an approximate two-nucleon propagator. In free space, the 
Bethe-Salpeter equation can be written schematically as ||Bj65| 


T = K 


KGGT 


(4.1) 


where K is the full two-body scattering kernel and G is the exact baryon propagator. 
(For illustration, we consider only nucleons here.) The four-dimensional integral implied 


in Eq. ([4.1|) can be reduced to a three-dimensional integral by dehning a unitarized, two- 
particle propagator go and a quasipotential U: 


T = U 


UgoT 


U = K + j K(iGG - 3o)U . 


(4.2) 

(4.3) 


These two equations are equivalent to (IT), and if one solves both of them, the results 
should agree with the Bethe-Salpeter results for any choice of go. In practice, however, 
only Eq. (|4.2|) is retained, under the assumption that the corrections from ( [f.3|) are small; 
moreover, the quasipotential U is approximated by keeping just the one-boson-exchange 
“ladder” kernel (which we call V), or the ladder terms plus some “box” and “crossed-box” 
amplitudes |[lVfa89| . 

To study the electromagnetic response of the two-nucleon system, the quasipotential 
is htted to reproduce the NN phase shifts and the properties of the deuteron, and then 
the coupling of the electromagnetic current is introduced. Since the quasipotential can be 
dehned in an inhnite number ways, it is important to test the sensitivity by utilizing dif¬ 
ferent approaches; currently, calculations using three different quasipotentials are available 
[Hu9CI|JDe93| , |Va95|| . Calculations of the deuteron charge and magnetic form factors are in 
good agreement with experiment out to three-momentum transfers of roughly q'^ 80fm“^, 

provided that meson-exchange currents are included. Moreover, calculated results for the 
tensor polarization (T 20 ) of the deuteron agree with experiment out to ~ 20fm“^. Signif¬ 
icant effort has been spent in deriving meson-exchange currents that guarantee electromag¬ 
netic gauge invariance, even when phenomenological form factors are used at the vertices 


34 

































Gr87|,Po94a||. Nevertheless, there are still some open questions regarding the consistent in¬ 


clusion of realistic off-shell nucleon form factors and the model dependence of the transverse 
parts of the meson-exchange currents. 

Meson-exchange currents are also important for describing pion production near thresh¬ 
old in the reaction pp —>■ A precise measurement of the total cross section for this 

process has recently been performed [|Me9CI|jMe92|| , but calculations that include only single¬ 
nucleon mechanisms severely underestimate the measured values ||Ko66|J5cB9|. In contrast, 


the inclusion of a scalar-meson exchange current (and smaller contributions from other 
mesons) increases the cross section by about a factor of hve and leads to excellent agreement 
with the data ||Le92|,[Me92|,|Ho94B| . The results are insensitive to changes in the potential that 


generates the two-nucleon wave function and to different choices for the phenomenological 
meson-nucleon form factors. Although some questions have recently been raised about pos¬ 
sible competing mechanisms, alternative calculations cannot reproduce the measured cross 
section ||(Jo9ti| . 


B. Two-Nucleon Correlations in Nuclear Matter 


The extension of the quasipotential approach to nuclear matter is known as Dirac- 
Brueckner-Hartree-Fock theory ||An83|,po84a|JBr8^JHo87| , P?e87| , |Ani92|| . The scattering ma¬ 


trix T in Eq. (|4.2|) is replaced by a reaction matrix T, which is determined by the quasipo¬ 
tential equation 


T = V+ VgT , 


(4.4) 


where V is the (ladder-approximated) quasipotential and 5 ^ is a unitarized, two-nucleon 


propagator that includes interactions with the surrounding medium |Pdo84a|JSe86|,|Ho87[| . In 


practice, one solves for matrix elements of T, so that the driving term on the right-hand side 
involves matrix elements of V] the new ingredient is that the Dirac wave functions must 
be determined self-consistently, since the nucleons have large scalar and vector self-energies, 
analogous to the mean helds studied earlier. 

If one writes the nucleon self-energy S(/c) in (the rest frame of) nuclear matter in terms 
of scalar (E'^), timelike vector (E°), and three-vector (S'") parts ||Se 86 || : 


E(fc) = T.\k) - + 7 • k S''(fc) 


(4.6) 


then the self-energy is determined by summing effective direct and exchange interactions 
between nucleon pairs, which may be written schematically as 

nk) = E l(kp|r|kp) - (kpiripk)]. (4.6) 

This is indeed a self-consistency condition, since the self-energy determines the Dirac wave- 
functions and the reaction-matrix elements that appear on the right-hand side, and it is in 
turn determined by these quantities. The S® and components are analogous to the mean 
helds {—gs(f>o) and {—g^Vo) studied earlier; S'" is a new contribution that arises because 
exchange diagrams are included. 
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The numerous approximations and procedures that go into solving the DBHF equations 

iHo84a|,P?84i B^ JHo87|,|Te8^ , |Am92|| 


are discussed at length in 


The most important ap¬ 


proximation is the use of a quasipotential, which implies that the four-dimensional integral 
in the box diagram that forms the fundamental unit of the ladder sum is reduced to a 
three-dimensional integral. There are many possible reductions that can be identihed as the 
relativistic extension of the nonrelativistic Bethe-Goldstone equation |[Wa95|| ; unfortunately, 
numerical results for the nuclear matter binding energy are sensitive to the reduction used. 
Even after this reduction, calculated results are sensitive to the high-momentum part of 
the loop integrals, implying important contributions from baryon transitions to states lying 
well above the Fermi surface. Phenomenological form factors inserted at the vertices signih- 
cantly reduce this sensitivity, but questions about the off-mass-shell and density dependence 
of these form factors remain to be studied. Moreover, since the self-energy E enters explic¬ 
itly in the self-consistent baryon spinors, the binding energy of nuclear matter is sensitive 


to the self-consistency condition | Le89a ]. At present, the construction of a self-consistency 
condition that leads to a conserving approximation when relativistic ladder diagrams are 
summed is still a controversial topic ||Po88| , [DeM| , |Hu93|| . 

It must be emphasized that many of the conventional approximations used in DBHF 
theory have no counterparts in nonrelativistic Brueckner-Goldstone theory, and thus the 
sensitivity to these approximations has never been systematically tested. Nevertheless, it is 
still possible to make three important qualitative statements about the effects of relativistic 
two-nucleon correlations in nuclear matter: 


1. Although the correlation corrections produce changes in the RMFT nuclear matter 
binding energy that are of the same order as the binding energy itself, the corrections 

Thus the DBHF 


to the large RMFT scalar and vector self-energies are small [Ho87 


self-energies are essentially the same size as the scalar and vector mean fields studied 
earlier. (Full DBHF calculations yield typical results of 0.55 < M*/M < 0.65 at 
equilibrium density ||iVfa89| , Pr90a|jDe9I| , |Am92|| .) Moreover, the momentum dependence 
(or state dependence) of E® and E° is small; these self-energies are essentially constant 
for occupied states in the Fermi sea. Finally, the new term E'^ is small. The conclusion 
is that the successful RMFT picture presented above involving large, constant, Lorentz 
scalar and vector helds persists when two-nucleon correlations are included. □ 

2. The depletion of the Fermi sea due to correlation effects is considerably smaller in the 
relativistic framework than in the nonrelativistic framework pa^] . 


3. The self-consistency condition ( |4.6| ) introduces a density dependence into the effective 
interaction F that goes beyond what is included in nonrelativistic Brueckner-Goldstone 
theory. Because of this extra density dependence, it is possible to simultaneously £t 
both the NN phase shifts and the nuclear matter equilibrium point at the two-hole- 
line level ||Ma89|,pr90a| , Pe9l[| . Although an accurate calculation of the nuclear matter 


^®The large scalar and vector self-energies are also consistent with recent analyses based on QCD 
rules I Co91a , Co92a , Fu92a , Ji93a , |Ji94a |. 


sum 
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equilibrium point is difficult, due to the sensitive cancellations, and there are still 
open questions about three-hole-line corrections, true many-body forces, and the role 
of the quantum vacuum, this successful result cannot he obtained in a nonrelativistic 


framework with two-body potentials |Pa83| , pa85 
to make this result more conclusive. 


Further investigations are needed 


C. Relation to Nonrelativistic Calculations 


To close this section, we briefly examine why the DBHF theory produces more favorable 
results for nuclear matter saturation at the two-hole-line level than corresponding nonrela¬ 
tivistic Brueckner-Goldstone calculations. It is well known that nonrelativistic calculations 
based on realistic NN potentials predict equilibrium points that have either the correct den¬ 
sity but too little binding energy, or the correct binding energy at too high a density; this 
produces the familiar “Coester line”. What is needed is an additional, density-dependent 
repulsive mechanism, and the new saturation mechanism discussed in Section |] is precisely 
of this type. Because of Lorentz covariance and self-consistency, as the nuclear density in¬ 
creases, the nucleon effective mass M* decreases, and the nucleon velocities increase. This 
weakens the attractive force, and the net result is an increased repulsion. 

Although this additional repulsion is often described as a purely “relativistic” effect, it 
is easy to see that in a nonrelativistic framework, some of it can be attributed to three- and 
many-body forces ||Br87| , pp^ . Rather than work with a complicated DBHF interaction, the 
basic point can be illustrated using the RMFT model of Section ^ for suitable choices of 
couplings and masses, simple scalar and vector exchange resemble quite closely the scalar and 
vector components of the DBHF interaction F. (We are not concerned here with numerous 
details like pi- and rho-exchange diagrams, the pion-exchange “tensor force”, etc. We will 
also set the nonlinear scalar couplings k, and A to zero to simplify the equations.) 

The basic idea is to expand the RMFT nuclear matter energy density of Eq. ( p.171) in 
powers of the Fermi momentum kp, and then to group terms together to isolate various 
powers of pa- One can then identify contributions from increasing powers of as aris¬ 
ing from two-body forces, three-body forces, etc., if one attempted to reproduce the same 
energy/nucleon with nonrelativistic potentials. We begin with the expansion of the self- 
consistency condition (|2.19|), which is given through order by ||Sef 


M* = M 
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(4.7) 


Here pp = as usual. 

Substitution of this result into Eq. (|2.17| ) and expansion through order k]} yields, for the 
energy/nucleon. 
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(4.8) 


The first term is the baryon rest mass, followed by the nonrelativistic Fermi-gas energy and 
the first few relativistic corrections, which are essentially negligible at equilibrium density. 
(These terms are the ones usually used to justify a nonrelativistic treatment of the nuclear 
matter problem ||lNe85|| .) The next two terms (proportional to p^) give the nonrelativistic 
limit of the potential energy coming from the vector and scalar mesons. The following term 
in brackets (with overall factor p^) is a relativistic correction to the scalar potential energy 
that arises from the Lorentz contraction factor in the scalar density, evaluated for nucleons 
of mass M. The final two terms (with overall factors of p\ and Pg) are also corrections to 
the scalar potential energy. These arise from the self-consistency condition ([4.7|) on M*; 
self-consistency implies M* < M, which increases the velocities of the nucleons and thus 
also increases their energies. These repulsive contributions are a signature of the velocity de¬ 
pendence inherent in a Lorentz scalar interaction. Nevertheless, the leading contributions to 
these terms (in powers of Pb) can be reproduced in a nonrelativistic calculation by including 
repulsive three-nucleon and four-nucleon potentials 
There are several relevant points to be noted: 


• Important “relativistic” effects in the RMFT of nuclear matter are equivalent to many- 
body forces in a nonrelativistic framework. 

• Although only the three- and four-body contributions are shown in Eq. ([4.8|) , all of 
the leading terms at each power of pe arising from the self-consistency condition are 
repulsive. Together with the Lorentz contraction contribution, they provide the new 
saturation mechanism discussed earlier. 


Since the scalar and vector mesons have relatively large masses (compared to the 
pion), the nonrelativistic many-body potentials used to generate these terms will be 
short-ranged ||Co95a . 


• With typical values for the RMFT parameters, one finds g^pB/m^M ^ 0.5 at equi¬ 
librium density. It follows that the three- and four-body repulsive terms in Eq. (|4.8|) 
are roughly 5 MeV and 2.5 MeV, respectively, which are significant on the scale of the 
nuclear matter binding energy. 


These relativistic (or many-body, if you prefer) contributions are inherent in the DBHF 
framework and allow a two-hole-line calculation to (essentially) reproduce the equi¬ 
librium point of nuclear matter, with an NN interaction fitted to phase shifts (see 
Fig. 10.13 in ||Ma89|| ). Although, to our knowledge, a relativistic three-hole-line calcu¬ 
lation of nuclear matter has never been performed, one might hope that the smaller 


38 







































depletion of the Fermi sea (as noted above) would lead to smaller three-body correc¬ 
tions than in the nonrelativistic case, implying that the two-hole-line results are fairly 
robust. Relativistic three-body-cluster effects remain to be investigated. 


5. EFFECTIVE FIELD THEORY (EFT) 


A. Introduction 


Quantum chromodynamics (QCD) is generally accepted as the underlying theory of 
the strong interaction. As we have seen, however, at energies relevant for most nuclear 
phenomena, hadrons are convenient and efficient degrees of freedom. In particular, meson- 
exchange models of the NN interaction accurately describe low-energy properties of the two- 
nucleon system, and relativistic mean-field theory based on QHD-I and QHD-II provides a 
realistic description of the bulk and single-particle properties of nuclei. 

Unfortunately, renormalizable QHD models have encountered difficulties due to large 
effects from loop integrals that incorporate the dynamics of the quantum vacuum [Co87, 
Pe87,Fu88,Fu89,We90,Li90a]. On the other hand, the “modern” approach to renormaliza¬ 
tion |[Le89| , po^,|Ge94| , |We95|| , which makes sense of effective, “cutoff” theoriesf^ with low- 
energy, composite degrees of freedom, provides an alternative. When composite degrees 
of freedom are used, the structure of the particles is described with increasing detail by 
including more and more nonrenormalizable interactions in a derivative expansion ||Le89| . 


Moreover, nonrenormalizable interaction terms between the boson fields allow us to describe 
the short-distance behavior of the underlying theory of QCD. So, based on the successful 
phenomenology we have seen so far, we would like to generalize the description using the 
modern viewpoint of effective field theory (EFT). 

The strategy underlying EFT relies on two basic observations. First, one argues that 
relativistic quantum field theory is simply the most general way to parametrize an S matrix 
(or other observables) consistent with analyticity, unitarity, causality, cluster decomposition, 

lWe95i 


and symmetries {e.g., Lorentz covariance, chiral symmetry, ...) 


With this view. 


there is no reason that relativistic quantum field theory should be reserved for “elementary” 
particles only. 

Second, one observes that in most problems in physics, the relevant phenomena are 
confined to a specific length scale, and thus it is not necessary to explicitly include dynamics 
at significantly shorter length scales ||Ge93|| . This implies that, at least formally, one can 
construct an effective field theory to be used at a given length scale by “integrating out” 
heavier degrees of freedom corresponding to shorter length scales; the effects of these heavier 
degrees of freedom will be implicitly contained in various coupling parameters in the low- 
energy, effective theory ||Le89| , I^95a|,Ma96 |. By fitting these parameters to experimental 
data, one can derive relationships between different observables within the dynamical regime 
of interest. 


^^Here “cutoff” means a regulator that maintains the appropriate symmetries. 
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Note that this strategy is opposite to the strategy of renormalizable held theory, where 
one argues that the fixed number of unknown parameters can be determined at any conve¬ 
nient length scale and then extrapolated to any other length scale using the equations of 
the “renormalization group”. Indeed, the effective theory is expected to contain numerous 
couplings of nonrenormalizable form when one integrates out the heavier degrees of free¬ 
dom. These nonrenormalizable couplings incorporate the “compositeness” of the low-energy 
degrees of freedom through their implicit dependence on short-distance physics. 

These considerations imply, of course, that the low-energy effective theory will generally 
contain an inhnite number of interaction terms, and thus one needs an organizing principle 
to make sensible calculations. First, one must hnd a suitable expansion parameter (or pa¬ 
rameters) that is small in the region of interest. Second, one assumes “naturalness”, which 
means that all of the unknown couplings in the theory, when written in appropriate dimen¬ 
sionless form (as discussed below), are of order unity. Thus one can estimate contributions 
from various terms by counting powers of the expansion parameter (s) and then truncate the 
lagrangian at the desired level of accuracy. 

Although the resulting framework is apparently structured less rigidly than the more 
familiar edihce of renormalizable quantum held theory, EFT nevertheless has its own “rules 
of the game”: First, one cannot simply limit calculations to the tree level; loops can and 
must be included, as this is the only way to correctly incorporate unitarity. Second, un¬ 
known parameters are to be determined either by explicitly integrating out the short-distance 
physics and “matching” the low-energy parameters to the results (if the underlying theory is 
tractable) or by htting to experiment (where, hopefully, the number of parameters is fewer 
than the number of data to be described). 

In the case of the nuclear many-body problem, it is still impossible to compute the desired 
low-energy parameters by working directly with QCD. Nevertheless, important constraints 
on the effective hadronic lagrangian are established by maintaining the symmetries of QCD. 
These include Lorentz invariance, parity conservation, isospin invariance, chiral symmetry, 
and electromagnetic gauge invariance. These symmetries constrain the lagrangian most di¬ 
rectly by restricting the form of possible interaction terms; if one is forced to £t unknown 
parameters to the data, one must include in the effective lagrangian all (non-redundant) 
terms that are consistent with the underlying symmetries. Moreover, the nature of the sym¬ 
metries may also dictate the appropriate low-energy degrees of freedom, or give relationships 
between some of the unknown parameters. 

Redundancy may arise because there is signihcant freedom in choosing the generalized 
coordinates of the effective lagrangian. In contrast to renormalizable theories, where there 
is (usually) a preferred choice of held coordinates in which the lagrangian is manifestly 
renormalizable, the choice of held variables in an ehective theory is motivated by the desire 
to make the description of the interactions as efficient as possible. (This is similar to the 
situation in classical lagrangian mechanics.) If a point transformation of the helds (subject 
to some mild constraints ||Co69|| ) allows one to eliminate certain interaction terms, they 
can be considered redundant ||Ba88| , [Ce^ , |Ba94b[] . The goal is then to hnd the best set of 
generalized coordinates (helds), so that the inevitable truncations are as accurate as possible. 

A well-known application of EFT is chiral perturbation theory (ChPT), in which one 
observes that the spontaneous breaking of chiral symmetry in QCD implies that Goldstone 
bosons (pions, ...) are the relevant low-energy degrees of freedom |[Le94||. Chiral symme- 
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try also implies that the pion interactions can be grouped order-by-order in the number of 
derivatives, so that there is a systematic expansion at low energies in powers of external 
momenta (and m^) ||Ga84|| . Moreover, because the loop expansion also proceeds in powers 
of momenta, one can systematically compute loop corrections. ChPT has been reasonably 
successful in describing scattering in the i? = 0 and B = 1 sectors of low-energy QCD 
Po92|,^ie93a)| . Studies of two- and many-nucleon systems have been initiated and are cur¬ 
rently under active investigation ||Ur92| , [Va9B| , |Ly93yUr96| , |Ka96|| . However, the prospects for 
extending ChPT to many-body calculations at hnite density are unclear at present. 

Thus we are motivated to consider alternatives. How can we exploit the ideas of EFT to 
develop a systematic, nonrenormalizable approach to nuclear structure? Let us enumerate 
the relevant concepts: 


1. Pions can be included by applying the framework of ChPT. 

2. Nucleon helds are necessary, since these are the observed fermionic degrees of freedom 
at low energies; that is, nucleons carry the conserved baryon number B. Moreover, 
nucleon compositeness is retained even at the “tree” level, if we include nonrenormal¬ 
izable interaction terms, as noted earlier. 

3. QCD constraints will be imposed through symmetries] all allowed (non-redundant) 
terms must be included. 


4. Since redehnitions of the held variables do not affect observables, we must strive for 
the most efficient variables and parametrizations of the interactions. 

5. We must identify suitable expansion parameters. Based on our earlier discussion, we 

know that the nuclear mean helds (or self-energies) <h = gsipo and W = are roughly 
several hundred MeV at ordinary densities. Thus we take as expansion parameters the 
ratios <h/M and W/M, which are small at normal densities, and also the ratios of 
gradients of the helds to M^, which are small in nuclei (|V<h|/M^ |VW|/M^ < 

0.1) |Fil95^ . If we then assume that the unknown coefficients are natural, we can 
truncate the ehective lagrangian. Of course, after the calculations are hnished and 
the parameters have been htted to empirical data, we must check that the naturalness 
criterion is satished. 


6 . What about non-Goldstone bosons, like a;, p, and our ehective scalar held 0? We 
argue that these are known to be useful degrees of freedom, since the mid-range NN 
interaction can be efficiently described by the exchange of these mesons ||Ma87|,[Ma89 


Moreover, the coupling constants in ChPT at 0{E^) in the meson sector can be repro¬ 
duced by a meson-resonance lagrangian applied at tree level, with the vector mesons 
playing the leading role ||Me88| , po89| , pilc89| , |Ec89a|| ; this is consistent with the well- 
known hypothesis of vector-meson dominance [^a69|| . Finally, we have already noted 
that two-pion exchange in the scalar-isoscalar channel can be efficiently simulated by 
a low-mass scalar held.^ Even though both the nucleon and scalar are composite ob- 


^^Recent results indicate that a low-mass scalar plays a similar role in ttN scattering |IV[e96|. 
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jects, recent work shows that the leading term [of 0 ($^/M)] in the interaction energy 
of a nucleon in a classical scalar held is model independent and determined solely by 
Lorentz covariance ||Wa95a|,[bi95| . 


7. An accurate description of pion interactions with few-body systems and nuclei requires 
that the A resonance be included. This can be done most efficiently by introducing the 
A as another effective degree of freedom ||Pe68|,pi8^ , Pe92^ , |We93|,P?a96|| . For simplicity, 


and because we are concerned primarily with bulk and isoscalar properties of nuclei, 
we will omit A interactions in the models discussed below. 


Thus our strategy will be to construct an EFT lagrangian containing nucleons, pions, 
and low-mass scalar and vector mesons. Heavier particles (with masses greater than roughly 
1 GeV) will be integrated out. 


B. Naive Dimensional Analysis 


There are still two important points to be addressed. First, we must understand how 
to extract the dimensional scales of each term in the lagrangian, so that the remaining 
dimensionless constants can be checked for naturalness. A naive dimensional analysis (NDA) 
for assigning a coefficient of the appropriate size to any term in the effective lagrangian has 
been proposed by Manohar and Georgi ||Ma84a|JGe93a|] . This allows for a determination of 
both the dimensional scales associated with each term and for the inclusion of an overall 
dimensionless constant that can be used to adjust the strength. The basic assumption of 
naturalness is that once the appropriate dimensional scales have been extracted, the overall 
dimensionless coefficients should all be of order unity. The NDA rules for a given term in 
the lagrangian density are: 


(1) Include a factor of I/Ztt for each strongly interacting held. 

(2) Assign an overall factor of 

(3) Multiply by factors of 1/M to achieve dimension (mass)^. 

(4) Include appropriate counting factors (such as l/n\ for 0"'). 


Here f-jr ~ 93MeV is the pion-decay constant, and the nucleon mass M is taken as the 
generic large-momentum cutoff scale, which characterizes the mass scale of physics beyond 
Goldstone bosons. In some cases ||Fu96b|| , it is more appropriate to use the (non-Goldstone) 
meson masses rather than M, but we usually do not distinguish here between the two. 

As noted by Georgi ||Ge93a|| , rule (1) simply assumes that the amplitude for producing 
any strongly interacting particle is proportional to the amplitude for emitting a Goldstone 
boson. This is a reasonable assumption, since /,r is the only natural scale. Thus, by dividing 
each held by /^r, we should arrive at a factor of 0(1). Rule (2) can be understood as an 
overall normalization factor that arises from the standard way of writing the mass terms of 
non-Goldstone bosons. For example, one may write the mass term of a scalar-isoscalar held 
cj){x) as 
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(5.1) 


1 n , n 1 .n , 

Jpj2 ’ 

where the scalar mass rris is treated as roughly the same size as M. By applying rule (1) 
and extracting the overall factor of the remaining ratios are of 0(1). Since all terms 

will have the same overall scale factor /^M^, higher-order terms or terms with gradients of 
helds will be suppressed by powers of 1/M relative to the leading mass terms, as a result 
of “integrating out” physics above the scale M. (A simple example is the low-momentum 
expansion of a tree-level propagator for a heavy meson of mass mn, which leads to terms with 
powers of 3“^/m\.) It is exactly because of these 1/M suppression factors and dimensional 
analysis that one arrives at rule (3). The origin of the combinatorial factors in rule (4) is 
discussed in |[Fu96a . 


Applying these rules to a generic term in the effective lagrangian involving the isoscalar 
helds and the nucleon held leads to (generalization to include the pion, rho, and photon is 
straightforward) |[hr9ti.l'u96a 


C 


1 I ( (t)\^ fVy (3 OI m^\P ^ ^ 


( 6 . 2 ) 


where is a baryon held, T is any Dirac matrix, derivatives are denoted generically by 3, 
and we have allowed for the possibility of chiral-symmetry-violating terms that contain the 
small parameter The product of all the dimensional factors then sets the scale in 

terms of the pion-decay constant and the nucleon mass M. The overall coupling constant 
g is dimensionless and of 0 ( 1 ) if naturalness holds. 

These scaling rules imply that a general potential for the scalar meson can be expanded 


as 


Ds = 


1 ^ ^ ^ 

2^3! M ^ 4! M2 ^ 


(5.3) 


Here we have included a factor of 1//^ for each power of cf); these factors are then eliminated 
in favor of Qs ~ M/ which is basically the Goldberger-Treiman relation ||Wa95|| . Factorial 
counting factors are also included, since the NDA rules are actually meant to apply to the 
tree-level scattering amplitude generated by the corresponding vertex ||We9U4F'u96a . 

The “naturalness” assumption states that after the dimensional factors and appropriate 


counting factors are extracted, the overall dimensionless coefficients [g in Eq. ( |5.2D and the 
K 3 , K 4 , ... in Eq. (|5.3|) ] should be of order unity. It should be clear, however, that the 
preceding arguments are not a proof of naturalness, since we know of no physical law that 
forbids large coefficients from appearing. Nevertheless, without such an assumption, it is 
basically impossible to construct an effective lagrangian with any predictive power.^ Until 


2 ^The assumption of renormalizability also leads to a finite number of parameters and well-defined 
predictions, but does so by imposing unnatural restrictions on the lagrangian, namely, that many 
parameters are identically zero in the absence of relevant symmetry arguments. 
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one can derive the effective hadronic lagrangian from QCD, the naturalness assumption must 
be checked by fitting to experimental data. 

The second point to be addressed is that we are actually fitting the QHD parameters 
based on calculations of finite-density observables (rather than scattering observables). One 
could, of course, fit the parameters to scattering observables, and then perform corresponding 
calculations of nuclear properties, but we must nevertheless understand how fitting parame¬ 
ters at the mean-field level is to be interpreted in the context of EFT. For this interpretation, 
we rely on the ideas of density functional theory, which we discuss in Section |^. 


C. Nonlinear Chiral Symmetry Revisited 


Before applying the ideas of EFT to the construction of a model lagrangian, we illustrate 
some basic features of nonlinear realizations of chiral symmetry by building on our earlier 
results. Return to the lagrangian of the chiral aoo model, which (after spontaneous symmetry 
breaking) is given in Eq. (|3.9|) . Recall that the nucleon mass M is related to the vacuum 
expectation value of the original scalar field (which we denote by sq) through 


M 

So — — . 

9tt 


(5.4) 


We now introduce an SU{2) matrix U defined in terms of the pion fields by 


U = exp f— T • tt) , 

^So ^ 


(6.6) 


as well as left- and right-handed baryon fields 


>(’1 = 5(1-75)'/', '/’R = 5(1+75)'/' . 


( 6 . 6 ) 


where V’ is the baryon held in Eq. (|3.9|) . Consider now the following phenomenological, 
nonlinear generalization C of Eq. 


C = i 


+ WvV^)il)i\ - 9 itSo(i -) f iPrU^l + '^PlU'^R 

J ^ So^ 

- \{dfj,ad^a) + \sl ti - V{U,d^,U]a) + +U^ -2) 




(5.7) 


Here V is a meson potential built from the indicated fields and their derivatives. For = 0, 
the new lagrangian is invariant under chiral S'[/(2)l x SU{2)r transformations of the form 


'/'H^fil/'R, U^LUB^, 


(5.8) 


this point, we free ourselves from the constraint of renormalizability. 
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where L and R are independent, global SU{2) matrices, as long as V is chosen to be invariant. 
(The a and helds are nnchanged.) 

It is now a simple exercise to show that in the limit sq —cxo, in which the nncleon mass 
M becomes very large, the two lagrangians are identical: 

C, ) sq ^ oo , (5.9) 


provided only that the potential in C is chosen to take the simple form 

V{U,d^U](T) = ^mla^ + 0(—) . 

Z ''So' 


( 6 . 10 ) 


Inspection of Eq. (p.9|) shows that the meson interactions in the auj model rednce to this form 
in the indicated limit. More generally, different choices for V lead to different descriptions 
of the spontaneons breaking of chiral symmetry, all of which preserve the SU{2)i^ x S'f/(2)R 
invariance of the lagrangian. 

To proceed farther with C of (|5.7|) , we implement the following clever change of variables, 
which leaves the particle content nnchanged [|D^: 

U = ^^, iVi = eVL, = (5.11) 


This change of variables takes the “sqnare root” of U and mixes the pions into the nncleons 
in a manner reminiscent of the chiral transformation in Eq. (|3.16| ). Straightforward algebra 
now shows that the part of C containing the fermion helds becomes (we leave the remainder 
of C in terms of U) 


C 


fermion 


= N 


+ iv^ + ig^Vf,) + 7 ^ 750 ^ - M + g^a 


N 


where 


-n, ^ ^ -ki'u - m') ■ 


(5.12) 


(5.13) 


^ Ml / ’ -f't 2 

The new form of the lagrangian is invariant nnder the following nonlinear chiral trans¬ 


formation Do92 


f(x) ^ L{{x)h\x) = h{x){{x)R^ , 
N{x) —>• h{x)N{x) , 


(6.14) 

(5.15) 


where h{x) is a local SU{2) isospin transformation that is dehned by Eq. (|5.14| ) and that 
depends on the pion held throngh This transformation implies 


—*• ha^h^ , —*• hv^h^ — ihd^fi^ . (5.16) 

By nsing the relation d^{h}h) = 0, one also obtains 

+ iv^)N ^ h[{d^ + IV^)N] , U ^ LUR} , (5.17) 

which means that the indicated derivative of the nncleon held transforms covariantly [that is, 
the same way as the nncleon held in Eq. (p.l5|) ], and the pion matrix U transforms globally. 
Thns the remaining parts of the lagrangian C remain nnchanged nnder the transformation 
(except for the symmetry-violating term proportional to m^). 

This new realization of the chiral symmetry has the following important properties: 
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• Parity is conserved. 

• The fermions (baryons) appear in isospin multiplets. 

• A baryon mass term is allowed in the lagrangian. 

• The Goldstone bosons (pions) arising from the spontaneous chiral symmetry breaking 
enter through the chiral matrices U and 

hence one reproduces all the soft-pion results implied by chiral symmetry. 


The pions are coupled to the nucleons with derivative couplings [Eqs. (5.12|) and (5.13 


• The isospin transformation h{x), which is an element of the unbroken SU{2)y subgroup 
of the full SU{2)i^ X S't/(2)R group, is local, it depends on x because the pion helds 
contained in ^ depend on x. 

• The scalar held a (which still appears in the theory) and vector held are chi¬ 
ral scalars. They can therefore be removed without destroying the invariance of the 
lagrangian, or alternatively, additional chiral scalars can be included. 


• One can rewrite the preceding results in conventional form by identifying 


"50 /tt • 


(5.18) 


In the elegant papers ||(Ja69| , |Co69|| , it is shown that any nonlinear realization of SU(2)2, x 
SU(2)j:j with these properties can be brought into the form of Eqs. (|5.14|) and (|5.15|) . 


D. A Nonlinear Chiral Model with Vector-Meson Dominance 


We now show how the ideas of EFT can be combined with the hadronic phenomenology 
discussed in Sections M by constructing a nonrenormalizable, ehective lagrangian that re¬ 
alizes chiral symmetry in a nonlinear fashion |[Fu95|jFu96a|| . Vector mesons are included in 
a manner consistent with both chiral symmetry and vector-meson dominance (VMD), and 
a light scalar held is introduced as before to simulate the exchange of two correlated pions 
between nucleons. Although the dynamics described by these non-Goldstone bosons could 
also be generated through pion-loop diagrams, the advantage is that we can now avoid (at 
least initially) the evaluation of complicated loop integrals in studying the many-body prob¬ 
lem. Moreover, the low-energy electromagnetic structure of the nucleon is described within 
the theory using VMD, so that ad hoc form factors are not needed, as in the calculations of 
Section p.B| . 

The general framework for nonlinear realizations of chiral symmetry is stated very com¬ 
pactly in the original work of Gallan, Goleman, Wess, and Zumino (GGWZ) ||Ga69|,[We96 


and here we paraphrase their discussion. (The results of the previous subsection illustrate a 
specihc example of these ideas.) We assume that G is a compact, connected, semisimple Lie 
group that has a continuous subgroup H. We denote by V) and Ai a complete, orthonormal 
set of generators of G, such that V sxe the generators of H. In some neighborhood of the 
identity of G, every group element g E G can be uniquely decomposed as 
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( 5 . 19 ) 


9 = 


^a-A ^(5-V 


where a ■ A = (3 - V = J2i(3iVi, and ae and Pi are real constants. 

A nonlinear realization of G that becomes a linear representation when g E H is given 
on the local held variables by the transformation 




where 


g^u^-A ^ gO^'.A ^u'-V 


!/>' = , 


(5.20) 


(5.21) 


uj' = uj'{uj;g), and u' = u'{(jj]g). Here uj represents the Goldstone boson helds, p denotes 
the other helds in the theory, and uj' is generally a nonlinear fnnction of uj. D{h), with 
h E H, is Si linear, nnitary representation of the unbroken subgroup if, which is assumed to 
be written in fully reduced form. In ||Co69|| it is proved that any nonlinear realization of G 
that is linear on H can be brought into the preceding form through a suitable redehnition 
of held variables. 

Unlike a linear realization of the symmetry, in which all particles must be assigned to a 
representation of the full group G, the non-Goldstone helds -0 in a nonlinear realization are 
in representations of the unbroken subgroup H only, just as in the example in the previous 
subsection (where the unbroken group is the SU{2)y of isospin). Information about the 
structure of the full group G is encoded in the Goldstone helds uj and the local transfor¬ 
mation with e“ G H. Thus a nonlinear realization eliminates the necessity of 

assigning particles to chiral multiplets, as in the linear a model and its extensions |[Se92b|| . 
These assignments are additional dynamical input that go beyond the assumption of group 


invariance alone 


We now proceed to construct our model following the strategy outlined at the end of 


Section b-A| . We include helds for the nucleon, pion, rho meson, and omega meson, together 
with an ehective scalar-isoscalar held to simulate two-pion exchange; the A will be omitted 
here for brevity, but its inclusion is discussed in 


Ta96|] . The Goldstone pion helds 7r“ 


Uc 


with a = 1, 2, and 3, form an isovector, which can be considered as the phase of a chiral 
rotation matrix: 


f (x) = exp(i7r(x)//^) , 7r{x) = n{x) . 


(5.22) 


Here the r“ are Pauli matrices and ~ 93MeV is the pion decay constant. The isospinor 
nucleon held is represented by a column matrix 


N{x) = 


p{x) 

n{x) 


(5.23) 


where p{x) and n{x) are the proton and neutron helds respectively. The rho helds p“(x) 
also form an isovector, and we use the notation | r • p^{x). 

A nonlinear realization of the chiral group SU{2)i^ x SU{2)ji is now dehned such that 
for arbitrary global matrices L E SU{2)i^ and R E S'17(2)r, there is a mapping 
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(5.24) 


i 0 B : (^, p N) p' N') . 


Because of the parity operation V, which produces the transformation 

V: Li -7r“(t,x) —>-7r“(t,-x) , f(t,x) — 

the chiral mapping (^.21| ) can be written as [compare Eqs. ( b.l4|) and (|5.15|) ] 

i\x) = Li{x)h\x) = h{x)i{x)Ii^ , 
p'J^x) = h{x)pj^x)h\x) , 

N\x) = h{x)N{x) . 


(5.25) 

(5.26) 

(5.27) 

(5.28) 


The second equality in Eq. (|5.26|) dehnes h{x) as a function of L, R, and the local pion 
helds: h{x) = h{L, R,7r{x)). It follows from Eq. (^.26|) that h{x) is invariant under the 
parity operation ( |5.25|) , that is, 


h{x) e SU{2)y , (5.29) 

with SU{2)y the unbroken vector subgroup of SU{2\ x S't/(2)R.[^ Equations ( |5.27|) and 
( b.28| ) ensure that the rho and nucleon helds transform linearly under SU{2)y in accordance 
with their isospins. Note that the matrix h{x) becomes constant only when L = R, in which 
case g G H = SU{2)y and h = L = R. The isoscalar helds V^{x) and (f){x) are chiral scalars 
and are unahected by both chiral and isospin transformations. 

For discussing purely pionic interactions, it is convenient to dehne the matrix [compare 
Eqs. (1^) and (|5.11|) ] 


U{x)=i (x) = exp(2i7r(a;)//^) 


(5.30) 


since the transformation law (|5.26|) then implies 

U{x) —^ U\x) = LU{x)R^ , 


(5.31) 


so that U{x) always transforms globally [see Eq. ( b.l7|) ]. Thus derivatives of U{x) transform 
the same way as U{x), and chirally invariant interactions involving pions alone can be 
constructed from products of U{x), U\x), and their derivatives. As is well known, these 
terms can be organized according to the number of derivatives, with the lowest-order term 


can express h in terms of the matrices L, R, and U [see Eqs. (5.30) and (5.31)] as 
h{x) = U'\x) L^Ufx) = \JRU^{x)L^ L^U{x). Given the decomposition fCoGSl] L = 


exp(iQ; • T)exp(i/3 • r), R = exp(— ia • T)exp(i/3 • r), and h = exp(i 7 • r), with a, (3, and 7 
real, the infinitesimal expansion of h is determined by 7 = (3 — {a x 7r)j2R + 0(q;^,/3^ 


^2, TT^i 


| Me93a| , pe^ 


48 













































A = ^ tr ([7r,a^7r]^) + 


(5.32) 


C 2 determines all multipion scattering amplitudes to second order in external momenta in 
terms of the single constant /t,. Terms with more derivatives (£ 4 , Cq, .. .) can be used to 
describe pion dynamics within the framework of ChPT ||Ga84|| . We will not need more than 
£2 in the model studied here. 

For describing the interactions of pions with other particles, U{x) is not convenient, be¬ 
cause other fields transform with the local h{x) of the unbroken isovector subgroup SU{2)y. 
It follows from the transformation laws given earlier, that interaction terms that are invari¬ 
ant under local isospin rotations will be invariant under global transformations of the full 
group S'[/(2)l X SU{2)ii. Thus, to form chirally invariant interactions involving pions and 
other fields, we need functions of the pion field that transform with h{x) only. 

The desired functions involving one derivative of the pion held can be written as [compare 
Eq. ifl^] 




Vf, = + =vl , 


(5.33) 

(5.34) 


where the hermiticity follows from = 0 = Under parity transformations, we 

have 


V : ^ , u^(f,x) —, (5.35) 

so that is an axial vector and is a polar vector. To leading order in derivatives, one 

hnds 

In 

+ 

Moreover, under a chiral transformation, Eq. ( b.26| ) implies 

^ , (5.38) 

^ v'fj, = hv^h^ — ihdfj,h} = hv^h} -|- i{d^h)h} . (5.39) 

Thus transforms homogeneously under the local SU (2)v group and can be interpreted 
as a covariant derivative of the pion-held matrix ^(a:). In contrast, the inhomogeneous 
transformation law for resembles that of a gauge held, so that allows us to construct 
chirally covariant derivatives of the other helds. For example, it is straightforward to verify 
that the covariant derivatives 

Df,N = + iy^)N , = df,p^ + i[v^, pj (5.40) 

transform homogeneously with h{x) under the full group: 


(5.36) 

(5.37) 
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(5.41) 


(D,N)' = h(D,N) , (D,py = . 

The covariant derivative of p can be used to construct the covariant held tensor 




(5.42) 


The antisymmetric combination of derivatives implies that the timelike components of 
the rho held have no conjugate momenta and are thus determined by equations of constraint, 
as appropriate for a massive vector held with three dynamical degrees of freedom. The hnal 
term in Eq. (|5.42| ) has the usual form for a non-abelian vector held and enables the p meson 

IKeSQI 


to couple to a conserved isovector current 
for the pion held as 


We can also construct a covariant tensor 




= SfVu - S^t!. + »[!'„. aj = a„] , 


(5.43) 


which transforms homogeneously with h, and which will allow us to produce an invariant 
piTTi coupling through an interaction of the form tr {p^^v^'^). 

Before exhibiting the lagrangian for the model, we consider electromagnetic interactions. 
As discussed in |[Fu96a|| , these can be included straightforwardly by dehning appropriate 
charge operators for the particles and by modifying the preceding covariant derivatives so 
that they remain covariant under the local U(l) transformations of electromagnetism. We 
will not present the intermediate steps here and simply show the contributions to the model 
lagrangian below. The electromagnetic interactions induce small violations of both the 
isospin and chiral symmetries. 

To write a general effective lagrangian, we need an organizational scheme for the in¬ 
teraction terms. We organize the lagrangian in increasing powers of the helds and their 
derivatives, as motivated by the principles discussed earlier. We assign to each interaction 
term an index 


n 


u = d + - + b , 


(5.44) 


where d is the number of derivatives, n is the number of nucleon helds, and b is the number 
of non-Goldstone boson helds in the interaction term. The hrst two terms in Eq. (|5.44|) are 


suggested by Weinberg’s work |[We90a|| . Derivatives on the nucleon helds are not counted 


in d because they will generally introduce powers of the nucleon mass M, which will not 
lead to small expansion parameters. The last term is a generalization that arises because a 
non-Goldstone boson couples to two nucleon helds. Equation ( b.44| ) is also consistent with 
hnite-density applications when the density is not too much higher than the nuclear-matter 
equilibrium density, as we will see below. 

The ehective lagrangian for the full model can be written as 


£ = £n + , (5.45) 

where each term will be truncated by considering the various values of n, as dehned above. 
Only non-redundant terms will be exhibited, and we will return to consider the redundancy 
problem after describing the lagrangian. 
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The part of the effective lagrangian involving nucleons can be written through order 
z/ = 3 as 


/:N(a:) = + iv^ + igpp^ + ig^Vp) + 5'a7^75S^ - M + g^cj) 

-IfiL No ^ NVp^af^'^N + ■ • • , 

AM - - AM - - 


N 


(5.46) 


where — d'^V^. Note that both vector and tensor couplings are included for 

the p and uj mesons, together with a Yukawa coupling for the effective scalar held 0; we 
explain the motivation for a simple Yukawa coupling below. The ellipsis represents terms 
involving ttN interactions that are not needed in the following discussion, as well as terms 
with derivatives on the nucleon held, which will be considered below. Four-nucleon contact 
terms are not included, since they can be represented by appropriate powers of other meson 
helds |[Fu96a|| , as discussed shortly. Additional terms with v = A are either redundant or 
tiny. 


We emphasize that because the nucleon held obeys the transformation law (p.28|) , the 
mass term in Eq. (^.46 ) is chirally invariant. Thus the nature of the spontaneous chiral 
symmetry breaking is outside the realm of our ehective theory, unlike a linear realization, 
where the symmetry breaking is implemented by the helds in the model [see Eq. ( p.8|) ]. In 
a linear model, the mechanism that generates the free-space mass M is also responsible for 
shifting M M* at hnite density, which leads to the problems discussed in Section ^ A 
nonlinear realization of the symmetry allows these two aspects of the nucleon mass to be 
treated independently. 

The mesonic part of the lagrangian is also organized in powers of z/. Keeping terms up 
to order z/ = 4, we hnd^ 


= j/pr(8„ra'‘C/t) + 1^2/2 + pt _ 2) + 1(^1 + 

- ) - l(l + “ 2 ^)UW'“' - 9^,^ tr (p vn 


-uy- ' A^ M 

,2 a2 


+ i(i + + AM(v,v>‘f 


M 2 M2 
9s(t) 


1 

4!' 


, A , 9s ( t >\ 2. ( 2,2/1 K3 5-80 , K4 


2a2 


M2 


(6.47) 


where we have included a small chiral-symmetry-violating term involving m^. Apart from 
conventional dehnitions of some couplings {gs and g^) and the masses, the parameters are 
dehned so that they are of order unity according to the naive dimensional analysis discussed 
earlier. This hypothesis will be tested by htting the parameters to nuclear properties. More¬ 
over, since the expectation value of the p held is typically an order of magnitude smaller 


2^Note that the leading-order coupling of the scalar-isoscalar field to pions has the form 
0tr 
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than that of the uj held, we have retained nonlinear p conplings only throngh order u = 3. 
Note that the ai and 02 terms involve derivatives of the meson helds and have u = 5, bnt 
these give contribntions to the nnclear snrface energy that are nnmerically of the same mag- 
nitnde as the qnartic scalar term, so we have retained them. Thus numerical factors such as 
1/n!, which are cancelled in scattering amplitudes, are relevant in deciding the importance 
of contributions to the energy. 

Higher-order self-couplings and derivatives involving meson helds alone [e.g., 0®, or 
or should be numerically small unless their coefficients are “unnat¬ 

urally” large. We show below that the parameters that have been retained indeed exhibit 
naturalness, so that the omission of these terms is justihed at the level of accuracy we can 
expect in comparisons with observables in hnite nuclei. 

Finally, the electromagnetic interactions are described by 






- NX'a^^^N - 

4M ^ ^ 


- - - - 2M2 


^2eflA>‘ tr (v^n) - tr 

zg^ o 


+ i3^^^Ts)Nd,F^^ 

+ ... 


(5.48) 


where is the electromagnetic potential and is now the usual held tensor. The 
lagrangian £em is invariant under the U{1) group of electromagnetism, and the resulting 
current [see Eq. ( b.3|) , below] is conserved, at least to 0(e).|^ The composite structure of the 
nucleon is included here through an anomalous moment [see Eq. (|6.1j) , below] and through 
terms that will generate a dependence in both the isoscalar and isovector 

electromagnetic form factors. Moreover, the coupling between the massive vector mesons and 
the photon generate contributions to the nucleon form factor in accord with vector-meson 
dominance. The end result is momentum dependence that resembles the empirical “dipole” 
form. Similar contributions will arise in the pion form factor due to the pnn coupling in 
Eq. (|5.47|) . We will return to the electromagnetic structure of the nucleon in this model in 
Section p. 

Turning now to “redundant” terms that have been omitted in the preceding equations, 
we emphasize that there is considerable freedom in the choice of generalized coordinates 
(helds) for the lagrangian. It is known that a wide class of point transformations of the 
helds do not change the on-shell scattering amplitudes ||Co69|| .p^ Thus the relevant question 
is which choice of coordinates leads to the most practical and accurate truncation scheme, 
and this is currently under active investigation. 

Some redundant terms that we have omitted through our choice of held variables are: 


1. Contact terms involving nucleon helds beyond bilinear order [e.g.^ {NN){ NN) or 
{ N'j^N)[ N'^fj^N)]. We observe that if one constructs the hamiltonian from the la- 


^^The fully U(l)-invariant lagrangian is discussed in |Fu96a|. 

assume that this is also true for finite-density observables, but we know of no proof. 
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grangian given above and then eliminates the meson helds using their equations of 
motion, contact terms involving products of nucleon bilinears will arise. Conversely, 
contact terms included originally in the hamiltonian can be eliminated in favor of the 
scalar and vector helds, if we allow products of helds to all orders. The issue then 
becomes one of efficiency, since one will always have to truncate in practice. Although 
both ways of representing these nonlinear interactions may turn out to be practical 
||Bo77| , pdM| , |N i92| , [hu96|| , the hts to nuclei presented below show that nonlinear meson 
interactions lead to an efficient and natural truncation. 

2. More complicated meson-nucleon couplings [e.g., gs{(j)) NN(j)]] these are often mo¬ 
tivated by the claim that they are necessary to incorporate the compositeness of 
the nucleon. Here we rely on the freedom to redehne the helds to rewrite compli¬ 
cated couplings in simple Yukawa form: rather than work with gs{4>) NN(j), where 
5fs(0) = 5's(l + Ci0 -|- C20^ + •••)) could dehne a new scalar held 0 by ^s0 = 5 's( 0 ) 0 - 
Inversion of this relation and substitution into the lagrangian would result in additional 
nonlinear interactions in powers of 0 that have the same form as those that have al¬ 
ready been included. Note that this procedure can actually be used on all possible 
scalar-isoscalar held combinations, so the following couplings are all redundant: 

FY0^ FY0^ NNd‘^4>, NNV^V^ . 


A similar observation holds for the vector-isoscalar and vector-isovector couplings. 
This discussion illustrates the important point that nucleon compositeness can be 
incorporated through nonlinear meson interactions. 


3. 


Couplings involving higher derivatives of the nucleon held [e.g., ( Nd‘^N)(t), NiV^d^N■, 
or (These are the most problematic, since derivatives acting on 
the nucleon held produce factors of the nucleon energies and since Ei/M k, 1 , 
this would spoil the expansion and truncation procedure outlined above. Fortunately, 
through partial integration, redehnition of the baryon held, and the use of the equations 
of motion in the construction of the hamiltonian, these terms can be recast in the form 
of the terms we have retained or can be shown to produce terms that should give only 
small contributions. Further discussion is contained in |[Ba88|,p'e91|,|Fu96a|| and in §7.7 
of ||We95|| . 


The mean-held equations and energy density resulting from the lagrangian (|5.45|) can be 
derived straightforwardly. For symmetric nuclear matter (7 = 4), one hnds through order 
z/ = 4 


^^Higher-derivative terms should of course be written in terms of covariant derivatives [see 
Eq. ( |5.40D ]. The contributions from the meson fields to these terms have already been classified; 
here we focus on the gradient pieces. 
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(5.49) 




L / SI/ 


$ r}2 


mtW'^ 


1/1 ^ ^3 4 > ^ fi :4 < 1 >^ 
9s 

1 


, 2^2 


rCoW" , 


2c/^v- ■ -M 2 M2/ '' 4!^2 

where $ = (7s0o and Ih = g^Vo are the scaled helds dehned earlier. One can also compute 


the bulk symmetry-energy coefficient ||Se 86 


04 = 


9p 


127r2m*2 ^ 


K + - 


kl 


6 Jkl + M *2 


(5.50) 


where the effective rho mass m* is defined by 


m* = m Al + gp^/M) 


(5.51) 


A comparison with Eq. ( p.l7 ) shows that the nuclear matter energy has been generalized to 
include additional nonlinearities that are not allowed in the renormalizable model QHD-I. 
The fields $ and W are again determined by extremization. 

The Dirac-Hartree equations for finite nuclei can also be derived using the procedures 
in Section |2.B| . The resulting equations are lengthy and will not be reproduced here; the 
interested reader is referred to ||Fu96a|| for details. One important result is that due to the 


additional nonrenormalizable interactions between the nucleon and the electromagnetic field, 
and also due to vector-meson dominance, the computed nuclear charge density automatically 
contains the effects of nucleon structure, and it is unnecessary to introduce an ad hoc form 
factor. 

It is clear from the discussion in Section that the present model has more than enough 
parameters to give an accurate reproduction of nuclear properties. The more important 


question is whether the parameters fitted to nuclei are natural. In |[Fu96a|| , the parameters 
were determined by calculating a set of observables nucleus and by adjusting 

the parameters to minimize the generalized defined by [|Ni92 


x^ = EE 

i X 


rAh) 

exp 


A 


(bn 

th 




(5.52) 


“exp” 


where i runs over the set of nuclei, A runs over the set of observables, the subscript 
indicates the experimental value of the observable, and Wx'^ are the relative weights. The 
weights are chosen to be the relative accuracy expected for the given observable in a good 
fit.|^ The nuclei chosen were ^°Ca, ^®Ca, ®®Sr, and 20®Pb. 

A total of 29 observables and their relative weights were taken as follows: 


2 ®In practice, a reasonable range of weights was tested, and the qualitative conclusions discussed 
below were always reproduced. Some of the considerations relevant in choosing the weights are 
discussed in Section |^. 


54 






































The binding energies per nucleon e/B, with a relative weight of 0.15% 

1 . 

The rms charge radii with a relative weight of 0 . 2 % 


The d.m.s. radii i?dms, with a relative weight of 0.15% 

The spin-orbit splittings AE^o of the least-bound protons and neutrons, with a relative 
weight of 5% for 15% for ^°®Pb, 25% for '^'’Ca and ^®Ca, and 50% for ®^Sr 

The proton energy T^p(lh 9 / 2 ) and the proton level splitting Ep{ 2 d 3 / 2 ) — Ep{lhii/ 2 ) in 
208pb^ with relative weights of 5% and 25%, respectively 

The surface-energy and symmetry-energy deviation coefficients 6 a 2 and ( 5 a 4 , each with 
a weight of 0.08. 


The so-called diffract ion-minimum-sharp (d.m.s.) radius of a nucleus is dehned to be ||hr82 


i?dms = 4.493/g(') , 


(5.53) 


where Q^q'^ is the three-momentum transfer at the hrst zero of the nuclear charge form factor 
F{Q) = -^chg(ci) with Q = |q|. The surface-energy and symmetry-energy deviation coeffi¬ 
cients 6 a 2 and ^04 are defined by htting the difference between experimental and calculated 
binding energies 6 ei = (ej)exp - (ei)th according to |3 


Sei = Sai Ai - 5a2 ("jy. _ jAi . 


(5.54) 


Here Ni and Zi are the number of neutrons and protons in the nucleus and Ai = Ni + Zi. 
(An exact £t to the energies would have Sai = 6 a 2 = = 0.) The deviations 6 a 2 and 

are included as separate terms in Eq. (|5.52| ) in the form [dai/WsaiY, with Wsai = 0.08. The 
motivation for this choice of observables is discussed more fully in ||Fu96a . 


The nucleon, a;, and p masses are taken to have their experimental values: M = 939 MeV, 
rUv = 782 MeV, and = 770 MeV. (Including the heavy meson masses as free parameters 
produces only minor changes in the hts.) The anomalous magnetic moments of the nucleon 
are hxed at Ap = 1.793 and = —1.913, and = 5.01 is chosen to reproduce the 
experimental partial width r(p° —> e’''e“) = 6.8 keV. The empirical free-space charge radii 
of the nucleon are used to £x /?s, and fp by solving Eqs. ( p. 8 |) and (|6.9|) , below. The 

remaining thirteen parameters Ps, S'v, 5 'p, hi) h 2 ) Vpi ^ 3 ) ^ 4 ) Co, /v, cq, and 0:2 for the 

V = A parametrization are then obtained by optimization of the generalized 

In Table ||, we show two parameter sets (G1 and G2) obtained from fits with roughly 
equal accuracy when all terms through order z/ = 4 are retained. The parameters have been 
displayed in such a way that they should all be of order unity according to NDA and the 


^®Note that phenomenological surface-energy and symmetry-energy coefficients are not used, so 
there is no direct input from nuclear matter to the fitting procedure. 
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TABLE II. Parameter sets from fits to finite nuclei, as described in the text. Note that sets 
W1 and Q1 include the same interaction terms as sets L2 and NLC in Table |. 


V 

VFl 

Cl 

Q1 

Q2 

G\ 

G2 

rUs/M 

2 

0.60305 

0.53874 

0.53735 

0.54268 

0.53963 

0.55410 

ffs/ 4 vr 

2 

0.93797 

0.77756 

0.81024 

0.78661 

0.78532 

0.83522 

5 v/ 47 r 

2 

1.13652 

0.98486 

1.02125 

0.97202 

0.96512 

1.01560 

Qp/An 

2 

0.77787 

0.65053 

0.70261 

0.68096 

0.69844 

0.75467 

Vi 

3 


0.29577 



0.07060 

0.64992 

ri2 

4 





-0.96161 

0.10975 

K3 

3 


1.6698 

1.6582 

1.7424 

2.2067 

3.2467 

K4. 

4 



-6.6045 

-8.4836 

-10.090 

0.63152 

Co 

4 




-1.7750 

3.5249 

2.6416 

Vp 

4 





-0.2722 

0.3901 

Ol 

5 





1.8549 

1.7234 

02 

5 





1.7880 

-1.5798 

/v/4 

3 





0.1079 

0.1734 

/p/4 

3 

0.9332 

1.1159 

1.0332 

1.0660 

1.0393 

0.9619 

/3s 

4 

-0.38482 

-0.01915 

-0.10689 

0.01181 

0.02844 

-0.09328 

Pv 

4 

-0.54618 

-0.07120 

-0.26545 

-0.18470 

-0.24992 

-0.45964 


naturalness assumption. This is seen to be the case.^ Most importantly, it is found that 
the accuracy of the £t and the contributions to the nuclear matter energy/particle are not 
driven by the last terms retained, as illustrated in Fig. 13. This result was further checked 


Fu96a| by including the u = 5 interactions 


n-.'s — — 


5! M3 


I ^ 3^3 
+ 777 % 


+ —Cl —Q^(V 

3!'''^M3 2 ^ ^ ’ 


(5.55) 


which do not improve the hts to the data; these contributions are essentially negligible unless 
the coefficients are unnaturally large, and no indication for such large parameters was found. 
Thus we conclude that NDA and the naturalness assumption are valid when applied to finite 
nuclei, and that the truncation procedure defined above is practical, at least for moderate 
densities. Moreover, although the parameters were obtained from a fit to a specific set of 
nuclei, one can now extrapolate to study other features of nuclear structure, such as nuclear 
deformations, isotope shifts in charge radii, etc. 

Also shown in Table are the results of fits with fewer parameters, which were obtained 
to verify that the best possible accuracy for this set of input data has been achieved by 
keeping terms through order z/ = 4. Roughly speaking, set W1 had ~ 1700, set Cl had 

~ 400, sets Q1 and Q2 had 100, while sets G1 and G2 achieved y^ 50. Including 


3°Natural parameters have also been obtained in a “point-coupling” model that describes the NN 
interaction through contact terms |Ni92,Fr96|. 
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FIG. 8. Contributions to the energy per particle in nuclear matter for parameter sets G1 and 
G2 from the n**^-order terms of the form where n = £ + m. The boxes are terms with £ = 0, 

the circles are terms with m = 0, and absolute values are shown. Results from set G1 are open 


and those from G2 are filled. The crosses are estimates based on Eq. (^). The arrow indicates 
the total binding energy cq = 16.1 MeV. 


the u = 5 parameters in Eq. (|5.55|) improved only slightly (< 2 units). Thus keeping 
terms through order z/ = 4 is essentially the best one can do, and in fact, the parameters 
are already underdetermined at this level, as is evident by the differences between sets G1 
and G2. 

Nuclear matter properties predicted by these parameter sets are given in Table fT|. Ob¬ 
serve that sets G1 and G2 yield similar results in spite of the differences in the parameters, 
implying that the nuclear matter properties are better determined than the parameters 
themselves. Note also that all sets including parameters through orders p = 3 or p = A pre¬ 
dict M* Ri 0.6 M at equilibrium, in agreement with our discussion in Section Although 
this result was obtained when spin-orbit information was included as input, fits without such 
information leads to similar values for M* |[Fu96aj|, as was first shown in |Prte86 . 


E. The Quantum Vacuum in QHD-I 

In any consistent relativistic field theory, one must ultimately consider loop diagrams. 
These contributions are an integral part of a fully relativistic description of nuclear structure, 
and as described in Section it is impossible to construct a meaningful nuclear response or 
consistent nuclear currents without including the negative-energy baryon states. Although 
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TABLE III. Nuclear matter equilibrium properties for sets from Table || and for the 
point-coupling model of Ref. [ Ni92[ (set PC). Values are given for the binding energy per nu¬ 
cleon (in MeV), the Fermi momentum fcp (in fm“^), the compression modulus K (in MeV), the 
bulk symmetry energy coefficient 04 (in MeV), M*/M, and 5 vPb (in MeV) at equilibrium. 


Set 

E/B - M 


K 

04 

M*M 

5vVb 

W1 

-16.46 

1.279 

569 

40.9 

0.532 

363 

Cl 

-16.19 

1.293 

304 

32.0 

0.657 

255 

Q1 

-16.10 

1.299 

242 

36.4 

0.597 

306 

Q2 

-16.13 

1.303 

279 

35.2 

0.614 

292 

G1 

-16.14 

1.314 

215 

38.5 

0.634 

274 

G2 

-16.07 

1.315 

215 

36.4 

0.664 

248 

PC 

-16.13 

1.299 

264 

37.0 

0.575 

322 


the MFT ground state is causal and consistent with Lorentz covariance and thermodynamics 
by itself, it is natural to ask about the role of contributions from the hlled Dirac sea. This is 
one of the motivations for constructing QHD-I as a renormalizable theory |[Wa74|J5e86|| . Our 
goal in this subsection is to determine if the simplest evaluation of these effects in QHD-I 
produces results that are consistent with NDA and naturalness. 

In Section we studied the consequences of the mean-field hamiltonian of Eq. (p.ll|) and 
its generalization to hnite nuclei. Let us now return to infinite nuclear matter and include 
the contribution from 6H in Eq. (|2.12|) . The inclusion of this term defines the so-called 
relativistic Hartree approximation or RHA. (This is also often called the one-baryon-loop 
approximation.) 

An inspection of 6H reveals that, even with the indicated vacuum subtraction, the sum 
still diverges. Since QHD-I is a renormalizable model, however, the sum can be rendered 
hnite by including counterterms in the lagrangian (|2.7|) . These counterterms also appear in 
the hamiltonian, and they can be grouped with 6H, resulting in a correction to the energy 
density of the form 


AS{M*) = -1 y [(k^ + M*2)V2 _ (k2 ^ ^2)1/21 _ y ^ 


^0 


(5.56) 


The counterterms enter as a quartic polynomial in 0o, and the (inhnite) coefficients are 
determined by specifying appropriate renormalization conditions on the energy. Following 
(Jh77|| and |[Se86|JWa95|| , we will choose the counterterms to cancel the hrst four powers 
of 00 appearing in the expansion of the inhnite sum. This is equivalent to dehning the 
renormalized parameters k and A to be zero. Although this procedure is not unique (and is 
also unnatural), it minimizes the contributions from this vacuum correction, and it is easy 
to verify that only the hrst four terms in this expansion produce divergent results. The 
divergences can be dehned by converting the sum to an integral and then by regularizing 
dimensionally ||Ch77| , |Se^ , |Wa95|] . 

After removing the divergences with the counterterms, the remaining terms are hnite, 
and one hnds (for spin-isospin degeneracy 7 = 4) 























AS{M*) = 


^ {M*Mn(M7M) + - M*) - -M^{M - M*f 


471-2 


+ —M(M - M*f - —(M - M*Y] 


$6 


$7 


^ 30M6 ^ 105M7 


4!(n-5)! 


n! 


M" 


}. 


(5.57) 

(5.58) 


where M* = M — q^c^q = M —$. A£^ is the hnite shift in the baryon zero-point energy that 
occurs at hnite density and is analogous to the “Casimir energy” that arises in quantum 
electrodynamics. Just as in the MFT, M* is determined at each pa by minimization, which 
produces the one-loop (RHA) self-consistency condition [compare Eq. (|2.19|) ] 


^2 1 

M* = M - ^ ps + - M) 

mi mi vr^ i- 


rUg TT 


^M{M* - Mf - y (M* 


Mf] . 


(5.59) 


Note that the solution to this equation contains all orders in the coupling g^. 

To discuss the size of the one-loop vacuum correction, we apply the NDA. Based on the 
scaling rules discussed above, a term of O( 0 o) should be scaled as 


M 2 


(5.60) 


and if this contribution is natural, any residual overall constant should be of order unity. 
However, if we perform a similar scaling on the leading term in Eq. (|5.58|) , we hnd 


M" 


5M® 


4 m2 


= 96 


M 2 

Wl 


(5.61) 


where we used 47r/,r ~ M and Qs ~ M/fn- Thus the one-baryon-loop contribution to 
the vacuum energy in QHD-I is roughly two orders of magnitude larger than naturalness 
requires. It is not hard to show from Eq. ( b.58|) that all higher powers of $ contain essentially 
the same large overall factor. 

Similar behavior occurs in the linear sigma model. If we rewrite the coefficients in 
Eq. (p.lO ) in terms of the nonlinear parameters in Eq. ( b.3|) , we hnd (in the chiral limit) 
K 3 = —K 4 = —3, so that the nonlinear parameters are natural at the mean-held level. 
However, if one includes the one-baryon-loop vacuum corrections, renormalized in a fashion 
that preserves the chiral symmetry |[[Vfa82| , [Fu93|| , one hnds unnatural corrections to the cubic 
and quartic couplings: Ak^ = 2M‘^k. 20, Afi :4 = —80. The quintic and 


higher corrections are exactly the same as in Eq. (|5.58|) . Thus the one-baryon-loop vacuum 
contributions again produce unnatural coefficients. 

These unnatural coefficients generate correspondingly large corrections to the MFT. If we 
consider QHD-I and adjust the model parameters from their MFT values to reproduce the 
desired nuclear matter properties in the RHA (equilibrium at = 1.30 fm~^ with a binding 
energy of 15.75 MeV), the baryon ehective mass at equilibrium becomes M*/M pe 0.73. This 
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FIG. 9. Contributions to the scalar potential per particle in nuclear matter from the n'''^-order 
terms of the form for the RHA model. The crosses are estimates based on Eq. ([5.21) . The arrow 
indicates the total binding energy cq = 15.75 MeV. 

translates into a change in the scalar potential <h from 430 MeV in the MFT to 250 MeV in 
the RHA, which is a large effect. In fact, the 0(<h®) term in the energy density forces <I> to 
signihcantly lower valnes, and the contribntions to the energy/nncleon from AS are mnch 
larger than what one wonld expect from this term in a natnral model, like the one discnssed 
in the preceding snbsection. This is illnstrated in Fig. P, where the RHA 0(<F^) contribntion 
for M*/M = 0.6 is as large as a typical 0(<I>^) contribntion in a natnral model. 

If one accepts the assnmption of natnralness, the conclnsion is that the treatment of 
the qnantnm vacnnm at the one-baryon-loop level in the renormalizable model QHD-I is, 
at best, inadeqnate.|^ Althongh higher-order corrections (some of which are discnssed in 
Section might rednce the size of the one-loop terms and nltimately yield a natnral size 
for the vacnnm contribntions, this can only occnr throngh sensitive cancellations against the 
one-loop terms. 

In contrast, in an effective held theory, the presence of nonrenormalizable conplings 
will generally canse the snm in Eq. (^.56|) to have divergences at all orders in 0O) and 




X natural EFT 
o RHA MVM=0.6 
□ RHA MVM=0.73 


o 


□ 


n 


emphasize, however, that naturalness is a strong assumption, just as the assumption of 
renormalizability is. In the latter, one sets all but a few of the model parameters to zero, while in 
the former, one maintains that all parameters are roughly of the same magnitude. 
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an infinite nnmber of connterterms must be added to produce finite results. (This will 
happen, for example, if M* = M — $ becomes a more complicated function of <h.) An 
explicit calculation of the counterterms is unnecessary, however, since the end result is 
simply an infinite polynomial in the scalar field, with finite, unknown, and presumably 
natural coefficients arising from the underlying dynamics of the QCD vacuum. To have any 
predictive power, one must rely either on the truncation scheme discussed above, so that 
only a small, finite number of unknown coefficients are relevant, or on some other dynamics 
to constrain the form of the renormalized scalar potential.^ There is, of course, no guarantee 
that an effective hadronic theory should be predictive, and until one can derive the hadronic 
theory from the underlying QCD, the range of predictive power must be determined by 
comparison with experiment. Moreover, since the non-Goldstone bosons are always off mass 
shell in nuclear structure calculations, it is likely that the ability to make realistic predictions 
depends on the choice of generalized coordinates (fields). 

In summary, even in an effective hadronic field theory, one must include loop contribu¬ 
tions that contain negative-energy baryon wave functions, since it is essential to maintain 
the completeness of the Dirac basis, which plays a crucial role in the field theory ||We95|| . The 
new ingredient, compared to renormalizable models, is that there is now an infinite number 
of counterterms. Since the divergent contributions must also respect the symmetries in the 
lagrangian, the appropriate counterterms to remove the divergences will always exist, but 
final, finite results will generally depend on an inhnite number of unknown parameters. Thus 
one must rely on truncation or on additional dynamical input to limit the parameters to a 
manageable number. We will discuss these observations further in Section p. 


6. RELATIVISTIC MEAN-FIELD THEORY 


The background for relativistic mean-field theory (RMFT) and its application to fi¬ 
nite nuclei^ is developed in detail in ||Se86|,^e9^ , |Wa95|| ; we have summarized this ma¬ 
terial in Section Important background references also include [Bo77,IIo81,Re86, 
Fu87,Ru88,Fu89c,Re89,Ga90,Fu91a]. The main conclusion from this material is that by 
using local relativistic Hartree equations for the baryon and meson fields with 

linear couplings and a minimal set of parameters fitted to the properties of nuclear matter, 
one derives the nuclear shell model. 

Moreover, by including phenomenological nonlinear meson couplings of the form 
(k/S!)^^ -f (A/4!)(;/)^, one can extend this model, fit the nuclear compression modulus and 
baryon effective mass ||Bo91| , Fh9^ , and obtain an accurate description of nuclear deforma¬ 
tions in light nuclei ||Bo77 , Bb54 , Re86| , |Fu87|JRu5^ , |Fi89|JSe92| , |Va92| . 


^^For example, a simple model is used in Pu95( | to show how the broken scale invariance of QCD 
leads to dynamical constraints on the scalar potential, and fits to the properties of finite nuclei 
also generate coefficients that are natural. 

^^We sometimes use the terminology “Dirac-Hartree theory” for the RMFT of finite nuclei. 
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A. Density Functional Theory 


To understand the success of the RMFT and to put the calculations discussed in this 
section in context, it us useful to consider some elements of density functional theory 
(DFT). A discussion of DFT applied to nonrelativistic systems can be found, for exam¬ 
ple, in ||Dr90|| . Speicher, Dreizler, and Engel |^p92|] extend the framework to study the 


relativistic many-body problem. These authors outline the density functional approach to 
the strong-interaction model of QHD and apply the Hohenberg-Kohn theorem, widely used 
in ab initio calculations of the structure of solids, to this situation. The nonrelativistic 
version of this theorem can be stated as follows: The ground-state expectation value of 
any observable is a unique functional of the exact ground-state density; moreover, if the 
expectation value of the hamiltonian is considered as a functional of the density, the exact 
ground-state density can be determined by minimizing the energy functional. 

Speicher, Dreizler, and Engel derive an approximate energy functional by using a gra¬ 
dient expansion of the noninteracting kinetic energy to order The energy functional 
includes the effects of four-vector meson exchange and of vacuum contributions, and the 
variational equations of the corresponding extended Thomas-Fermi model are discussed. 
A great deal of work has been done on the extended Thomas-Fermi model, where it 
is relatively straightforward to generalize the results of inhnite nuclear matter calcula¬ 
tions to hnite systems. The extended Thomas-Fermi approach is discussed in detail in 


In more general terms, the central object in a DFT formulation of the relativistic nuclear 
many-body problem is an energy functional of scalar and vector densities (or more precisely, 
vector four-currents). Minimization of the functional gives rise to variational equations 
that determine the ground-state densities. By introducing a complete set of Dirac wave 
functions, one can recast these variational equations as Dirac equations for occupied orbitals; 
the single-particle hamiltonian contains local scalar and vector potentials, not only in the 
Hartree approximation, but in the general case as well.Q Rather than work solely with the 
Dirac wave fnnctions and the resulting densities (as in |[Ni92|,|5c^ , [Fr96|| ), one can introduce 


auxiliary helds corresponding to the local potentials, so that the energy fnnctional depends 
also on classical meson helds. The resulting DFT formulation produces held equations that 
resemble those in a Dirac-Hartree calculation, but correlation ehects can be included, if the 
proper functional can be found. 


The procedure described above is analogous to the well-known Kohn-Sham |[Ko65|| ap¬ 
proach in DFT, which is based on the following theorem [pr90|| (generalized here to rela¬ 
tivistic systems): 

The exact ground-state scalar and vector densities, energy, and chemical poten¬ 
tial for the fully interacting many-fermion system can be reproduced by a col¬ 
lection of (quasi)fermions moving in appropriately dehned local, classical helds. 


^^Note that the Dirac eigenvalues do not correspond precisely to physical energy levels in the 
general case [pr90| , |Sp^ . 
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In the QHD case, the local scalar and vector helds play the role of (relativistic) Kohn-Sham 
potentials, and by introducing nonlinear couplings between these helds, one can implicitly 
include additional density dependence in the single-particle potentials, as well as the com¬ 
posite nature of the nucleon ||Fu96a|| . Thus, even though the Dirac nucleons in an RMFT 
calculation move in local, classical potentials, this does not preclude an exact description of 
the observables mentioned in the theorem. 

The exact energy functional has kinetic-energy and Hartree parts (which are combined in 
the relativistic formulation) plus an “exchange-correlation” functional, which is a nonlocal, 
nonanalytic functional of the densities that contains all the other many-body and relativistic 
effects [|Mu96|| . Rather than try to construct the latter functional from the lagrangian using 
explicit many-body techniques ||Se86|,^e9'^ , ^c95|| , the basic idea behind the RMFT approach 
is to approximate the functional using an expansion in classical meson helds and their deriva¬ 
tives, based on the observation that the ratios of these quantities to the nucleon mass are 
small, at least up to moderate density.p^ The parameters introduced in the expansion can 
be htted to experiment, and if we have a systematic way to truncate the expansion, the 
framework is predictive. 

Thus a conventional RMFT energy functional htted directly to nuclear properties, if 
allowed to be sufficiently general, will automatically incorporate ehects beyond the Hartree 
approximation, such as those due to short-range correlations. Future work can then be 
focused on the explicit inclusion of higher-order many-body ehects (as discussed later in this 
work), to examine the accuracy and limitations of the relativistic Kohn-Sham approach. 

Why should we expect an approximate, mean-held functional to work well? We observe 
that while the mean scalar and vector potentials $ and W are small compared to the 
nucleon mass, they are large on nuclear energy scales |[Bo91 , Fu95 |. Moreover, as is illustrated 
in Dirac-Brueckner-Hartree-Fock (DBHF) calculations |pio87| , [Te57| , [Ma89[ , the scalar and 
vector potentials (or self-energies) are nearly state independent and are nearly equal to 
those obtained in the Hartree approximation. Thus the Hartree contributions to the energy 
functional should dominate, and an expansion of the exchange-correlation functional in 
terms of mean helds should be reasonable. This “Hartree dominance” also implies that it 
should be a good approximation to associate the single-particle Dirac eigenvalues with the 
empirical nuclear energy levels, at least for states near the Fermi surface ||Dr90| . 


We also observe that the nuclear properties of interest include: 1) nuclear shape prop¬ 
erties, such as charge radii and charge densities, 2) nuclear binding-energy systematics, and 
3) single-particle properties such as level spacings and orderings, which rehect spin-orbit 
splittings and shell structure. Since the Kohn-Sham approach is formulated to reproduce 
exactly the ground-state energy and density, and the Hartree contributions are expected 
to dominate the Dirac single-particle potentials, these observables are indeed the ones for 


Since the meson fields are roughly proportional to the nuclear density, and since the spatial 
variations in nuclei are determined by the momentum distributions of the valence-nucleon wave 
functions, this organizational scheme is essentially an expansion in kp/M, for kp corresponding to 
ordinary nuclear densities. Here the nucleon mass M is the generic large mass scale characterizing 
physics beyond the Goldstone bosons. 
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which meaningful comparisons with experiment should be possible. Nevertheless, because 
the Dirac eigenvalues do not correspond precisely to observed single-particle energies (except 
exactly at the Fermi surface), we should not expect to reproduce spin-orbit splittings at the 
same level of accuracy as rms charge radii and total binding energies. s 

As discussed in Section E3 an RMFT energy functional of the form in Eq. ( p.33D , 
extended to include meson self-interactions as in Eq. (|5.49|) , successfully reproduces these 
nuclear observables with parameters of natural size [|Fu96a|| . This justihes a truncation of 
the energy functional at the hrst few powers of the helds and their derivatives, as is evident 
from Fig. Moreover, the full complement of parameters is underdetermined, so keeping 
only a subset does not preclude the possibility of a realistic £t to nuclei. Both the early 
RMFT calculations mentioned above and the newer calculations discussed below should be 
interpreted within the context of this Kohn-Sham approach to DFT, since they typically 
involve truncation at low powers of the helds and include only a subset of the possible 
parameters. We also emphasize that the Dirac-Hartree approach to nuclei is not really a 
Hartree approximation in a “strict” sense, in which one would determine the parameters in 
the lagrangian from other sources (NN scattering, for example) and then solve the mean-held 
equations for nuclei with the same parameters. The DFT interpretation implies that the 
model parameters htted to nuclei implicitly contain ehects of both short-distance physics 
and many-body corrections. 


B. Nuclear Structure 


Recent applications of these concepts exist for a wide variety of nuclei. Properties of light 
nuclei with neutron halos are examined in ||Ta92|JZh94|| . The single-particle structure of odd- 
A nuclei is discussed in |[Fu89c|,Pe94|,|Wa94|,|Ne95|| . “Islands of inversion” in neutron-rich Ne, 


Na, and Mg nuclei are studied in Pa91|. Exotic nuclei near Z = 34 and the proton drip line. 


which play a role in the nuclear r-process and in astrophysics, are studied in [^h93b|,|Ga94a 


The shapes of nuclei with N = Z and 20 < A < 48 are discussed in [[Pa93b|JMa96aj|, those of 


superdeformed Hg isotopes, in [Pa94aj| , and super deformat ion for 140 < A < 150, in ||Af96 
The shapes of neutron-dehcient Pt isotopes are discussed in |^h9 


similar discussions exist 


for Pt, Hg, and Pb in |[Yo94|| , for Sr and Zr in ||Ma92a|| , for Sr and Zn isotopes near the proton 
drip li] 

[La96| 


Hexadecapole moments of Yb isotopes are examined in 


in 

P^ 

. Exotic Ba isotopes 


,La96a 

. Fission barriers 


in heavy nuclei are discussed in |[Ru95| 


Light nuclei are also examined in |[Pa93a|| , and the ehects of pairing are included in 

and 


Pa93c|| . Kinks in the isotope shifts of charge radii near Z = 40 are examined in ||La95 


for the Pb isotopes, in ||Sh93a||. They appear to come out quite naturally in the RMFT. 


The general conclusion of this body of work is that the relativistic mean-held theory 
provides an economical means of describing much of the structure of observed nuclei and a 


^®These arguments are relevant to the choice of weights for the fitting procedure described in 


Section 5.D 
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relatively reliable way to extrapolate to new regions of nuclear structure. 

The charged mesons (tt, p) hrst enter these calculations at the relativistic Hartree- 
Fock (HF) level, where the exchange interaction is included. Systematic studies at this 
level are contained in |[Vo92aj , ^h93| , |Zh9^ , ^u93| , po94| , [Ha9^JSc95aj| . In nonrelativistic Hartree- 
Fock calculations with Skyrme interactions, one makes a phenomenological £t to nuclei 
using a contact NN interaction containing various powers of the density. (A computer 
code to carry out such calculations now exists in the literature |[Re91|] .) Comparisons be¬ 
tween the RMFT and Skyrme calculations of various nuclear properties are presented in 

The effects of retardation and of medium modihcations to the Dirac wave functions are 

The relationship to Landau Fermi-liquid 


examined in more detail in Mi92.Zh91,Za92,Zh92 


theory is studied in ||Ue92| , fra93[| . Extensions to include other types of phenomenological 
nonlinear self-couplings of the meson and baryon helds are discussed in ||Gr92| , F£9^ . 


As noted above, an important goal is to relate calculations of nuclear properties more 
directly to parameters determined from NN observables by explicitly computing many-body 
contributions. Because the relevant equations must be solved self-consistently, this is a 
very difhcult problem for hnite systems, and an important advance was made by Gmuca 
Gm92|,[Gm92a|| . Here the DBHF calculations of nuclear matter |[Br904Pe91|] are parame¬ 


trized by htting the RMFT with scalar and vector nonlinear self-interactions to the DBHF 
results for the energy/nucleon and the self-energies $ and W over a range of densities. (The 
nonlinear parameters used correspond to K 3 , ^ 4 , and Co in Eq. (^.49|) . The important advance 
is that it is more efficient to fit the nucleon self-energies than the effective NN interaction 
(or G matrix), as in |[Br92|,[Bo9^,|Ha94| .) The effective interactions thus obtained are used in 


RMFT studies of the structure of and '^^Ca nuclei without the introduction of additional 
free parameters. The calculated binding energies, single-particle spectra, and charge radii 
agree reasonably well (although not completely satisfactorily) with experimental data and 
present an improvement over the nonrelativistic Brueckner-Hartree-Fock approximation. 
This approach provides a framework for relating the RMFT to the DBHF calculations of 
nuclear matter in a quantitative manner, and thus, ultimately, to the free NN interaction. 
The relationship between the RMFT calculations and DBHF studies of nuclear matter is 
also explored in ||Zh92|,[Za^,|Su94| . 


Lenske, Fuchs, and Wolter have also made an important contribution ||Le95| , [Fu95b|| . Here 
a fully covariant, density-dependent hadronic field theory is obtained in which nonlinear ef¬ 
fects are described through a functional dependence of meson-nucleon vertices on the baryon 
field operators. Rearrangement self-energies arise in the baryon field equations from the vari¬ 
ational derivatives of the vertices. Solutions are studied in the Hartree limit and compared 
to the local-density approximation to DBHF theory. Parametrizations of nonlinear effects 
in terms of the scalar density or baryon (vector) density are discussed. Hartree calculations 
for nuclei between and show that rearrangement corrections simultaneously im¬ 

prove the description of binding energies, root-mean-square radii, and density distributions. 
This approach provides an alternative to the method used by Gmuca for connecting nuclear 
observables to the NN interaction. 

In Section |3.A| , we noted that a linear realization of chiral symmetry, with the usual form 
of symmetry breaking, cannot produce successful nuclear phenomenology at the RMFT level. 
This was the primary motivation for constructing a model with a nonlinear realization of 
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the symmetry in Section ^.D] , which includes a light scalar to simulate the exchange of two 
correlated pions between nucleons. The failure of the linear models, however, lies more with 
the form of the scalar potential responsible for spontaneous symmetry breaking than with 
the linear realization of the symmetry. 

In |[He94a) , |Ca96|| , nuclear matter and hnite nuclei are studied in the RMFT with a chiral 
lagrangian that generalizes the linear a model and also accounts for the QCD trace anomaly. 
A logarithmic meson potential that involves the a and tt helds and also a heavy glueball 
held (p is used to spontaneously break both scale invariance and chiral symmetry. The 
scale-invariant term that leads to an uj meson mass after spontaneous symmetry breaking 
is strongly favored to be of the form by the bulk properties of nuclei; they also 

rather strongly constrain the other parameters. A reasonable description of the closed- 
shell nuclei oxygen, calcium, and lead can be achieved, and the results are improved by 
including a quartic omega self-interaction term in the lagrangian. These results are consistent 
with the discussion at the beginning of this section, since this linear chiral model contains 
meson self-interactions corresponding to the K 3 , /t 4 , and Co terms in Eq. ( b.49|) . (See also 
parameter set Q2 in Table ^.) The important point is that the use of a more general 
(he., nonrenormalizable) interaction potential provides more freedom to adjust the self¬ 
interactions, and the number of free parameters is sufficient to describe nuclei, in contrast 
to the usual (renormalizable) sigma model, as discussed after Eq. (|3.9|) . 

In ||Pr94|| , a single, adjustable scalar-nonlinearity parameter is included in the energy, 
together with the zero-point energy AT discussed in Section |5.E| . Although observed neu¬ 
tron star masses do not constrain this parameter, it can be chosen to provide a reasonable 
description of bulk nuclear properties, when the scalar mass is roughly 600 MeV. In this case, 
however, the coefficients of the nonlinear and terms are unnatural; in particular, the 
quartic coefficient is unnaturally large and negative. This is necessary to cancel the effects 
of the unnaturally large, positive coefficients introduced by AT [see Eqs. (|5.58|) and (|5.61|) ]. 


C. Electroweak Interactions in Nuclei 


1. Electromagnetic currents in QHD 

In QHD-I there are no charged mesons. One can, however, introduce an effective elec¬ 
tromagnetic current, to be used in lowest order, that incorporates some of the nucleon’s 


internal structure 3e86.Se92,Wa95 


J), = ^ , 

Q = 2(1 + '^3) , A' = Ap -(1 Ta) A[j -(1 - Ts) . 


( 6 . 1 ) 


Here Ap = 1.7928 and A(, = —1.9131. This current is covariant and local, and it is conserved 
by virtue of the QHD-I field equations. It also contains the correct anomalous magnetic 
moment. To include the spatial extent of the nucleon, one can introduce a single overall 
form factor 

= [ . . ( 6 . 2 ) 


l-g7(855MeV)2 
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where is the four-momentum transfer, to be used in all matrix elements of the current 
(3.1). This is equivalent to replacing the photon propagator 1/g^ with the effective Mpller po¬ 
tential This current can be used consistently with the RMFT solutions in QHD- 

I, such as the Dirac-Hartree wave functions. Many applications exist ||Se86| , ^e9^JWa95| . 

The use of a single function /s.n.(Q'^) assumes that the momentum dependence of the 
charge and anomalous form factors is the same. This assumption breaks down at large 
since it is actually the Sachs form factors that scale similarly. This observation can be 
implemented easily with a simple change to /s.n.(?^) and ||Wa95| . 


Alternatively, one can attempt to calculate the single-nucleon structure by starting with 
the electromagnetic current in the renormalizable model QHD-II, which contains charged 
p and TT mesons, and then by evaluating quantum loop diagrams. This is discussed in 
|5i86i 


The structure of the electromagnetic current in QHD-III 


I, a chirally invariant. 


renormalizable extension of QHD-II, is currently under investigation ||Pr96| . 

In the context of effective held theory, the composite structure of the particles is de¬ 
scribed with increasing detail by including more and more nonrenormalizable interactions in 
a derivative expansion. For example, the electromagnetic current obtained from Eq. (|5.48|) 
by taking 6C/6{eA^) is, after some partial integration. 


j*'= Iw + r 3 h'-iy + iVAV-iV) - ^ 


a2[Ar(/?(o) +^(1)7-3)^/^ iV] 


+ -(d,pr + + 2/pr(u'‘Ts) . 

9'y 


(6.3) 


Note that the photon can couple to the nucleon either directly or through the exchange of 
neutral vector mesons (rho or omega). 

We can determine the tree-level electromagnetic form factors of the nucleon from the 
current (^) and the lagrangian (5.4(j|). For spacelike momentum transfers = —q^, the 


isoscalar and isovector charge form factors are 

_ 1 /3'°’ Q"- 




Sv 








2 M2 

3g-y + ml 


9p 

2 M2 

2g^ Q2 + ^2 

are 


Xp + X'n 

fv9v 

2 

3g^ Q2 + ml 

a;-a; 

fp9p ^ 


+ ••• 




+ ■ ■ • , 


(6.4) 

(6.5) 

( 6 . 6 ) 

(6.7) 


The corresponding mean-square charge radii are 

, 2g, 


(r=)f = 


3g^ml 
4 /v^v 
K + K 9'i'ml ’ 






6 


fp9p 
K - K 9^rnl 


( 6 . 8 ) 

(6.9) 
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The form factors have a contribution from vector dominance and a correction from the 
intrinsic structure of order Q^, that is, to second order in a derivative expansion. This 
correction is adequate for most applications to nuclear structure; thus, when the model 
of Section b-D| is applied to finite nuclei, an ad hoc form factor need not be introduced.^ 
In practice, as the values of and Qp are varied to £t nuclear properties, I3^^\ and 
fp are chosen to reproduce the empirical isoscalar and isovector charge radii and isovector 
anomalous radius of the nucleon. (The parameter is determined from nuclear properties, 
since the nucleon’s isoscalar anomalous radius is poorly known.) 

We note, however, that although the derivative expansion for the current is adequate 
for most RMFT calculations, applications involving large energy-momentum transfers will 
require additional terms of higher order in Q^. The utility of this expansion for ^ is 
an open question that is currently under active investigation. 


2. Weak currents in QHD 

For the effective weak currents, one can proceed analogously. The effective weak current 
to be used in lowest order in QHD-I is given as the sum of a polar-vector and an axial-vector 
part by |^e86| , |Wa95[| 


_l_ 7± 

fj, fi ‘ fib 1 

where the charge-changing, polar-vector current is dehned by 

(a; - A') 


( 6 . 10 ) 


= '07/.T-±V’ + 


2 M 




( 6 . 11 ) 


with = |(ti ± iT 2 )- Similarly, the charge-changing, axial-vector current is given by 


= Fa{ 0 ) g. 


rrii -I- 92 


dpd" . 


( 6 . 12 ) 


The weak current so dehned is covariant, it satishes PCAC (partial conservation of the 
axial-vector current), it contains a nonlocality implied by pion-pole dominance, and it gives 
the correct result for semileptonic weak interactions on a free nucleon. 

The effective weak neutral current for QHD-I is dehned to have the symmetry properties 
of the standard model ||Se86| , |Wa95| 


yu _ tu I tu 
fj, ~ ^ "'/x5 


2 sin^ Oy^Jp 


(6.13) 


where and are obtained from Eqs. (|6.11| ) and (|6.12| ) by the replacement t± —> ^Ts, 
and Jp in the hnal term is the electromagnetic current. Several applications of these currents 
The weak current in the ehective theory of Section |5.D| remains to be 


exist [3e86,Wa95 


studied. 


form factor for the pion also arises directly from vector-meson dominance [Fu96a|. 
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3. Electroweak exchange currents 


Within the framework of a consistent hadronic held theory, one can also calcnlate two- 
body exchange currents. Early work on electromagnetic exchange currents is summarized 
and referenced in |^e 86 || , where explicit expressions are given for the long-range pion and 
pair currents. By the conserved-vector-current theory (CVC), these also give the exchange- 
current contributions to the weak vector current. Exchange currents play a crucial role 
in the accurate reproduction of the high-momentum-transfer behavior of electromagnetic 
observables for light nuclei, as discussed in Section |4.A . 

One can also compute the exchange-current contributions to the weak axial-vector cur¬ 
rent. Instead of strict current conservation, one must now respect chiral symmetry and 
PCAC ||Wa95|| . While the forms of the relativistic one-body and two-body axial-vector cur¬ 
rents are model dependent, a simple and useful version arises from the chirally transformed 
(nonlinear) sigma model discussed in Section |3.B|, where a projection operator of the form 


in Eq. ( |6.12|) appears directly in the current ||An96|| . The long-range pion-exchange-current 


contribution to the isovector axial charge operator in this approach is given by 


J05(xi,X2,x) = 


dvr 


[r(l) X t(2)] [5(xi - x)cr 2 (5(x2 - x)cri] ■ f (1 -F x^) 


exp (-x^) 


Xt 


(6.14) 


where = m 7 r|r|, with r = Xi—X2, and f = g.,^/2M. The Goldberger-Treiman relation, 
which follows from PCAC, states that MF^iS^) = —g^fn, where ~ 93MeV is the pion 
decay constant ||Wa95|| . Many applications of the above expression, including additional. 


shorter-range, hadronic-exchange contributions, also appear in the literature ||To95|| .P^ 


4. Recent developments in electromagnetic interactions 

The simplest and most informative electromagnetic process to study is electron scattering 
(e, e'). Relativistic analyses of scattering in the region of the quasielastic peak, corresponding 
to single-nucleon knockout, are found in |Pe92| , PT92|JEi94|jni9^JEr94| , pi95|| . The systematics of 
the location and shape of the peak, from low energies to the high-energy data from SLAC on 
®®Fe, is examined in |pe92| , |Fr9^ . In the former work, the roles of M* and exchange currents 
at high are emphasized. In the latter, it is shown that one must include momentum 
dependence in the relativistic self-energy S(fc) to understand the data; the RMFT gives a 
constant S. 

The failure of the (e, e') data to satisfy the Coulomb sum rule throughout the periodic 


table has long been one of the signihcant puzzles in nuclear physics ||Wa95|| . The correct 
form of the relativistic sum rule, the effects of nuclear binding, and the use of the off-shell 
current are analyzed in ||Ee94|,[KoU^. An important recent contribution ||Jo95|| indicates that 


the solution to this problem may now be in hand. In this work, the world data on inclusive 


^^Equation (|6.14 ), based on pion-pole dominance of the relevant amplitudes, contains an extra 

wm - 


factor of E 4 relative to that used in 


69 








































1 


1 




1 


1 


1 



p„(MeV / c) 

FIG. 10. The (e, e^p[3si/2]) reduced cross section of for high missing momentum with 

incoming kinetic energy Ei = 487 MeV and the binding energy Eh = 10 MeV. The solid and 
dotted lines are the results of Kim and Wright ||Ki96 |, who use an RMFT model with a spectroscopic 
factor fitted to low-missing-momentum results. The solid line includes Coulomb distortion in an 
approximate way, while the dotted line is the plane-wave result. The dashed line is a nonrelativistic 


result by the Belgium group [ Va96 |, which includes only the one-body current. The data are from 
NIKHEF Po94b|| . 


quasielastic electron scattering have been used to separate the longitudinal and transverse 
response functions of ®®Fe and The resulting longitudinal response functions lack the 
“quenching” that has been such a problem. There are still some inconsistencies in the total 
collection of world data, however, and this is a subject that cries for a complete experiment 
with all the kinematic flexibility of CEBAF for dehnitive resolution. 

With coincidence capabilities at extreme kinematics now available at Bates, Mainz, 
NIKHEF, and CEBAF, one can use the reaction (e, e'X) to probe, among other things, 
the high-momentum tail of the momentum distribution in the nucleus. There are many 
effects that will contribute at high momentum transfer and at high missing momentum, in 
(e, e'p) from a hlled orbital. For example, in ^°^Pb(e, e'p) from the 3si/2 level, one must 
consider long- and short-range correlations, meson-exchange currents, isobar currents, the 
spin-orbit interaction, other relativistic effects, etc. One of the prime motivations for such 
experiments is to disentangle and study these effects. The RMFT analysis of such a pro¬ 
cess, where the initial wave function is the solution to the Dirac-Hartree equations and the 
hnal wave function is generated in the RIA optical potential, has the distinct advantage 
that much of this physics, as in nucleon-nucleus scattering (Section p.Q) , is incorporated 
in the starting approximation. Figure 0, for example, shows recent calculations of Kim 
and Wright ||Ki96|| for this transition. This calculation, with no free parameters, accurately 
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reproduces the existing data. It is only by pushing to more extreme kinematics, with preci¬ 
sion experiments, and with many nuclear transitions, that one will be able to unambiguously 
disentangle the role of the various additional contributions. 

The subject of (e, e'X) is also discussed within the framework of QHD in [Pi92| , |Ji93i|Ga94| . 


It is emphasized in the hnal paper that one must understand relativistic hadron dynamics 
before drawing conclusions about new phenomena, such as color transparency. 

Relativistic analyses of the photonuclear reactions ( 7 ,p) and ( 7 , 7 r“p) are discussed in 
[Lo92|jHe^JJo94|| . Other phenomena such as {auj'y) mixing, {poj) mixing, charge-symmetry 
breaking, and the production of (e’''e“) pairs to probe the quark-gluon plasma are discussed 
in |K?9^ , |Cl95a| , [Cl9lq . 

In a recent review article on Electromagnetic Response Functions in Quantum Hadrody- 
namics ||We93|| , elastic, inelastic, and quasielastic electron scattering over a wide range of 
energy and momentum transfers are discussed as a probe of the electromagnetic response of 
nuclei. Electromagnetic response functions obtained with different QHD models of nuclear 
structure and at different levels of approximation are compared with data. It is shown that 
RPA correlations are important and have different effects in different kinematical regimes. 

A second recent review of Nuclear Response in Electromagnetic Interactions with Com¬ 
plex Nuclei is given in ||Bo93|| . Here the response of nucleons and complex nuclei to an ex¬ 
ternal electromagnetic probe at intermediate and high energies is illustrated by considering 
both inclusive and semi-inclusive electron scattering. Form factors for elastic scattering and 
structure functions for inelastic scattering are derived, and several examples are discussed, 
including polarization observables. 


5. Recent developments in electroweak interactions 

One can obtain information about the single-nucleon matrix elements of the weak current 
by studying quasielastic neutrino scattering (z/, u') from nuclei. The importance of the 
single-nucleon self-energies and of a realistic (RPA) description of the nuclear excitation 
spectrum in the extraction of the axial-vector form factor of the nucleon F^iq^) is discussed 
in |pi)n93| , |KiU^ . 

The single-nucleon matrix element of the weak axial-vector current contains a one-pion- 
exchange contribution leading to the induced pseudoscalar coupling Fp(g^), and the weak 
axial-vector exchange current also has a one-pion-exchange contribution (|6.14| ); one might 
expect this longest-range contribution to be modihed in the nuclear medium. A relativistic 
hadronic held theory of the nucleus (QHD) allows one to estimate these corrections. A great 
deal of recent work has progressed in this direction ||To92|,|Ga92|,Pa93|,[tz94|JPaM|,|Gi95|| , much 
of it carried out within the RMFT of nuclear structure. The basic conclusions are that weak 
axial-vector exchange currents can play an important role if one looks in the right place, 
for example, in hrst-forbidden /d-decay, and that there can be signihcant nuclear-structure 
effects in the weak axial-vector coupling to the nucleus. 

In |P?o92||, the “enhancement factor”, dehned as the ratio of the axial-charge matrix 


element in hrst-forbidden /3-decay to its one-body value, is calculated in a meson-exchange 
model for both light and heavy nuclei. Pion-exchange processes are computed for a chirally 
symmetric, phenomenological lagrangian and compared with results obtained in the soft-pion 
approximation [Eq. ( |6.14 )]. Heavy-meson pair graphs are included, with coupling constants 
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determined from the Bonn NN potential |[Ma89|| . Nonrelativistic reductions are performed 
to both leading order and next-to-leading order. The results are sensitive to the choice of the 
short-range correlation function used in conjunction with harmonic oscillator wave functions. 
A ratio of transition matrix elements, however, is less sensitive to this choice. The main 
result in this paper is the comparison of a transition in a heavy nucleus with a transition 
in a light nucleus: r = (5mec(^ = 208; l 5 f 9/2 ^ 0 / 19 / 2 )/(5mec(^ = 16; lsi /2 ^ Opi/ 2 ) = 1.38. 
This result is smaller than the value of 1.58 ± 0.09 deduced by Warburton ||Wa94a|| in a £t 
between one-body shell-model matrix elements and experiment. 


In |pr’a94||, shell-model matrix elements of the axial-charge exchange-current operator are 


calculated through next-to-leading order in heavy-fermion chiral perturbation theory. It is 
found that the loop corrections to the one-soft-pion-exchange contribution in Eq. ( |6.14|) are 
small (roughly 10%) and have no signihcant dependence on the nuclear mass number or on 
the valence-nucleon orbits. 


Radiative /i-capture is studied in the RMFT in ||he92 


6. Parity violation in (e, e') 

One of the most important areas of future research at CEBAF will be the study of parity 
violation in electron scattering. This arises from the interference between the amplitudes for 
photon exchange and for exchange, which involve couplings to the electromagnetic and 
weak neutral currents, respectively [|Wa95|| . Feasibility has been demonstrated in pioneering 
experiments at Bates and Mainz. 

Two particularly interesting contributions are [ A193 , paM| . In [|A193|| , the impact of pionic 
correlations and meson-exchange currents in determining the (vector) response functions for 
electroweak, quasielastic lepton scattering from nuclei is discussed. The Fermi-gas model 
is used to maintain consistency in treating forces and currents (gauge invariance) and to 
provide a Lorentz-covariant framework. Results obtained in first-order perturbation theory 
are compared with HF and RPA calculations and are found to provide quite successful 
approximations for the quasielastic response functions. The role of pionic correlations is 
investigated in some detail, and meson-exchange currents are shown to provide a small, but 
non-negligible contribution to the vector response. 

In the second paper, parity-violating quasielastic electron scattering is studied within the 
context of the relativistic Fermi-gas model and its extension to include pionic correlations 
and meson-exchange currents. The work builds on previous studies using the same model; 
here, the part of the parity-violating asymmetry that contains axial-vector hadronic currents 
is developed in detail, and a link is provided to the transverse vector-isovector response. 
Various integrated observables are constructed from the differential asymmetry, and the 
most favorable observables for studying pionic correlations and the strangeness form factors 
of the nucleon are determined. Comparisons are also made with recent predictions based on 
the RMFT. 

Nuclear-structure effects pertinent to the extraction of the weak, neutral, axial-vector 
form factor of the nucleon, a quantity of particular interest because of the potential role 
played by ss quark-antiquark pairs, are emphasized in ||Ho93| , po93a . 


Parity-violation in the structure of heavy nuclei is examined in the relativistic HF ap¬ 


proximation in |Ho94|. The role of nuclear structure in atomic parity-violating experiments 
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is discussed in the definitive work Po92a 


D. Strangeness in Nuclei 


The addition of strangeness adds another dimension to nuclear structure. Properties of 
hypernuclei with S = —l are investigated in the RMFT in ||(Jo93|JGTUB| , po95| , [Lo95|| . The two 
primary systematic features of A hypernuclei are the relatively small depth of the central 
potential and the small spin-orbit splittings. The scalar and vector couplings in a (a, a;) 
model can be adjusted phenomenologically to describe the central potential, but the spin- 
orbit splittings have been a more vexing problem. 

It was observed in |Te^ , PoMl that a tensor coupling of the oj to the hyperon, with 
the appropriate sign, can indeed produce a small spin-orbit splitting. The systemat- 
ics of the interaction of hyperons with this additional tensor coupling are investigated in 
|Co92| , Po94| , |Co95| , [Ma94| , |Ma95|| . In particular, it is shown in ||Co94| , |Co95|| that a relativistic 


optical potential with a tensor coupling can describe the data. 

If one views QHD as an effective hadronic field theory, such a coupling will exist in 
the lagrangian. Even within a renormalizable framework, one has an induced coupling of 
this form. While motivation for a tensor coupling appears to exist within a quark-model 
framework, the theoretical challenge is to quantitatively explain the difference in size between 
spin-orbit effects in the hyperon and nucleon sectors. 

One interesting goal of the study of strangeness in nuclei is to extrapolate to large [S'!; the 
RMFT provides a convenient basis for this extrapolation, and this problem is examined in 
|Sc92|,Po9^[Ma93^,|Sc93|,|Ga95|| . A large class of bound, multi-strange objects is one distinct 


possibility raised by these studies. 

The has a relatively weak coupling to nucleons, and thus it provides a useful hadronic 
probe for viewing the nuclear interior. scattering and production from nuclei are studied 
within the framework of the RMFT in ||Ca92| , pa92|,[Ru92|,|Ko95a)| . In [[Ko95b|| , it is shown that 
K~ quasielastic scattering from nuclei can be explained within the framework of the local- 
density RPA and isoscalar correlations contained in QHD-I. Subthreshold kaon production is 
discussed in |[Li95|| , and strong-interaction effects in atoms, in ||Ma95a|| . Other properties 
of hypernuclei are examined in ||Lo94| . 


On a more basic level, the properties of a hyperon in nuclear matter are examined within 
the framework of QCD sum rules in 


7. MATTER UNDER EXTREME CONDITIONS 
A. Relativistic Transport Theory 

One of the principal thrusts of nuclear physics has been, and will continue to be, 
the use of relativistic heavy-ion reactions to study the properties of nuclear matter un¬ 
der extreme conditions of density, temperature, and flow. A great deal of work has 
been done within the framework of QHD using relativistic transport theory to describe 
heavy-ion collisions. For example, the foundations of relativistic transport theory are dis¬ 
cussed in ||Ko87| , pi55a] , [Ko88| , p.ji89a| , [Ca9U| , pVfa92|jMa93| , pVfo94|JSc94|| and further developed in 
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Pe91|JDe9^JMe93yNi93| , |Mi^ , [Fu95a| , Po95|| . In this section, we start the discussion with a 
simple introduction to transport theory. 

The Vlasov-Uehling-Uhlenbeck (VUU) model describes the transport of the distribution 
function in phase space. It is a one-body transport model that includes the effects of a long- 
range mean held on the one-body dynamics, together with a short-range, two-body collision 
term. The mean held is calculated from the distribution function, and when parametrized 
as a function of the density, it probes the equation of state of the medium. The dynamical 
description is classical, except that collisions are treated stochastically, and a Pauli-blocking 
factor is included for the hnal states in the collision term. The approach can be generalized 
to include inelastic channels. The basic ideas behind various transport models are described 
in the review article by Greiner and Stocker [pt86|| , and a computer code to carry out such 

|Ha93a|| 


calculations is now available 


We review some of the basic concepts. 


Consider the microcanonical ensemble that consists of a collection of identical, randomly 
prepared, microscopic systems (particles). The ensemble can be characterized by the distri¬ 
bution in phase space: 


dA^ = Number of members of the ensemble in 


d^p d^q 
(27r)3 


= /(p,qG) 


d^p d^g 
(27r)3 


(7.1) 


The probability to hnd a particle in this region of the ensemble is the probability of picking 
such a member of the ensemble at random; this is equal to dN/N, where N is the total 
number of members in the ensemble, obtained by integrating Eq. O) over all phase space. 
This probability can be used to compute expectation values over the ensemble. 

As a function of time, a particle moves from (po, qo) to (p, q) in phase space. Liouville’s 
theorem (see, for example, |[Wa89|| ) states that with hamiltonian dynamics, the volume in 
phase space is unchanged with time. Since the number of particles diV in this volume is 
also conserved, one concludes that the distribution function is unchanged along a phase 
trajectory: 


/[p(i).q(<).*l = /[Poiqo.<ol = constant. 


(7.2) 


Differentiation with respect to time yields 


dt dt dp dt dq dt 


(7.3) 


Now Hamilton’s equations state 


dq dH dp dH 

dt ^ (9p ’ dt ctq 


(7.4) 


and insertion into Eq. Q gives 


dt dq'dp dp'dq ^ 


(7.5) 
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The last equality identifies the Poisson bracket. In equilibrium, there is no time dependence 
to the distribution function, and this expression must vanish. The solution to this condition 
is that / = f{H), and / is a constant of the motion. 

A short-range, two-body Boltzmann collision term is now included. Go to the canonical 
ensemble, which consists of identical, randomly prepared assemblies (collections of systems). 
Particles are now scattered into and out of the region in phase space in Eq. (|7.1|) . Assume 
that 


dt 



collisions 


(7.6) 


Momentum is conserved in the collisions, so that Pi -l- P 2 = p) + P 2 - Detailed balance then 
states that Rate(z —>■/) = Rate(/ —*• i), or v^a = Vii 2 '<j', where Vij is the relative velocity.^ 
The number of transitions per unit time in the direction i —>■ f in the assembly is given by 


number of transitions 
time 




a 


/(P2,qG) 


d^P2 d^q 
(27r)3 


(7.7) 


The first factor is the incident flux, and the last factor is the number of target particles, 
both in the appropriate momentum interval. The Boltzmann collision term at the point q 
can thus be written as [^t86| 



collision 


= J d^P2 dVi dV2 (p + P2 - pi - pi)c^G2 

X [/(pi, q, t)f {P 2 , q, t) - /(p, q, t)/(p 2 , q, t)] • 


(7.8) 


The last term counts particles scattered out of this region in phase space, and the first term 
counts those particles scattered in. This relation will be abbreviated as 

^ (^ ~ fh](^vi2 + P2 - pi - p' 2 ) • (7.9) 

This is the result of classical transport theory; with zero-range collisions, the angular distri¬ 
butions will be isotropic. 

If the one-body hamiltonian has the form 

77 = R + c/(q) . (7,10) 

where 7/ is a mean field, then the transport equations take the form 

(7.11) 


^^Note that the number of final states in momentum space = d^pf/{2 tt)^ has now been 

explicitly removed from the definition of a. 
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These are nonlinear, integro-differential equations for the distribution function /. In equi¬ 
librium, with no time dependence, the collision term must vanish, which implies = // 2 , 
or f{E[)f{E 2 ) = f{E)f{E 2 ). Since energy is conserved in the two-body collision cross 
section, one has E + E 2 = E[ + E'^. If these relations are to hold for all E and 
E 2 , then f{E)f{E 2 ) = g{E -|- E 2 ). Differentiation with respect to E and E 2 in turn 
leads to f'{E)/f{E) = constant = —l/k^T, which yields the Boltzmann distribution 
f{H) = exp {—H/ksT) as the equilibrium solution to these kinetic equations. 

In molecular dynamics calculations [with no f/(q)], one simply follows the classical equa¬ 
tions of motion of all the members of the canonical ensemble numerically. In contrast, 
in the Boltzmann-Uehling-Uhlenbeck approach [still with no f/(q))],Qone includes a Pauli- 
blocking term for identical fermions that prevents them from scattering into occupied states. 
An occupied state is assumed to consist of one particle in the unit cell d^pd^g/(27r)^ in phase 
space. Thus one makes the replacement 

if'ji - ff2] ^ if'jni - /)(i - /2) - //2(i - /:)(! - m . (7.12) 

The equilibrium solutions to these kinetic equations are the familiar Fermi distribution 
functions. 

In the VUU approach, one also includes the additional long-range mean field D(q) 
produced by the local particle density p, which is obtained in turn by integrating the dis¬ 
tribution function over d^p/(27r)^. Parametrization of the p dependence of U then probes 
various equations of state. 

What advantages does QHD have for studying these transport equations? 


It provides a covariant description of the nuclear many-body system. 

It thus provides a basis for a relativistic treatment of the transport equations. 

It allows for a consistent treatment of all hadronic channels. 

In the RMFT, it provides an excellent hrst approximation to the nuclear mean held, in¬ 
cluding an energy dependence in the equivalent nonrelativistic optical potential ||Se86 


A great deal of work on relativistic transport theory, mostly in connection with relativistic 
heavy-ion reactions, has been done in the past few years. A covariant Boltzmann-Uehling- 
Uhlenbeck (BUU) approach, the basic ideas of which have been presented above, is developed 
B188a| , [ck^[Ma92| , [Ma94a)| and applied in ||B189| , |B19'T|jKo9CI||Ko91|| . Relativistic transport 


m 


coefficients are discussed in [|Ha93b| , |Ay94| , |Mo9'4 


The connection of the scalar and vector mean helds to the underlying relativistic two- 
body theory is explored in ||E192|JFu9^ , |We92|| . Effective cross sections in the medium are stud¬ 
ied in ||Li93|,^Ta94b|| . Additional momentum-dependent scalar and vector potentials, which 


provide a more accurate description of the optical potential, are introduced in ||We92| , [We93a 


Shock waves are discussed in ||Mo95|| . The role of the Dirac sea in such collisions is discussed 


40 


This is sometimes called the Nordheim-Uehling-Uhlenbeck approach. 
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in 


Ju92|. The production of kaons is examined in [Fa93,Fa93a,Fa94| and of antinucleons, in 


fMlfI^j Li9^ 


Interesting extensions include studies of a classical version of QHD |[Bu93| , Pu93a| , |Bu95 
and of a transport theory for quarks and mesons ||Zh92a|| . Experimental implications of a 


relativistic, mean-held, two-huid model are explored in ||lv94|,|Ru9^ . 


B. Extrapolation and Connections to QCD 


The original motivation for the development of QHD was to hnd a theoretical framework 
that allows extrapolation of the properties of observed nuclear matter to extreme conditions 
of density and temperature, while retaining general principles of physics, such as quantum 
mechanics, Lorentz invariance, and microscopic causality. The only consistent theoretical 
framework we have for describing such a strongly interacting quantum many-body system 
is relativistic quantum held theory based on a local lagrangian density. The hrst model 
attempt using hadronic degrees of freedom, QHD-I, consisted of a renormalizable theory 
with baryons and neutral scalar and vector mesons. When solved at the level of mean- 
held theory (RMFT), one hnds an equation of state with a minimal number of parameters 
that can be htted to the equilibrium point of nuclear matter, which then provides a simple 
equation of state at all densities pe and temperatures T. These results have been discussed 
in Section |. 

As noted earlier, when viewed as an ehective hadronic held theory describing the under¬ 
lying QCD, all possible couplings consistent with the relevant symmetries should be included 
in the lagrangian. Without some organizing principle, one soon loses predictive power. The 
organizing principle presented here is that while the mean helds are large on the scale of the 
nuclear binding energy, the dimensionless ratios ^/M and W/M are still relatively small, 
and one can sensibly expand in these. This observation allows one to understand the suc¬ 
cess of the RMFT treatment of nuclear structure. While providing a systematic basis for 
extrapolating away from the properties of observed nuclear matter, this approach clearly 
limits the extent over which one can reliably carry out this extrapolation. 

A major goal of the extrapolation to high baryon density pe is to describe the properties 
of neutron stars. Extrapolations within the framework of RMFT with effective couplings are 
studied in [|Ro92| , |GT9^ , |Ba94a| , |Mo944[Fr94| , [Su9^ , |Mu96|| . In ||Mo94a|| , cold nuclear matter is 
investigated in nonlinear, mean-held, scalar-vector models that include density-dependent 
meson parameters. This dependence can be both explicit and implicit through the nu¬ 
cleon effective mass. Interactions between the scalar and vector helds are included, and 
the properties of neutron stars are investigated using the resulting equation of state. The 
implications of the model of Zimanyi and Moszkowski ||Zi90| 
plings are examined in | G192a , Fa94a , Pe95 


with derivative scalar cou- 
The connection between the RMFT results and 

Dirac-Brueckner theory is studied in ||Li92| . 

A study of the uncertainties in extrapolation is contained in |[Mu96|| , where the proper¬ 
ties of high-density nuclear and neutron matter are computed using a relativistic mean- 
held approximation to the energy functional. Various types of nonlinearities involving 
scalar-isoscalar (a), vector-isoscalar (a;), and vector-isovector (p) helds are introduced to 
parametrize the density dependence of the energy functional. The model parameters are 
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calibrated at equilibrium nuclear matter density, and it is possible to build different mod¬ 
els that reproduce exactly the same nuclear matter properties, but which yield maximum 
neutron star masses that differ by as much as one solar mass, even when the nonlinear pa¬ 
rameters are restricted to be of natural size. Moreover, with enough nonlinear couplings, 
one can reduce the predicted maximum neutron star mass to Mmax/d^© ~ 1-6, which is 
only 10% larger than that of the most massive observed neutron stars. Implications for the 
existence of kaon condensates or quark cores in neutron stars are discussed. 


Properties of neutron stars are examined within a similar RMFT framework in ||Ma92c 
and within DBHF theory in |[hlu94 


clear matter into a 
Lo92a|JLo94a| . 


Teierls’ type’" 


The interesting possibility of the transition of nu- 
of periodic structure is investigated in the RMFT in 


Another goal of QHD is to provide a reliable description of the hadronic phase of matter 
as one approaches the transition to the quark-gluon plasma (QGP) ||Se92|,[Wa^ , assuming 


that such a transition exists. The role of the QGP in the neutron matter equation of state 
is examined in |[Ro92|| and more generally in ||Ma93b|,pr95| , [Ka95| , [Ri95| . 


The less spectacular, but more easily accessed liquid-gas phase transition is studied 
within the framework of RMFT in ||So93| , pa94a) , |So94|JMu95|| . In ||Mu95|| , a RMFT model of 


nuclear matter with arbitrary proton fraction is studied at hnite temperature. An analysis is 
performed of the liquid-gas phase transition in a system with two conserved charges (baryon 
and isospin) using the stability conditions on the free energy, the conservation laws, and 
Gibbs’ criteria for phase equilibrium. For a binary system with two phases, the coexistence 
surface (binodal) is two-dimensional, and thus the liquid-gas phase transition is continuous 
(second order by Ehrenfest’s dehnition) rather than discontinuous (hrst order), as in familiar 
one-component systems. Using a mean-held equation of state calibrated to the properties 
of nuclear matter and hnite nuclei, various phase-separation scenarios are considered, and 
the model is applied to the liquid-gas phase transition that may occur in the warm, dilute 
matter produced in energetic heavy-ion collisions. In asymmetric matter, instabilities that 
produce a liquid-gas phase separation arise from huctuations in the proton concentration, 
rather than from huctuations in the baryon density. 

The description of nuclear matter at higher temperature is also of great interest; for 
example, very hot nuclear matter existed in the early universe. A similar substance will 
be produced at the Relativistic Heavy-Ion Gollider (RHIG), and supernovae also involve 
a large range of pe and T. An analysis of hot nuclear matter within the framework of 

mm 


QHD is developed in 


General principles of covariant thermodynamics and thermo¬ 


dynamic consistency are discussed, and these principles are illustrated by computing nuclear 
matter properties in an arbitrary reference frame, using the RMFT of QHD-I. The results 
are shown to be Lorentz covariant, and thermodynamic consistency is demonstrated by 
proving the equality of the “thermodynamic” and “hydrostatic” pressures. The mean-held 
results are used in a simple hydrodynamic picture to discuss the phenomenology of heavy- 
ion collisions and astrophysical systems, with an emphasis on new features that arise in a 
covariant approach. Govariant Feynman rules for going beyond the RMFT are discussed in 
Fu91|JFu9Tb[ . 


The nucleon mean free path is considered in [|Ha92|,|Ha94b 


thermal huctuations and 
quantum corrections are examined in |[Ku91| , ^924IQi93|,|Ag94| , ^u95a|| , and screening is dis¬ 


cussed in [|Ga94c||. Supernovae and neutron stars are studied together in ||Su92|| . Statistical 
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properties are examined in |[Vo92| , pa93|| , and collective modes at finite T in ||Ni93|| . Various 


aspects of relativistic heavy-ion physics are discussed in [Pu95 


A novel 1/N expansion of the theory is pursued in |[la924|He93| , [la93|JHe95|| , where, for 
isospin, one has SU{N) with N = 2. The thermodynamic potential with exchange and RPA 
corrections to the RMFT is calculated within this framework. 

A key element of the RMFT of nuclear matter in QHD-I is that the baryon acquires a 
density-dependent effective mass M* due to the classical scalar field, which in turn arises 
from the baryons themselves. M* must be determined self-consistently at all densities and 
temperatures, and in this model, M*/M goes to zero at high and high T. The successful 
description of the location and shape of the quasielastic peak in (e, e') at the RMFT level 
gives additional evidence for the model values of M* /M in the nuclear medium [^e86| , [Wa95 


There is no shift in the meson masses at this level of approximation, although meson masses 
are modified when one includes polarization insertions |^e86| , |Wa95[| . The modification of 
meson masses is crucial, for example, for deciding whether pion or kaon condensation takes 
place in the nuclear medium, and at what density. The modification of the masses of vector 
mesons can be studied experimentally by looking at in-medium lepton-pair formation; there 
are plans in place to do so at CEBAF. 

A great deal of work in the past few years has gone into the investigation of the mod¬ 
ification of hadron masses in the nuclear medium | Ga94b , We94 , [?a95a , [Ko95cJLi95c , [Sa95 
The effect on the nuclear force is examined in 
Ko95dJLi95c|| . In 


|Ga95a|| and on relativistic transport in 


aspects of chiral symmetry and their implications for changes in 
hadron structure in the nuclear medium are summarized and discussed. This includes issues 
such as the density dependence of the chiral (quark) condensate, the related appearance of 
strong scalar mean fields in nuclei, the stability of the pion mass against compression in 
dense matter, and recent explorations of s-wave kaon-nuclear interactions. Kaon condensa¬ 
tion is also discussed in |[Ma94d| , |Sc94a| , ^I^ . Fortunately, there is an extensive survey article 
available on this subject in ||Ad93||. Here various scenarios for chiral symmetry restoration 


and deconfinement at finite temperature or density are studied assuming universal scaling 
relations for some hadron masses. 

An interesting extension of QHD focuses on the incorporation of the broken scale invari¬ 
ance of QCD. Reference ||Ja93a|] examines various chiral lagrangians in which the QCD scale 
anomaly is implemented by introducing a dilaton field representing the gluon condensate. 
The lagrangians are used to study the chiral phase transition and the propagators of both 
scalar and vector mesons in dense baryonic matter. A hybrid model is proposed that allows 
for an unambiguous definition of meson masses and coupling constants in dense matter. 

At the level of QCD, there is a quark condensate, that is, a vacuum expectation value 
(0|gg|0), whose presence signals the spontaneous breakdown of chiral symmetry and leads to 
a baryon mass. QCD sum rules combine hadronic amplitudes, quark and gluon condensates, 
and asymptotic freedom to obtain constraints on the properties of hadrons. This quark 
condensate will be modified in the presence of other baryons; using sum rules, one can relate 
the change in baryon (or more generally, hadron) masses in nuclear matter to the underlying 
QCD. Important papers in this regard include ||Co91a| , ^I^ , |Ce92a| , |Co92a| , |Fu92a|JCe93|,|Ji93a| . 


In |[E191|| , chiral and scale symmetries of QCD are used to describe the interaction between 
these condensates and hadrons. The resulting equations are solved self-consistently in the 
RMFT approximation. For these QCD condensates to be driven towards zero at high density. 
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their coupling to scalar and vector mesons must be such that the masses of these mesons 
do not decrease with density. In this case, a physically sensible phase transition to quark 
matter ensues. 

In ||Co91a| , po92a| , |Fu92a|| , the self-energies for quasinucleon states in nuclear matter are 
studied using QCD sum rules. A correlator of nucleon interpolating helds, evaluated in 
the hnite-density ground state, is calculated using both an operator product expansion and 
a dispersion relation with a spectral Ansatz. This approach relates the nucleon spectral 
properties (such as the quasinucleon self-energies) to matrix elements of QCD composite 
operators (condensates). With increasing nucleon density, large changes in Lorentz scalar 
and vector self-energies arise naturally; the self-energies are found to be comparable to 
those arising in RMFT phenomenology. The most important phenomenological inputs are 
the baryon density and the value of the nucleon “sigma term” (un) divided by the average 
current mass of the light quarks. The successful comparison to RMFT phenomenology 
is, however, sensitive to assumptions about the density dependence of certain four-quark 
condensates ||Ji93a[pT9^ . 

One extension of QHD involves developing model held theories where the mesons couple 
directly to quarks. Three important papers in this regard are [|Ja92| , pa94| , ^a94b| | . In ||Ja92|| , 


chiral symmetry restoration at hnite baryon density is studied in a quark model involving 
both scalar and vector interactions. The presence of vector interactions makes chiral sym¬ 
metry restoration more difficult. On-shell masses and coupling constants are calculated for 
the Tz,u,p, and ai mesons. An attempt is made to relate the quark-meson interactions to 
the observed nucleon-meson coupling constants. 

In summary, relativistic quantum held theories of the many-hadron system, with local 
lagrangian densities, appear to provide the appropriate framework for extrapolating the 
properties of ordinary hadronic nuclear matter to all temperatures and densities up to the 
transition to the quark-gluon phase. We can understand this based on our earlier discussion. 
Since the truncation of an ehective lagrangian involves an expansion in powers of fcp/M, 
then even at 10 times equilibrium density, this ratio is only 0.6. Moreover, the relevant 
temperatures are T < 200 MeV, which are small compared to the baryon and heavy meson 
masses. Nevertheless, interaction terms that are small at normal densities and tempera¬ 
tures, and thus difficult to calibrate, can become important at densities and temperatures 
relevant to the phase transition, thus producing signihcant uncertainties in the extrapolation 
Mu96|jMu96a|| . Finding appropriate ways to calibrate these interactions is a major challenge 


for future investigations. 


8. LOOPS IN QHD 

We have seen that the mean-field approximation to QHD gives a concise and highly suc¬ 
cessful nuclear phenomenology. We have also noted, based on ideas from density functional 
theory, that a mean-field energy functional fitted to nuclear properties implicitly includes 
some effects of higher-order, many-body corrections. Nevertheless, mean helds are insuf- 
hcient for a detailed understanding of nuclear structure and its relation to the underlying 
NN, NNN, ... interactions. Moreover, various experimental observables probe aspects of 
nuclear dynamics that go beyond a single-particle description, which can tell us about the 
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modification of the NN interaction inside a nucleus. We must therefore develop reliable, 
relativistic techniques to extend the RMFT calculations discussed above. Useful tools in 
this endeavor are Feynman diagrams, dispersion relations, and path integrals, as discussed 
where a historical development is presented, and in ||Se92| 


m 


Se86 


which updates the 

1986 volume. References to the original literature are given in these review articles. Here 
we briefly summarize some of those results and discuss some of the work that has been done 
since that time. A more detailed discussion of the theoretical background is contained in 
the recent text [|Wa95| . 

Corrections to the RMFT will generically involve “loop diagrams”, which are important 
for including several different aspects of the quantum nature of the system. For example, loop 
diagrams must be included to ensure the unitarity of scattering amplitudes. Baryon loops 
are necessary for incorporating familiar many-body effects, like the exchange of identical 
nucleons or the summation of ladder and ring diagrams |[Fe71|| . Moreover, loops introduce 


effects arising from the modihcation of the quantum vacuum in the presence of valence 
nucleons, and meson loops in particular generate contributions to the extended structure of 
the nucleon. The relevant question is how to treat these loops most efficiently, consistent with 
the notion that QHD is intended to be a large-distance, hadronic model of the underlying 
QCD. We begin the discussion by considering pion loops and then turn to the more difficult 
question of heavy-meson and baryon loops. 


A. Pion Loops 


Because the pion has a small mass, pion loops make signihcant contributions to both 
scattering and nuclear-structure observables. Moreover, since the energies and momenta of 
interest are typically of the order of the pion mass (or even larger), pion-loop contributions 
generally involve nonanalytic functions of the external four-momenta. Thus these contribu¬ 
tions cannot be absorbed in a local effective lagrangian, which, by dehnition, contains only 


hnite powers of derivatives 


They must be computed explicitly, whether one 


is using a renormalizable model or an effective, nonrenormalizable lagrangian. 

A highly developed framework for systematically including pion loops in scattering pro¬ 
cesses is chiral perturbation theory (ChPT) ||Ga84|,^e93a|| . Chiral symmetry implies that 


pion interactions in the lagrangian can be grouped order-by-order in the number of deriva¬ 
tives and powers of (as in the nonlinear models discussed earlier), so that there is a 
systematic expansion for scattering amplitudes at low energies. Moreover, an expansion in 
the number of loops also proceeds in powers of momenta |[We67| , |We90a|] , so one can sys¬ 
tematically include loop corrections. ChPT has been successful in describing scattering in 
both the R = 0 and B = 1 sectors of low-energy QCD ||Ga84| , ^e88| , |Me93a[] . Studies of two- 
and many-nucleon systems have been initiated and are currently under active investigation 
|0792| , [Va9^ , P796| . 

Properties of the baryon also arise through pion-loop integrals in QHD. For example, the 
vertex diagram consisting of the emission of a pion, its interaction with the virtual electro¬ 
magnetic held, and its reabsorption by the baryon contributes to the baryon’s anomalous 
magnetic moment. The two-pion contribution gives the low-mass, or long-distance, part of 
the spectral weight function P 2 {cr‘^) for the anomalous magnetic form factor F 2 {q‘^): 
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Assume that the two-pion contribution arising from this vertex diagram dominates the 
spectral weight function everywhere. This yields a qualitative description of the isovector 
anomalous magnetic moment and its mean-square radius |^e86|JWa95|| .P] 

Pion loops also give hadronic contributions to vacuum polarization. For example, a spec¬ 
tral analysis of the strong-interaction contribution to electromagnetic vacuum polarization 
shows that the spectral weight function starts at 4m^: 


n”(g) = {,q^,qy - q ) , 


W) = - 


p{a^ 


TT JAml — q^ 


da' 


( 8 , 2 ) 

(8.3) 


In the complex plane, n(g^) is an analytic function with a branch cut running along the 
real axis from 4m^ to infinity. The discontinuity across that cut for 4m^ < O'' < 9m^ comes 
from the electroproduction of two real pions. Thus the low-mass singularities of propagator 
and vertex functions are most efficiently expressed in terms of hadronic variables. The 
higher-mass singularities are more complicated, and hadrons are less efficient. Nevertheless, 
if one is interested in studying the low-momentum behavior of the vacuum polarization, one 
can emphasize the low-mass part of the spectral integral by making several subtractions 
from the integral and by determining the unknown coefficients empirically ||Do96| . 


The exchange of two correlated pions between nucleons also involves multiple pion loops, 
which can be treated with dispersion relations ||Ja75| , Pu77| , Pu80| , |Li89| , |Li90| , |Ki94|| . The resnlt 
is a strong NN attraction that can be simulated by introducing an effective scalar-isoscalar 

and B.D. 


field coupled directly to the nucleon, as discussed in Sections 


More detailed 


studies of the scalar-isoscalar NN interaction can be performed by returning to the descrip¬ 
tion with explicit pion loops ||Du93| , [Ao95|| . Similar observations can be made in exchange 
channels with vector meson quantum numbers. In particular, a tree-level effective lagrangian 
with vector mesons is found to be essentially eqnivalent to ChPT in the pion sector at one- 
pion-loop order | Me88 , rTo89 , Pi)c89 , |Ec89a |. Thus these non-Goldstone (heavy) bosons can 
be introduced into an effective lagrangian containing baryons to efficiently acconnt for the 
intermediate-range NN interactions and to conveniently describe nonvanishing expectation 
valnes of nuclear bilinears {e.g., NN and N'j^N), as in Section ES 

Moreover, as noted in Section p.C| , the sum of ttN ladder diagrams with nucleon ex¬ 
change can be investigated with partial-wave dispersion relations, leading to a resonance in 
the A channel. Thus the A(1232) can be included in an effective hadronic lagrangian to 
incorporate this important physics ||De92a|JWe93| , P?a96| . Alternatively, the dynamical model 
of the resonance provides a means to investigate many interesting questions concerning the 
behavior of the A in the nuclear many-body system, such as its binding energy in nuclear 
matter, its optical potential, and the modification of its electroweak properties. 


^^The isoscalar anomalous moment vanishes in this approximation; experimentally, it is indeed 
very small. 
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Virtual pion effects have been studied extensively, 
examined in 


Pa92|| using a coherent-state approach. Reference 


Pion loops and pion dressing are 


studies in depth the 


contribution of the virtual pion cloud to the electromagnetic properties of the nucleon. 


B. Loops in Renormalizable QHD Theories 


By starting with the Feynman rules for the Green’s functions in a renormalizable QHD 
theory, one can in principle go beyond the RMFT, compute observables in terms of a hnite 
number of couplings and masses, and then compare with experiment. In practice, this 
program is extremely difficult, since QHD is a strong-coupling relativistic quantum held 
theory. Moreover, renormalizable theories explicitly include contributions from all length 
scales, and the QHD couplings get stronger at short distances, since these theories are not 
asymptotically free. 

Nevertheless, just as in nonrelativistic many-body theory, one can use intuition to sum 
selected inhnite sets of diagrams, determine the renormalized coupling constants by rehtting 
nuclear matter properties, and then see whether the RMFT results are stable under the 
inclusion of these additional contributions, while investigating new physical phenomena. All 
of the extensions we discuss involve loop corrections to the RMFT of one sort or another. 
These corrections include familiar many-body effects, where the loop momenta are typically 
of the size of (or at most, several times /cp), and also corrections from the dynamical 
quantum vacuum, which involve shorter distance scales. 

In the RMFT, the baryon Green’s function can be written as 


Se86 


GW = (7,r'‘ + M-){p^ 

= Geik) + Goik) . 


+ 


+ ie E 


^^s[ki-E'{kmk,-M)] 


(8.4) 


in the rest frame of the nuclear matter. Here k*^ = (fco — QvVo, k) is the kinetic four- 
momentumQ and E*{k) = y/k^ -|- M*^, with M* = M — 5's</>o- The hrst term Gpik) is the 
Feynman propagator for a baryon of mass M*, and the second term is the contribution 
arising from baryons already present at hnite density; this latter contribution reproduces 
the RMFT results. In discussing the following extensions, we shall frequently distinguish 
results obtained with the full baryon propagator G{k) from those obtained with just the 
second, or “density-dependent”, contribution Goik). Since the three- and four-momenta are 
constrained in Goik), loop integrals over this second term give well-dehned, hnite results 
that are direct analogues of the terms arising in nonrelativistic many-body theory. 

For example, relativistic Hartree theory is obtained by self-consistently summing the 
tadpole graphs in the baryon self-energy |^e86|| . Retention of in the tadpoles gives rise 
to the MFT, while the full G (together with appropriate counterterms SC) produces the 


“^^Note that in closed-loop integrals, such as those involved in computing the ground-state energy, 
a simple shift of integration variables allows one to eliminate the dependence on ^vH) ||Se86,Fu89|. 
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RHA, as discussed previously.P] A characteristic result of QHD-I is that the Lorentz scalar 
and vector self-energies are large on the scale of nuclear energies; these contributions cancel 
in the binding energy but add in the spin-orbit interaction. 

Hartree-Fock (HF) theory is obtained by including the meson emission and reabsorption 
(“exchange”) graphs in the baryon proper self-energy ||Se86| , pe^ . Retention of Gd only 
in this HF theory leads to exchange terms that are the direct relativistic generalization 
of those arising when Slater determinants are used to determine the best single-particle 
wave functions in nonrelativistic many-body theory. The inclusion of the exchange terms 
does not qualitatively alter the size of the large scalar and vector self-energies found in 
the RMFT, and thus the RMFT is stable with respect to these corrections. In fact, after 
renormalization to equilibrium nuclear matter properties, the binding energy curves in the 
RMFT and HF approximations are almost indistinguishable ||Se86|| . Moreover, the equation 
of state approaches that of the RMFT at high baryon density. Further discussion is contained 
in 

Relativistic HF calculations that include charged mesons and that examine the relation 
of the scalar and vector self-energies to DBHF calculations in nuclear matter are carried out 
Bo87|,|Ce9^Pe93|,|Fr93a|JBo94a|, Fr94a|,^a94c|,|Zh94a[Pe96| 


m 


The existence of large Lorentz 
scalar and vector self-energies appears to be a hrm conclusion. (Recall our discussion of 


DBHF results in Section m 

voia 


The connection to semiclassical approximations is given in 
The general structure of the self-energy in matter is analyzed in ||He91|jRu95a| . 


The fully self-consistent HF theory retaining the complete G{k) and meson retardation 
is complicated |pe92|| ; it has not yet been successfully solved. The summation of exchange 


diagrams only (he., HF theory at zero baryon density) is discussed in 

The sum of fermion ring diagrams involving G{k) is equivalent to the random-phase 
approximation. We shall refer to the calculation of the rings that keeps only terms with at 
least one factor of Gd as the RPA. This includes loops that are at least linear in the density 
(particle-hole parts and admixtures between hlled valence states and those in the Dirac 
sea) and is the direct relativistic extension of the RPA in nonrelativistic many-body theory. 
The calculation that also includes the modihcation of the strong vacuum polarization in the 
nuclear medium due to the shift M —> M* will be called RRPA. 

We emphasized the need 


The RPA as applied to hnite nuclei was discussed in Section 2.D 


to include contributions from negative-energy states in order to maintain the conservation 
laws in the theory. Some important results from RPA studies of nuclear matter are as follows. 
The scalar and vector propagators mix in nuclear matter, and at high density, vector meson 
exchange dominates in QHD-I. The excitation spectrum of nuclear matter in the RPA is 
that of zero sound, where the sound velocity Cq approaches the speed of light from below 
as the baryon density gets large (cq —> 1“ as kp —>■ oo). This implies that signals in the 
medium cannot propagate faster than the speed of light, in accord with special relativity. 
There are other branches in the excitation spectrum corresponding to meson propagation. 

We turn now to calculations involving strong vacuum polarization (RRPA). This polar- 


^^One can also solve H-l dimensional QHD exactly for finite systems in the one-loop approxima¬ 
tion; this is carried out in |Fe93|. 
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ization does not have any explicit density dependence, but depends on it implicitly through 
the baryon mass M*. These vacuum loops correspond to the one-baryon loops contained in 
the RHA (Section ^.Ep , and just as we found unnaturally large contributions in the RHA, 
analogous large effects appear in the RRPA. In consequence, in nuclear matter, Landau 
“ghost” poles appear in the meson propagators at zero frequency go = 0 and hnite wave 
number |q| 7 ^ 0; the value of this wave number is a few times the nucleon mass in QHD-I 
||(Jo87| , Pe57| , |Tu88|jWe9U| , |Li9Ua)| . Such poles imply an instability of the system against density 
fluctuations of the corresponding wavelength and are a manifestation of the lack of asymp¬ 
totic freedom in the theory. Apparently, a description of the quantum vacuum by summing 
simple baryon loops is inadequate. 

Similar problems occur when one extends the RHA to include two-loop contributions 
[hu89|| , which consist of a closed baryon loop with an internal meson line. When the full 
propagator G is used in this calculation, one incorporates the many-body modifications to 
both the strong vacuum polarization and the baryon self-energy that arise from the Pauli 
exclusion principle and from the shifted baryon mass M*. One again finds unnaturally large 
contributions to the nuclear matter energy. 

The problem in all of these calculations of vacuum effects is that an expansion in powers 
of loops is basically an expansion in the dimensionless coupling constants, which are large 
in QHD. The quantum corrections are correspondingly large, the series is not converging, 
and the RMFT is not stable against this perturbative loop expansion. Clearly, an alternative 
procedure must be found to systematically and reliably calculate the vacuum-fluctuation 
corrections to the RMFT in renormalizable QHD, if that is indeed possible. 

One idea is to include corrections that sum diagrams to all orders in loops, such as 
the “ladder” summation discussed in Section |4-B| . Alternatively, one observes that in a 
theory with a vector boson coupled to a fermion, the vertex form factor is highly damped 
at large spacelike momentum transfers < 0 ; this is analogous to the result hrst derived 
by Sudakov in quantum electrodynamics |[[Vfi91||. Such a form factor would decrease the 


sensitivity to high-momentum (or short-distance) contributions to loop integrals and would 
provide a favorable situation for QHD. Moreover, it is essential to include vertex corrections, 
as these reflect the internal hadron structure present in renormalizable QHD. Remember that 
baryons are complicated objects in the full held theory; they are surrounded by a cloud of 
virtual mesons and baryon-antibaryon pairs. It was an initial hope that the theory might 
contain vertex functions that would damp the contributions from loop integrals before one 
reaches distance scales where the lack of asymptotic freedom becomes manifest.^] 

In ||A192|| , vacuum polarization is studied in QHD-I (without the scalar meson). The 
lowest-order (one-baryon-loop) polarization produces a “ghost” pole when summed to all 
orders in the RRPA, as discussed above. It is hrst verihed that the infrared structure of 
the meson-baryon vertex in this model produces an on-shell proper vertex function that is 
strongly damped at large spacelike momentum transfer. When the model vertex function 
is approximated by its on-shell form and combined with the lowest-order polarization, the 


^^Note that in a theory with a vector field coupled to a conserved fermion current, Ward’s identity 
implies that the structure of the vertex function is unrelated to the running of the coupling constant. 
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vacuum contributions are significantly reduced. The resulting RRPA meson propagator has 
no “ghost” poles and is finite at large spacelike momenta. 

In |^e95|| , it is shown that a similarly damped form factor arises for the scalar-nucleon 
vertex in QHD-I due to the vector-meson dressing. The on-shell approximations to the 
vertex functions are used to investigate the two-loop contributions to the properties of 
nuclear matter, and it is found that they are greatly reduced by the inclusion of vertex 
functions. Similar results are found using ad hoc form factors in |[Fr92| , [Fr9^ . 

RRPA calculations are performed with ad hoc vertex cutoffs in |Pr92a)| and with the 
on-shell model described above in ||la93a |. A two-loop calculation that also consistently 
treats the electromagnetic interaction to the same level of approximation is developed in 
Extensions of RPA and RRPA to include the self-consistent sum of ring diagrams 

Be93a|| 


Be92 


with additional A-hole propagation are examined in 


The relation of the density 


dependence of zero sound to the renormalization prescription is discussed in ||Ca95| . 

The fully off-shell vertex function is complicated in any field theory. An attempt to 


include the off-shell vertex in the RRPA calculation for QHD-I has been initiated in ||A195 
no concrete results exist yet for this very difficult calculation. A simplified approach to the 
baryon self-energy, using an ad hoc off-shell form factor and dispersion relations, is contained 
in [^fg. 

Nuclear Schwinger-Dyson equations, which provide a basis for an analysis of the rel¬ 
ativistic field theory content of the nuclear many-body problem in terms of propagators 
and vertices, are developed in |[Na91| , |<o9^ , [Na94|,|Na94aj| . Truncated applications usually use 
some sort of spectral representation or dispersion relation [P?a91|| to eliminate the Landau 
“ghost” poles in the meson propagators discussed earlier. 

Pion propagation in the nuclear medium, which involves baryon-antibaryon loops, is 
studied in [|He92a||Ra92|| ; again, the “ghost” poles are generally removed through one of the 
mechanisms described above. Vector meson propagation is examined in ||Je94|| and vector 
meson mixing through baryon loops is examined in |[Pi93| , |Pi93a|| . (Note that here the strength 
of the mixing is finite without any renormalization.) Meson modes in nuclear matter are 
also examined in ||Ja93| , paM^ . 

A chirally symmetric, renormalizable, anomaly-free theory that contains baryons and 
(tt, cr, a;, p, ai) mesons is developed in |^e92b|| ; we refer to this theory as QHD-III. This 


model provides a hadronic description of strongly interacting matter that includes isovector, 
pseudoscalar and vector fields in addition to the isoscalar, scalar and vector fields of QHD- 
I. Although this theory provides a consistent, self-contained description of nuclear physics, 
including loop processes, its phenomenology remains to be investigated |[Wa95| . 


To summarize, the systematic calculation of loop corrections to the RMFT is an impor¬ 
tant goal in renormalizable QHD. At present, vacuum contributions evaluated at various 
orders in loops are all unnaturally large, signaling the inadequacy of these approximations 
for representing the vacuum dynamics. Although these large effects can be reduced by in¬ 
troducing vertex form factors external to the renormalizable theory, the question of whether 
such form factors can be generated internally, and indeed, whether calculations with off- 
shell vertices are even feasible, is still open. At best, even if renormalizable QHD provides 
a realistic description of vacuum dynamics, the results will involve sensitive cancellations 
between unnaturally large contributions. We are thus motivated to consider other options. 





























































C. Loops in Effective QHD Theories 


The results discussed above present significant evidence that the requirement of renor- 
malizability is too restrictive, and that it is more appropriate to consider QHD as an ef¬ 
fective, nonrenormalizable hadronic field theory. In this approach, short-distance effects 
and vacuum modifications are parametrized through nonrenormalizable interaction terms in 
the lagrangian, with couplings that are determined by fitting to data.0 Nucleon and non- 
Goldstone boson fields are still needed in the lagrangian to account for the valence nucleons 
and to treat the large mean fields conveniently. Loop diagrams involving these degrees of 
freedom must also be included to ensure the unitarity of scattering amplitudes and to in¬ 
corporate many-body, density-dependent effects. The relevant question is how to separate 
the short-range contributions, which are to be absorbed in the model parameters, from the 
long-range, many-body effects. 

The baryon propagator in Eq. ( |8.4|) separates into an explicitly density-dependent part 
Gd and a part Gp that has no explicit density dependence. The familiar many-body contri¬ 
butions arise from loop integrals involving Gp alone (as in the RMFT or HE approximations) 
or combinations of Gp and Gp (as in the DBHF approximation or the RPA); these are to 
be computed explicitly. One must also retain various density-dependent contributions that 
contain negative-energy states, but that are finite, such as those needed to preserve the 
conservation laws in the RPA. The remaining contributions involve divergent loop integrals 
that contain vacuum dynamics and short-distance behavior; these are to be absorbed in the 
parameters of the lagrangian that are fitted to empirical data. Some of these contributions 
will have explicit density dependence (like the strong Lamb shift), while some will not (vac¬ 
uum polarization). We emphasize that this is a more difficult problem than the inclusion of 
loops in ChPT, since the loop corrections do not produce an expansion in powers of external 
momenta, and thus corrections at a given order in loops renormalize the parameters at all 
lower orders in loops. 

In a renormalizable theory, the divergences can be removed by cancelling them against 
a finite number of counterterms. In a nonrenormalizable theory, however, the number of 
divergences is, in principle, infinite. Nevertheless, these divergences can always be formally 
cancelled by the appropriate counterterms, since all possible interaction terms consistent 
with the symmetries are already included in the effective lagrangian. For example, in infinite 
nuclear matter, any vacuum contribution to the energy arising from loops containing only 
factors of Gp must appear in the form of a polynomial in the scalar field 0o-0 As another 
example, contributions from the nucleon’s Lamb shift can be (formally) cancelled by coun¬ 
terterms of the form or N{i'^^d^ — The important issue thus becomes 

how to limit the required counterterms to a finite and manageable number. Moreover, due 


^^These couplings also implicitly include the effects of heavier degrees of freedom not included in 
the model lagrangian. 

Although nonanalytic functions of (f)Q may arise in principle, we assume that any function 
involving heavy mesons can be expanded in a Taylor series in the relevant density regime. 
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to the freedom one has to redehne the held variables, the most efficient way to write the 
connterterms is still an unsolved problem. 

One way to constrain the parameters is to rely on the broken scale invariance of QCD, 
which restricts the interactions in the scalar-isoscalar channel. In ||Fu95|| , an effective rela¬ 


tivistic hadronic model for nuclear matter that incorporates nonlinear chiral symmetry and 
broken scale invariance is applied at the one-baryon-loop level. The model contains an ef¬ 
fective light scalar held that is responsible for the mid-range NN attraction and that has 
anomalous scaling behavior. One-loop vacuum contributions in this background scalar held 
at hnite density are constrained by low-energy theorems that rehect the broken scale invari¬ 
ance of QCD, so that the scalar ehective potential contains only three free parameters. The 
resulting mean-held energy functional for nuclear matter and nuclei has only a hnite number 
of parameters and yields good hts to the bulk and single-particle properties of nuclei. 

These results are consistent with the discussions in Sections || and |, where it was ob¬ 
served that only the hrst few terms in an expansion in powers of the helds (and their 
derivatives) are relevant up to moderate densities; in particular, only three parameters are 
needed in the scalar potential [m^, ks, and K 4 in Eq. ( |5.47|) ]. This result obtains because 
the empirically htted parameters are of natural size. Thus one may expect that the as¬ 
sumption of naturalness is sufficient to limit the unknown parameters in the lagrangian to 
a manageable number. This expectation remains to be conhrmed by calculations beyond 
one-baryon-loop order, which are currently in progress. 

To summarize this section on loop processes in QHD: although there are many qualita¬ 
tive and even semi-quantitative insights, and applications of the ideas of nonrenormalizable 
effective held theory look promising, there is still no consistent, reliable, practical approach 
to the relativistic nuclear many-hody problem that includes all loop terms. 


9. SUMMARY 

Our goals in this paper are to describe consistent microscopic treatments of the relativistic 
nuclear many-body problem based on hadrons (quantum hadrodynamics) and to summarize 
work in this held from early 1992 through 1995. Although QCD is the underlying theory 
of the strong interaction, the QCD couplings are strong at large distances, and hadrons are 
more efficient degrees of freedom. But since it is still impossible to derive the low-energy 
hadronic lagrangian directly from QCD, we must rely on more indirect information, such 
as the symmetries of the QCD lagrangian, and on well-known nuclear phenomenology to 
guide us in the construction of the low-energy theory. The framework is based on local, 
Lorentz-invariant lagrangian densities, as this is the most general way to parametrize ob¬ 
servables consistent with the desired constraints of quantum mechanics, special relativity, 
unitarity, causality, cluster decomposition, and the intrinsic QCD symmetries. Historically, 
the hadronic lagrangian was required to be renormalizable, so that one could calibrate the 
theory and then extrapolate without the appearance of new, unknown parameters. There are 
now strong indications, however, that the constraint of renormalizability is too restrictive, 
so it is important to generalize our viewpoint to include nonrenormalizable, effective held 
theories, which can still provide a consistent treatment of the relativistic nuclear many-body 
problem. 







The simplest model that we study (QHD-I) contains protons, neutrons, and neutral 
scalar and vector mesons, and is renormalizable. At large enough densities, the meson held 
operators can be replaced by their expectation values, and the result is a relativistic mean- 
held theory (RMFT), which must be solved self-consistently. In the original version of the 
model, which omits scalar self-couplings, one obtains a simple, two-parameter description of 
the equilibrium properties of nuclear matter that can be extrapolated to arbitrary density, 
temperature, and proton fraction. The classical meson helds are large (several hundred 
MeV), and nuclear saturation occurs because the nucleon ehective mass M* decreases as 
the density increases, so that the attractive forces saturate, and the binding curve develops 
a minimum. If this RMFT is applied to hnite nuclei (with the addition of classical helds for 
the neutral rho meson and Coulomb potential), then from this minimal set of parameters 
htted to the properties of nuclear matter, one derives the nuclear shell model. Moreover, 
when cubic and quartic scalar self-interactions are included, one can tune the equilibrium 
nuclear matter properties to provide a realistic description of nuclear charge densities, charge 
radii, binding energies, single-particle spectra, and quadrupole deformations throughout the 
periodic table. 

To explicitly include pions, we need a lagrangian consistent with global chiral S'[/(2 )l x 
SU{2)ii symmetry, which is a symmetry of QCD in the limit of massless u and d quarks. It 
is simplest to use a linear representation of the symmetry, in which the helds enter as chiral 
multiplets. If one also demands renormalizability, one is led to the linear sigma model, with 
an additional neutral vector held. The neutral scalar held plays a dual role as the chiral 
partner of the pion and the mediator of the attractive NN force. One hnds, however, that 
the mean-held approximation to this model is unable to provide a realistic description of 
nuclei, primarily due to the strong constraints on the scalar self-interactions arising from 
the scalar potential used to induce spontaneous symmetry breaking. This is evidence, even 
at the level of the RMFT, that the simultaneous constraints of linear chiral symmetry and 
renormalizability are too restrictive. 

We observe, however, that it is possible to make a chiral, point transformation of the 
helds and to recast the lagrangian so that the symmetry is realized nonlinearly. The new 
scalar held is a chiral singlet that can be decoupled from the theory by taking its mass to 
be very large, thus removing the unwanted nonlinear interactions. (In the limit of inhnite 
scalar mass, the model becomes nonrenormalizable.) Although the mid-range NN attraction 
is apparently destroyed by this procedure, it can be restored by including correlated two- 
pion exchange in the scalar-isoscalar channel. Moreover, this two-pion exchange can be 
efficiently and adequately simulated by introducing a new scalar held into the theory with a 
mass of roughly 500 MeV.Q The new held is also a chiral scalar, so the lagrangian remains 
chirally invariant, and the RMFT of this chiral model is identical to that of QHD-I. Thus 
we draw the important conclusion that the scalar held in QHD-I is to be interpreted as an 
effective held that incorporates the pion-exchange contributions that are the most important 


observe with some amusement that the possibility of introducing such an additional scalar 
into a chiral theory was already noted in the early work of Coleman, Wess, and Zumino on nonlinear 
realizations of chiral symmetry. See footnote 9 in [Po69|. 






for describing the bulk properties of nuclear matter. 

To get a deeper understanding of the large mean fields and to connect the RMFT de¬ 
scription to the underlying NN interaction, we consider corrections from nucleon exchange 
and short-range correlations in nuclear matter. This can be done in the Dirac-Brueckner- 
Hartree-Fock (DBHF) framework. Here an NN quasipotential htted to two-nucleon data is 
used to determine a self-consistent NN interaction in the medium, which includes the effects 
of the Pauli exclusion principle and the modihcations to the single-nucleon Dirac wave func¬ 
tions at hnite density. Three important conclusions from this work are: the nucleon scalar 
and vector self-energies (which are analogous to the scalar and vector helds in the RMFT) 
are essentially the same size as the scalar and vector mean helds studied earlier, and the state 
dependence of the self-energies is small; the depletion of the Fermi sea due to correlations 
is considerably smaller in the DBHF theory than in nonrelativistic Brueckner-Goldstone 
theory; the DBHF effective NN interaction contains density dependence that goes beyond 
what is included in nonrelativistic Brueckner-Goldstone theory, and it is therefore possible 
to simultaneously ht both the NN phase shifts and the nuclear matter equilibrium point 
at the two-hole-line level. This last result can be understood by considering an expansion 
of the RMFT energy/particle in powers of the Fermi momentum; one hnds that even two- 
body interactions involving a Lorentz scalar held lead to terms that would be interpreted as 
many-body forces in a nonrelativistic framework. 

Motivated by the appearance of an ehective scalar held in our discussion of chiral symme¬ 
try, we generalize the QHD framework to embrace the ideas of nonrenormalizable, ehective 
held theory. We still retain a local, Lorentz-invariant lagrangian, but now the lagrangian 
must contain all possible (non-redundant) interaction terms consistent with the symmetries 
of the underlying QGD. The coefficients of these terms parametrize the short-distance dy¬ 
namics that will be modelled by the ehective hadronic lagrangian. Since we cannot calculate 
these coefficients directly from QGD, they must be determined by htting to data, and we 
have the important freedom to choose our generalized coordinates (helds) to make this ht¬ 
ting as efficient as possible. Moreover, since the ehective lagrangian has, in principle, an 
inhnite number of terms, we must have an organizing principle to retain some predictive 
power, and this is also inhuenced by the choice of dynamical variables. 

Based on the successes of the hadronic description of both the NN interaction and the 
nuclear-structure observables, we choose as degrees of freedom the nucleon, pion, and low- 
mass scalar and vector helds. By retaining the “heavy” non-Goldstone bosons, we can 
describe the NN interaction without the explicit calculation of multi-pion loops, and we 
can efficiently describe the expectation values of nucleon bilinears using mean helds. We 
then construct a lagrangian consistent with the underlying symmetries of QGD, in which 
the chiral symmetry is realized nonlinearly and electromagnetic interactions are introduced 
both through minimal-coupling terms and through derivative couplings of the nucleon to 
the photon, which allows us to incorporate the nucleon electromagnetic structure directly in 
the theory. 

To organize the lagrangian, we rely on naive dimensional analysis, which allows us to 
extract the dimensional scales of any term, on the assumption of naturalness, which says that 
the remaining dimensionless coefficient for each term should be of order unity, and on the 
observation that the ratios of mean meson helds gs4>o and g^Vo (and their derivatives) to the 
nucleon mass M (which generically represents the “heavy” mass scale) are good expansion 


90 



parameters. Since the meson fields are ronghly proportional to the nnclear density, and 
since the spatial variations in nnclei are determined by the momentnm distributions of the 
valence-nucleon wave functions, this organizational scheme is essentially an expansion in 
kp/M, where fcp is a Fermi wavenumber corresponding to ordinary nuclear densities. If 
naive dimensional analysis and the naturalness assumption are valid, we can expand the 
lagrangian in powers of the fields and their derivatives, truncate at some finite order, and 
thus have some predictive power for the properties of nuclei. The model parameters fitted 
to bulk and single-particle nuclear properties show that this is indeed the case. 

We also observe, based on concepts from density functional theory, that RMFT pa¬ 
rameters fitted to nuclear properties implicitly include not only the short-distance effects 
parametrized in the lagrangian, but also long-range, many-body effects arising from cor¬ 
rections to the RMFT. For evidence that mean fields should provide a reasonably accurate 
way to parametrize these many-body effects, we rely on the DBHF results, which show that 
exchange and short-range-correlation corrections do not significantly modify the mean fields 
obtained at the RMFT level. Thus we conclude that densities “large enough” to justify 
the use of mean-field theory are present in the interiors of medium and heavy nuclei, since 
the Hartree contributions dominate the nucleon self-energies already at these densities. The 
numerous successful RMFT calculations of nuclei throughout the periodic table should be 
interpreted within this density functional context. By explicitly calculating higher-order 
corrections, one can improve the description of the density functional, and by using the 
properties of nuclei to determine the parameters, one will ultimately obtain the values that 
should appear in the effective lagrangian. 

This observation leads to another important conclusion. No matter how well one can 
compute the many-body contributions, there will always be additional terms in the en¬ 
ergy functional that depend (for example) on powers of the scalar field (or scalar density). 
The coefficients of these terms, which parametrize the short-distance behavior of QCD, are 
unknown and are presently impossible to calculate. Moreover, these contributions will be 
significant on the scale of the nuclear binding energy when their coefficients are of natural 
size. Unless one has a concrete argument for why these coefficients should be anomalously 
small, they will always be present to adjust the energy functional in a way that allows one 
to fit nuclear matter saturation or nuclear properties. The goal of predicting the nuclear 
matter equilibrium point from an NN interaction fitted to two-body data therefore loses 
some of its importance. In the conventional nonrelativistic language, these additional terms 
correspond to the short-range parts of the three-body (and many-body) forces that cannot 
be calculated and that must be fitted to the desired observables.^ 

The problems of extrapolation to extreme conditions of density and temperature, or 
to large energy-momentum transfers, and the systematic computation of loop corrections 
still pose challenges for QHD. Although one can argue that the effective lagrangian can be 
truncated at a few terms for studies of neutron stars or the transition to the quark-gluon 
plasma, interaction terms that are small at normal densities and temperatures, and thus 


^®The situation has been succinctly described by Jackson | Ja92a |: 
too little dirt and too much rug under which to sweep it.” 


“The problem is that there is 
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difficult to calibrate, can become important under extreme conditions. This incomplete 
calibration leads to significant uncertainties in the extrapolation, even for parameters of 
natural size. As for loop calculations, there is now signihcant evidence that it is insufficient 
to represent vacuum modifications by simple loop integrals in renormalizable QHD theories. 
Although it may be possible to reduce the size of these corrections by incorporating vertex 
modihcations, calculations involving off-shell vertex functions computed entirely within the 
context of the theory remain to be done. 

Nonrenormalizable effective theories appear more promising, since short-range and vac¬ 
uum effects are absorbed into the parameters of the lagrangian, and long-range, many-body 
effects are (typically) straightforward relativistic generalizations of their nonrelativistic coun¬ 
terparts.^ Moreover, calculations of pion loops based on chiral perturbation theory or dis¬ 
persion relations can be performed systematically. Nevertheless, a fully consistent, practical, 
relativistic many-body description of nuclei based on effective hadronic field theory remains 
to be formulated. 

In summary, nuclear physics is the study of the structure of strongly interacting baryonic 
matter {B > 1), and the only consistent theoretical framework we have for describing such 
a relativistic, interacting, quantum-mechanical, many-body system is relativistic quantum 
held theory based on a local lagrangian density. Although QCD of quarks and gluons pro¬ 
vides the basic underlying theory, lagrangians based on hadronic degrees of freedom (QHD), 
which are the particles observed in the laboratory, provide the most efficient description 
of the physics in the strong-coupling, nuclear domain. While it is not surprising that ini¬ 
tial QHD attempts to model the system within the subset of renormalizable lagrangians 
appears to be too restrictive, some results of the initial simple models remain robust. In 
particular, the notion of strong isoscalar, Lorentz scalar and vector mean helds remains 
valid in a density-functional framework where correlations and other higher-order effects 
are incorporated in the functional. In the modern effective held theory approach to QCD, 
one incorporates only the underlying symmetries of QCD in the hadronic lagrangian. With 
the crucial observation that while the scalar and vector helds in nuclei are comparable to 
the rest mass of the nucleon, the ratios gs(j)o{x)/M and g^VQ{x)/M still provide expansion 
parameters less than unity, one is able to understand the multitude of successful applica¬ 
tions of the relativistic mean-held treatment of nuclei (RMFT). Applications of QHD to 
nuclear structure, electroweak interactions with nuclei, the hadronic region of the nuclear 
phase diagram, relativistic heavy-ion reactions, and to many other phenomena, now abound. 
The challenge for the future is to further understand the successful applications, the failures 
where they occur, the full role of hadronic loops, and the deeper connection to QCD. 


^^There are, however, some nontrivial aspects to these generalizations from a nonrelativistic to a 
relativistic description. For example, the dynamical nature of the mesons introduces retardation 
effects that have not yet been adequately studied. 
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